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Hypo-Elasticity and Elasticity 
BARRY BERNSTEIN 


Communicated by C. TRUESDELL 


1. Introduction 


When TRUESDELL [1955, 6] introduced the equations of hypo-elasticity, he 
stated that it was his intention to find a new concept of elastic behavior, mutually 
exclusive with the theory of finite strain, which reduces to the classical linearized 
theory of elasticity under the assumptions appropriate to the latter. However, 
he remarked that Nor had shown him privately, and published simultaneously 
[1955, 2] a proof that every isotropic Cauchy-elastic material for which the stress- 
strain relations are invertible is hypo-elastic. 

A question then arises whether a given hypo-elastic material is elastic. And 
this inquiry is the theme of our discussion below. In seeking to answer this 
question I was impeded by a difficulty: The equations of hypo-elasticity in 
themselves were not sufficient to define a material well enough to enable me to 
proceed with my inquiry. This hurdle was overcome by the means described 
in §3 below, in which are set down the stipulations which, in addition to the 
equations of hypo-elasticity, I feel are required in general to define a hypo- 
elastic material. 

In the second section we introduce our notation for discussing motions of 
materials and define elasticity in the sense of CAUCHY and of GREEN. In the 
fourth and fifth sections we derive criteria, based on the equations of hypo- 
elasticity, for a hypo-elastic material to be elastic in the senses of CAUCHY and 
of GREEN. These criteria are obtained by means of seeking integrability con- 
ditions. 

The sixth section is devoted to showing that unless there exist motions of 
a given hypo-elastic material for which the initial and final configurations are 
the same and the work done by the stresses is negative, the given material must 
necessarily be elastic in the sense of GREEN. The existence of a result of this 
type was suggested by the work of CaprioLt [1955, 1] and of BERNSTEIN & 
ERICKSEN [1958, 1]. 

Inthe final section, § 7, we prove that a material which is both hypo-elastic 
and Cauchy-elastic is isotropic with respect to a given reference configuration 
if and only if the stress is a uniform hydrostatic pressure for that reference 
configuration. We then indicate that this result implies that there exist an- 
isotropic Green-elastic materials which are not hypo-elastic. 

A summary of our results is given by BERNSTEIN [1960]. 

Arch, Rational Mech. Anal., Vol. 6 7 


90 BARRY BERNSTEIN: 


2. The constitutive equations of elasticity 
We follow generally the notation of TRUESDELL [1952], letting Xe AA 2s 
denote material or Lagrangean coordinates, x', 1=1, 2,3, stand for spatial or 
Eulerian coordinates, and ¢ for time. A (kinematically possible) motion is defined 
by a continuous, piecewise smooth mapping 


x = *'(X4, 2), (2.1) 


where carats, bars, efc., over a letter indicate that it stands for a function, a 
distinction which will be made when the context requires it. Marks appearing 
under or to the right of a kernel letter will be taken to be part of the kernel 
letter. Henceforth all spatial coordinates will be considered cartesian. Compo- 
nents of tensors in spatial coordinate systems will be written with lower case 
Latin kernel letters and lower case Latin indices, while those in material coordinate 
systems will be written with upper case Latin kernel letters and upper case 
Latin indices. Since all spatial coordinate systems are to be cartesian, cartesian 
tensor notation will be used in these systems, and lower case Latin indices may 
be transferred between the upper and lower positions at will. A material point 
will be denoted by X, Xx, etc., when no reference is made to coordinates. 


The derivatives of (2.1) with respect to the X%, called the displacement 
gradients, will be written 


N= ae (2.2) 
It is assumed that (2.1) may be solved for the X“, yielding 
KAKA At), (2.3) 
which leads to the definition of the inverse displacement gradients, X7 by 
$A 
5 ee 
ax?’ 24) 
and the subsequent relations 
Kevg seein enti Xf pian aX det |x/,|. (2.5) 
By differentiation of (2.5) we get 
Bue! + ia|> ead A ax} A WL 


When £',, %',, etc., denotes an assignment of the displacement gradients as 


functions of given variables, X4, X¢ etc. will respectively denote the corresponding 
inverse displacement gradients determined by (2.5) as functions of the same 
variables. 

The symbols 6; and 6% appearing in (2.5) are the usual Kronecker deltas. 
We also adopt the Kronecker symbols 6‘, and 64, which take the value unity 
when A and 7 are numerically equal, and take the value zero otherwise *. 


* See [1954], footnote on page 14. 
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The mass density, @, is given from the principle of conservation of mass by 
Qo det |x44| =, (2.7) 
where 9 depends only on X and in fact is the density for any value of x such 
that det |x!,|=14. 
If at each point of a smooth manifold 3 of material points* the stress, ¢’’, 
is given as a function of the displacement gradients, 
fi = fi (od, X¥), (2.8) 


we say that we have a representation of a material which is elastic in the sense 
of CaucHy, or Cauchy-elastic. These representations form equivalence classes 
if we say that, for fixed IN, two representations 


{it — £7 (a, xo )S and: t= £7 (a, Ga) 


are equivalent if there is a transformation of material coordinates 


ye XX?) (2.9) 
such that bs z 
1 (xq, X?) = #7 (xh, X*) (2.10) 
where 
bir ax sala axa 
ax® ax 


An equivalence class of representations of Cauchy-elastic materials will be called 
by us a Cauchy-elastic material. Such a material is well defined by one of its 
representations. 

The possible relations (2.8) are restricted by the requirement that material 
properties be invariant under rigid motions. This requirement has been given 
the name of the principle of zsotropy of space, or objectivity [1955, 2] [1958, 2]. 
It has also been formulated as invariance under Euclidean transformations, and 
it has been shown [1958, 3] that (2.8) will be consistent with this principle if 
and only if it takes the form 


#7 = [det | xh, |] 24, xy T4 (CEE, x4), (2.11) 
where 
CEE zx XP0", (2.42) 
and T4® are scalar functions with respect to time-dependent orthogonal trans- 
formations of spatial coordinates. 
In case (2.8) has the form 


(2.13) 


where 


is, at each X, a scalar function of the displacement gradients, then we say that 
the material defined by IN and (2.13) is elastic in the sense of GREEN, or Green- 


* MM is the same as the body Z of [1958, 2]. 
7% 
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elastic. (If (2.8) has the form (2.13) for one representation, it has that form for 


all equivalent representations.) 
The principle of isotropy of space is fulfilled, for materials elastic in the sense 


of GREEN, if and only if Sis a function of the CX” and X [1955, 2]. 
3. Hypo-elastic materials 
The constitutive equations of hypo-elasticity may be written [1955, 2, 4, 5] 
f= ty, jn t+ tyr Opn + Asjniltyg) dy, (3.1) 


where, if a quantity / is expressed as a function of X and t#, f/=/(X, ¢), then 
} =0 flat is called the material derivative of f; where 


jn @ (Up Up, Ae d= & (2+ Uz, ae (3.2) 
v; denoting the velocity, which is defined, using (2.1), by 
ax’ 
a ot? (3.3) 


and where Aetna) is an invariant or isotropic tensor function of the stress 


a 


and possesses the symmetries Aj jar = Ayia = 4; -siz- Using (3.3), (2.1), (2.2) and 
(2.4), we may write the velocity gradients as 


te XA (3.4) 
Substitution of (3.4) into (3.1), using (3.2), yields 
bey = Beyarlly) XP, (3.5) 
where a Va 
Bignr = ¥ [bi Opn + ty Ose — tie Oj1— ty Oi2] + Aijar- (3.6) 


In the spirit of differential geometry, we shall assume that es jr and hence 


Ba nv 1S differentiable as many times as required, and that we have enough 
continuity conditions to insure existence and uniqueness of solutions to the 
initial-value problem for (3.1) or (3.6) in which the velocity gradients, or the 
displacement gradients, are assigned as continuous piecewise smooth functions 
of time, stress is assigned initially, and it is required to solve for stress as a 
function of time at each X. Our theorems are local unless obviously global. 

We do not feel that equations (3.1) or (3.5) suffice to define a hypo-elastic 
material. The following example will bring out the difficulty which we have 
in mind. 


Consider the particular case of (3.1) given by 


t 
by j= bp @jnt by, Oi, — 0; 2" dy, 
1.€., 3 


t 
ee p 
ijn — 6;; 2 dy. 


3 
One set of solutions to (3.7) is given by 


t;=—x00,;, (3.8) 


Hypo-Elasticity and Elasticity 5] 


where x is an arbitrary constant. Hence if initially the stress is a hydrostatic 
pressure with a given assigned value, x may be adjusted so that the right-hand 
side of (3.8) takes that value initially and, by uniqueness, (3.8) must be the 
solution. It follows that if ever the stress is a hydrostatic pressure, it is always 
a hydrostatic pressure, and if the stress is once not a hydrostatic pressure, it 
may never become’a hydrostatic pressure during any motion. Thus it would 
seem that (3.7) entails at least two essentially different types of materials, one 
for which the shear must vanish and one for which the shear cannot vanish. 
And we shall see (in the closing comments of §§ 4 and 5) that this difference 
is even more profound. 
What is clear from the foregoing discussion is that in some way the possible 
initial conditions must be specified. And we now proceed to formalize this idea. 
Let I be a smooth manifold of material points as in § 2. A function £%, (X?) 

over J, 7.e., an assignment of displacement gradients over IM, subject to the 
compatibility conditions 

OR ae 

a a 6.9) 

Ox Ox 
will be called a configuration (of Wt). An assignment of both a value of stress 
at each material point in § and a configuration will be called a stress-configuration 
pair and will be denoted by 


{t;;, 4}. (3.10) 
Now (3.10) is a function over J. Its value at one point X of Yt will be denoted by 
[tij, 4] (X). (3.41) 

The symbol 
[t,;, x4] (3.42) 


will denote a possible set of values of stress and displacement gradients at a 
material point. It will be called a stress-displacement gradient pair. The ex- 
pression (3.11) is necessarily associated with a given stress-configuration pair, 
while (3.12) is not. 


Given two stress-configuration pairs 
eat 2 2 
k k 
{t,;,x4} and {t;;, x4}, 


the following statement may or may not be true: There exist continuous piece- 
wise smooth functions 


lee gy ASG es eit See 2h 
satisfying (3.5) and (3.9) such that | 


he 2 ae 
[ij X), 4 (4, X)] = [Fj %a] (KX), = 1,2 
for all X € Me. 
If this statement is true, we shall say that the two stress-configuration pairs 


are equivalent 
| {ij AB} ~ flap BB, 
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and if not, then the stress-configuration pairs are not equivalent. It is readily 
verified that this is an equivalence relation. And an equivalence class of stress- 
configuration pairs will be called a stress-configuration class. A stress-configu- 
ration class for a hypo-elastic material is, then, the class of all stress-configuration 
pairs which may be reached from a given stress-configuration pair. Further, 
these stress-configuration classes themselves may be separated into equivalence 
classes as follows: Two stress-configuration classes are equivalent if one may 
be obtained from the other by a transformation of material coordinates. 


In other words, given two stress-configuration classes, say Jj and 13, J; is 
equivalent to Jj if there exists a coordinate transformation 


X4 = X4 (x4) (3.13) 
such that if 
{t;;; tatehy, 


{ti;, xg} EL, 


then 


and conversely, where 


From equation (3.5) it is readily verified that 


ts; — By; ei (tya) Xie t3— Bzjniltpg) Xf ak, (3.14) 


where 


PMG ELS 
ie ax! ax ars 

Thus if any stress-configuration pair from J] may be transformed into a stress- 
configuration pair from J} by a transformation (3.12), it is immediate from 
(3.13) that every stress-configuration pair in J} is a transform under (3.12) of 
a stress-configuration pair in J,, and conversely. Hence, necessary and sufficient 
that Jj and J; be equivalent is that one stress-configuration pair in J] may be 
transformed into some stress-configuration pair in J} by a transformation of 
material coordinates. 


By a hypo-elastic material we shall mean the assignment of a set of equations 
(3.1) or (3.5) and a corresponding equivalence class of stress-configuration classes. 
By a representation of a hypo-elastic material we shall mean an assignment of a 
set of equations (3.1) or (3.5) and a corresponding stress-configuration class. 


. It may be noted that in a representation of a hypo-elastic material which 
is also Cauchy-elastic there is only one value of stress for each value of the 
displacement gradients at each X, and hence a representation of such a material 
consists of an assignment of stress as a function of x, and X. Also the equi- 
valence classes of stress-configuration classes must be those obtained by con- 
sidering classes of equivalent representations of an elastic material. Hence our 
definitions of a representation and a material are consistent for the case of a 
material which is both Cauchy-elastic and hypo-elastic. 
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4. Elasticity in the sense of Cauchy 


Suppose that we are given a set of hypo-elastic equations (3.1) or (3.5), and 
we ask whether any hypo-elastic material corresponding to them is elastic in 
the sense of CaucHy. If (2.8) defines such a material, then by substitution of 


(2.8) into (3.5) we see that me must satisfy 


for arbitrary motions. But since x4, and «4 are in general independent, (2.8) 
must satisfy (4.1) for arbitrary «4 and hence must satisfy 


BSB INXS, (4.2) 


Conversely, if we have a solution (2.8) of (4.2), then (4.1) and hence (3.5) is 
satisfied by this solution and, because of uniqueness of solution of the initial 
value problem for (3.5) as stated in § 3, the stress-configuration class determined 
by any single stress-configuration pair satisfying (2.8) must be identical with 
the totality of stress-configuration pairs satisfying (2.8). And thus this stress- 
configuration class, together with (3.5), determines a representation of a material 
which is both hypo-elastic and Cauchy-elastic. Thus our question reduces to 
that of whether (4.2) admits a solution. We therefore apply the standard method 
of answering this question as outlined in textbooks on differential geometry *. 


The first integrability condition for (4.2) is obtained by writing 


Ci; Ct;; 
(Selee ee: 
axe 0x4 Ox 4 axe (4 3) 
A A A a a 1G, 
= 5 Ai ay Beal: ipl =e ee 
Oly, OXF Ox Oly, O%%4 OX 4 


By substituting (4.2) for the derivatives of t,, and using (2.6),, we get a set of 
relations involving X# and ty, From these X} may be eliminated using (2.5)o, 
so that the first integrability conditions for (4.2) take the form 


oe rspq see BiB; ng Opi t Bijp1 Ong = 9- (4.4) 

For ¢,; to satisfy (4.2) it is necessary that its components be related in such 

a way as to satisfy (4.4). If (4.4) is not satisfied identically in 4;; and does not 
consist of more than six independent relations, we may differentiate (4.4) with 
respect to x4,, substitute the right-hand side of (4.2) for the derivatives of #;,, 
and again eliminate X? by the use of (2.5), to obtain the second integrability 
conditions as relations among the 4,;. If the second integrability conditions 
introduce new relations not contained among the first integrability conditions, 
and if the total number of independent relations among the first and second 
integrability conditions does not exceed six, we obtain the third integrability 
conditions by applying to the second integrability conditions the same process 


* See, for example, § 23 of [1947]. 
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that we applied to the first. And we continue this process till either we obtain 
more than six independent relations, or at some stage no new relations are 
introduced and the total number of independent relations is y<6. In the former 
case there is no solution of (4.2). In the latter case there are solutions for which 
6 —r components of ¢;; may be assigned arbitrarily for one value of x*, at each X. 
And for each such assignment the solution is unique. 

In case (4.4) consists of more than six independent relations there is, accord- 
ing to the above discussion, no solution of (4.2). On the other hand if and only 
if (4.4) is satisfied identically in ¢,; is there a solution of (4.2) for any arbitrary 
assignment of ¢;; for one value of x" at each X. Thus, because of the uniqueness 
of solutions to the initial value problem for (3.5), we have the following: 

Theorem I. Necessary and sufficient that every hypo-elastic material cor- 
responding to (3.1) or (3.5) be Cauchy-elastic is that (4.4) be satisfied identically 
mM t;;. 

As an example to illustrate the foregoing discussion, we consider the equations 
(3.7), for which the totality of independent integrability conditions may be 
written 
= ik ms 


“i= aq? 


(4.5) 


showing that a hypo-elastic material FB) i to (3.7) is elastic in the 
sense of CAucuy if and only if the stress is a hydrostatic pressure. 


5. Elasticity in the sense of Green 


Given a material which is both hypo-elastic and Cauchy-elastic, we now 
ask for a criterion that it be elastic in the sense of GREEN. We ask, in other 


words, if there exists a function D (x4, , X*) satisfying (2.13). We rewrite (2.13) as 
as Q A 
Sais eee (5.1) 
x4 Q 
Since we are assuming that (2.8) is satisfied, a necessary and sufficient condition 
for the existence of a solution 2 to (5.1) is 


as as a 


20 7A | é e@ L 
= =o t;, XF = = 
On BL MO ae Oe net i axi, | 9 big Pe 
From (2.7) and (2.6) we get 3 
—, £0 — 2 x4. (5:2) 
Ox 4 0 a 


Since we are also assuming that (3.5) is satisfied, we get from (5. 
1 
(3.5) and (2.6), 8 (5.1), using (5.2), 


Bean wigs Xe janie XG + 4,; Xk XP — 4, X4¢ XP 2 KA Xj tj; + XPXA th; = 0, 


which, using (2.5), may be written 


Be Boe See jet tp Oj; — tp Oi + ty; 6); = 0. (Se) 
Using (3.6), we see that (5.3) may also be written 


Aj jpi tb; O41 = Ap ijit ty10;;: (5.4) 
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We state our results now as two theorems: 


Theorem II. Necessary and sufficient that a given hypo-elastic material which 
ts Cauchy-elastic be also Green-elastic is that (5.4) be satisfied for all values of 
stress that occur. 


Theorem III. Necessary and sufficient that all hypo-elastic materials cor- 
responding to (3.1) or (3.5) be elastic im the sense of Green is that (4.4) and (5.4) 
be satisfied identically 1n t,;. 

We see that (3.7) satisfies (5.4) for Cauchy-elastic solutions, 7.e., when (4.5) 
is satisfied. Hence Theorem II implies that hypo-elastic materials correspond- 
ing to (3.7) are elastic in the sense of GREEN if ever the stress is a hydrostatic 
pressure, which is to be expected since they are perfect fluids in that case. If 
ever, on the other hand, the stress is not a hydrostatic pressure, then the material 
is not elastic even in the sense of CAucHy. Thus the stress-configuration class 
can be very important in determining the nature of the material corresponding 
to a given set of hypo-elastic equations. 


6. An energetic condition for elasticity in the sense of Green 


The work done by the stresses* in any material during a motion occuring 
in the time interval 4; </<#, is given by 


by 
Lfovsarar hips 20 dVy at, (6.1) 


t Vo 


where V is the volume of material at time ¢ and M) is the volume for some fixed 
reference configuration, for which the density is @). In this section we shall 
prove the following: 


Theorem IV. I/f for a given hypo-elastic material the work done by the stresses 
is non-negative for all (kinematically possible) motions for which the initial and 
final configurations are the same, then this material is elastic in the sense of Green. 

Before proving the theorem, we shall make some observations and establish 
some lemmas. 

Consider a given representation of a hypo-elastic material and a given material 
point x, Let {’t;;,'x4} be an element of the corresponding stress-configuration 


class. Suppose that |t;;, x4] is a stress-displacement gradient pair for which 


there exist piecewise smooth functions a (t), xt (t), 4) StSty, satisfying (3.5) and 
such that rs 

[#5 (ta), 24 (4) = ['6;.'"4] (X), 

and " 


On OEaane oie 


Then we may find an element {'"t;;, "x x". of the stress-configuration class cor- 


responding to our given representation Sieh takes the value [""4,;, xk] at Xx, 
as we shall now demonstrate. 


x § 27 of [1952]. 
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Let '%*,(X) be the values of the displacement gradients corresponding to 
{t,,, x4}, and let Xe be a set of numbers defined by 


xe th (X) = 68, 


i.e., the inverse displacement gradients at x corresponding to {’t;;, 'x4}. Then 
define the function 
HA (X,1) =3h() XP (X), 4SiSh (6.2) 


We see that (6.2) satisfies (3.9), since ‘t!,(X) does by assumption; that 


and that 


for all X. Hence (6.2) defines a kinematically possible motion in which the 


values of the displacement gradients are prescribed at x. And to find the stress- 


configuration pair {’t,;, ’*4} we have but to put (6.2) into (3.5), take 1 Tae 
as initial conditions at {=t,, solve (3.5) for 4: <t<t#,, and determine the stress- 
configuration pair at t=?,. 

For a given representation of a hypo-elastic material, we may speak of the 
allowable stress-displacement gradient pairs at each material point as the possible 
values that the stress-configuration pairs may take there. And, from the fore- 
going discussion, these allowable stress-displacement gradient pairs will be pre- 
cisely all those which may be reached from any one of them along a solution 
of (3.5). 

Moreover, an examination of (6.2) will show that an assignment of Ey (¢) at 
x such that %4 (,) = x4 (ts) may be imbedded in a motion of the entire set of ma- 


terial points for which the initial and final configurations coincide. Thus if it were 


possible to find a solution [i;;(é), 4 (0)], 4<tSty, of (3.5), where [f,, (4), #4 (4)] 
is an allowable stress-displacement gradient pair at some given material point 
x for some given representation of a hypo-elastic material, such that “w <0, 


ea BS 
bh 


then for some material volume containing x it would, by (6.1) and our foregoing 


where 


6S 
o 


dt, (6.3) 


FDW 


remarks, be possible to find a motion for which, though the initial and final 
configurations are the same, the work done by the stresses is negative. Hence 
we may make the following 


Remark 1. Suppose given a representation of a hypo-elastic material for which 
the hypotheses of Theorem 1V hold, henceforth called a representation 1V. Then 
the quantity ~w defined in (6.3), and henceforth called the work density, must be 
non-negative for any set of continuous, piecewise smooth functions es xt (2), 
4) StSty, satisfying (3.5) for which |t,;(t), xt (t,) | ¢s an allowable stress-displacement 
gradient pair at any given material point X such that X*, (4) =x (to). 
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And this shall be our starting point to establish Theorem IV. 
Note that (6.3) may be written 


dx’, 


by 
ae x4 Saar at, (6.4) 


the total anne replacing the material derivative because, since X is fixed, 
x', becomes a function of t alone. We now proceed to our lemmas. 


Lemma 1. /} [t,t Bs (t)| satisfy (3.5), ees and if f(t) ts a smooth 


function of t, ty St<t,, dt/dt +0 with, say, i(t 1) = by, E(t) = te, then t, AU t), i (t) 
defined by 5 


St 
Ss 
a, 

eK 
~—” 

| 
oy) 
Ss 
os 
sSeet 
Rx 
bs 
——s 
~ 
™—" 
I 
&) 
no 
— 
~ 
Sa 
xo 
os 
* 
ae 
~ 


satisfy (3.5) for min (4, t) St < max (4, te). Also the corresponding values of the 
work density, namely 


ty 
Pale Bl a at, 
and (6.5) 


ve fost eee 


are related by 


oY 


Proof. This lemma is readily verified by substitution in (3.5) and (6.5). g.e.d. 


It follows from Lemma 1 that, without loss of generality, we may assume 
that 4;=0 and ¢,=1 where convenient. We proceed to 


Lemma 2. Consider a representation IV and any material povnt x. et 


[z;,;(), Ki (t)|, 4StSt,, be any solution of (3.5) with [4 (4), A (t,) | an lsouBle 


stress-displacement gradient pair, and % (t,) =2 (t). Then “w=0, where “w is 
the corresponding value of the work density. 


St y 
SS 
| aoe 
xo 
™— 
Re 
ho 
a 
xo 
——" 
| Sees | 


Proof. Using Lemma 1 we may obtain a second solution of (3.5), [4; 
hSts t,, obeying the hypotheses of Lemma 2, by writing 
ys es 


Let *w be the value of the work density corresponding to this second solution. 
Then by Lemma 1 w= — *w. But, because of Remark 1,20 and *w 20. Hence 
w=—*o=—0. g.e:d. 

Lemma 3. Consider a representation IV. Let [tig, ¥4] and [tigs 2A] be two 


allowable stress-displacement gradient pairs at a given material point. Suppose 
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that AGE x (t) | and [t,,(t), x',(t)] are two solutions of (3.5), 4StSth, us CET, 
a = = = k ZR 

such that [fs (h), Fh (6) ]=Ltags hu] = [ees(ta), 24 (ta)] Ba (4) = 2h = Fh (ta). Let the 

corresponding values of the work density be w and “w respectively. Then w=~w. 

Proof. By Lemma1 we may, without loss of generality, take 4 =1=0 and 


ty=T,—1. Then we define another solution of (3.5), ed ee PES RSS. 
as follows: 


That this is a solution in each of the intervals O<’t<1 and {S782 may be 
verified immediately with the use of Lemma1. Also it is continuous at be 
because, by the hypotheses of our lemma and the definition of this solution, 
[4;;(0), #4 (0)] = [4;;(0), #4 (0)] =[4, (1), ¥4 (1)]. Hence, since ¥4 (0) = #4 (1) = 2% 
i* (41) =X4 (2), Lemma 2 yields 


ax’ 
“w= f det |x] 4; X/ EB (6.6) 
0 
But from Lemma 1 
il 
ose, Ue 
—k A A 1 as 
f det [x tip Ap sae Aion 
i (6.7) 
det ele ax- a d't="w 
jf 2 [Xp ESA = she — 


and hence the lemma. g.e.d. 


Remark 2. Suppose that we consider a given stress-configuration pair, {ti 7 xa} 
corresponding to a given representation IV. Then, by Lemma}, we may define a 
function E(x, X) as follows: for each fixed X, 3S ('x%, X) is the value of the work 
density for any solution of (3.5) for which [t, > #4 | (X) constitute the initial conditions 
and for which the final values of the displacement gradients are 'x',. 

We cannot differentiate }’ with respect to x', by ordinary means since Dy 


is defined by an integral whose integrand has not been shown to be a function 
of x4 and X. Nevertheless, we shall now prove 


Lemma 4. The partial derivatives of 2) (x4, X) with respect to the displacement 
gradients exist and are given by 


~ 


2A = det |xh| 4, X4, (6.8) 
Ox, 


where, for given x", t;; has any value such that [t;;, x4] is an allowable stress- 
displacement gradient pair at X. 
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Proof. Let tae | be any allowable stress-displacement gradient pair at X. 
Then there exists a solution of (3.5) [é;;(t), #4 (t)], OS¢S1, such that 


By definition 


gl 3 
ne Ce ee A 
Exh X) = f det |x|, Xf FA at. 
0 


Let us extend the function %% () beyond ¢=14 in any one of the following eighteen 
ways; write 
#4 (t) =#4+ e(¢—1) 6904, 15¢S1+4h, (6.9) 


where # and L each take one of the values 1, 2, 3, ¢ takes either the values +1 
or —1, and h is some positive real number. Let 4,,(¢) be the corresponding 
solution of (3.5). Then 


rae De F ae if] a = 

2 (x4 + en 6; 04, X) — 2i(x4, X) = f det |#% (t)| t,; (6) X7 (0) edt 
1 ia i (6.10) 

= det | #5 (t)| t,;(t) XF(t) e(n — 1), 


where 1<7<1-+h, and ft is some value between 1 and », the last inequality 
following from the law of the mean for integrals. Using (6.9), (6.10) becomes 


~ 


zt 
fae) 


= det |#%(z)|t,;(t) F(t),  1<t<n. (6.14) 


Since the limit of the right-hand side of (6.10) exists as 7->1, the limit of the 
left-hand side must also exist. Since %%4 (7) = 4 for 721 unless t= and A=L, 
the limit of the left-hand side of (6.10) as 7->1 represents the right or left-sided 
partial derivative of » with respect to xf according to whether ¢ is +1 or —1. 
But both of these limits are equal and, in fact, are given by the right-hand side 
of (6.11) evaluated at t=1, since solutions to (3.5) must be continuous in ¢. 
Hence the partial derivatives of » with respect to the displacement gradients 
exist and are given by (6.8), where ¢;; is determined as in the hypotheses of 
Lemma 4. q.e.d. 

With the help of our lemmas, Theorem IV follows immediately. From 
Lemma 4, in particular (6.8), and from (2.5),, a stress which will make [f,;;, xt] 
an allowable stress-displacement gradient pair at any given X, for a given repre- 
sentation IV, must satisfy 

OCA PRe 5; 02 (ey. X) 


lag SR teria ee, 


(6.12) 


and thus there can be but one such stress. Moreover, comparison of (6.11) with 
(2.13) shows that a material obeying the hypotheses of Theorem IV must be 
elastic in the sense of GREEN. Hence Theorem IV is established. 
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7. Isotropy 
Consider a representation of a Cauchy-elastic material in which the material 
coordinates are one possible set of values of the spatial coordinates and are 
cartesian. Let O04 and 04, A, B=1, 2,3, be a set of constants such that 


4p OC OF = 0cp- 


Then 
XA Oe Xe (7.1) 
denotes an orthogonal transformation of material points so that 
ax" i OB eS: Sa BOA 
ly “_ = 45,0 % AS ee . Lee 
XA 3’x4 BVA ax B (7.2) 
If in (2.8) we have Sai. 
ty; (x4, 'X”) = t,;(%4,X"), (7.3) 


when X4 and 'X“ are related by (7.1), we say that the given Cauchy-elastic 
material is isotropic with respect to the given representation. 


It is well known [1952] [1955, 3] that the above definition of isotropy requires 
that the stress 4, be a uniform hydrostatic pressure when x',=6%. But the 
converse of this statement is not true, as can be verified by considering a Green- 
elastic material for which 


Ree [(O_: 4 Xe — O4p) HA, (7.4) 


where H1!1=1, H?2—2, H33—1 and H4?=0 if A+B. That (7.4) is consistent 
with the principle of isotropy of space may be checked using the last paragraph 
i. 8-25 
Using (2.13), the stress for (7.4) becomes 
fis a (Sgr +4 4g — O4n) HA? of HE™ ah, (7.5) 
0 


Consider the following two special cases of (2.1); 


, 


be Mae I, ae IA, gh NE (case 1) 
and 
El ae feo, CRY tae (case 2) 


which differ by the orthogonal transformation 


PCS, CD, Cw, CD, Ce, CF 
We obtain from (7.5) 


pia 6, PAS 8 6 On eo (case 1) 
#2=12, £8=6, B8=3, fFia0, 7 HF (case 2) 


so that the Green-elastic material given by (7.4) is not isotropic with respect 
to the given representation although from (7.5) it is readily seen that the stress 
vanishes, and hence is certainly a hydrostatic pressure, when x4 = 8. 
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We shall now establish the following: 


Theorem V. Given a representation of a Cauchy-elastic material which is also 
hypo-elastic, this material is also isotropic with respect to this representation if and 
only if the stress is a uniform hydrostatic pressure when x',=6',. 


Proof. Necessity has been established elsewhere, as stated above. We have 
to prove sufficiency. Assume, then, that the stress is a uniform hydrostatic 
pressure when x)= 064. 

Let (2.8) correspond to our representation and define 


t, (24, x) = t,,('"4, X, (7.6) 
using (7.1) and (7.2). Then, by (7.2), 


Bay ) = G (Ont On x LON, (7.7) 
and by substitution into (4.2) we find that iss is a solution. But when x/;=064 
we have t;;= — f0,;;, where p is independent of X. Hence, using (2.11) and (2.12), 
we see that 


— p 51 =17(d5, X°) = by Ob T4” (Of Oy 0°, X°) = bi Oy T*? (*", X9), (7.8) 
from which a 
Pe alee ae? face TO A), V2) 


the last equality following since # is independent of X. Thus, putting x/,= 6, 
into (7.6) and using (7.7), (7.8) and (7.9), we obtain 


i7 (6%, X°) = 14 (680%, "X°) = 52 OF 66 08 T4” (af! OF BN OF 6?8, 'X°) 
= d,0% 66 0§ T48 (5X*,'X°) = — p dz OF 0503 64” 
= — p bh & FS = — p d= Hd, X%. 


Since the two solutions of (4.2), i and fi, coincide for one value of x‘, at each X, 
they must, because of uniqueness of solutions to (4.2), coincide for all values 
of x), at each X, 1.¢., 


i (‘ach X°) =F (xh, X°) = Bia, X9), 


and hence the theorem is established. 

We remark, then, that the material (7.4) is elastic but not hypo-elastic, since 
if it were hypo-elastic it would have to be isotropic because of Theorem V, and 
it is not. Thus hypo-elasticity does not include all of elasticity. 
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Inversion of Lambda-Matrices and Application 
to the Theory of Linear Vibrations 


P, LANCASTER 


Communicated by A. M. OSTROWSKI 


1. Introduction 
We consider the problem of inverting the matrix 


D® (A) =(A)A+A4,N14---+4),A+4)), (1) 


where Ay, A,,..., A, are (n Xm) matrices of real elements, and J is a scalar. 
The special case /=1 is comparatively well known, for (assuming A, to be non- 
singular) it can be reduced to the problem of finding the inverse matrix of 
(IA— A), known as the resolvent of A (see [1]). The resolvent of A can be 
viewed as a series of partial fractions with denominators (A—4,)/, where A; 
are the latent roots of A. In the same way the inverse of matrix (1) may be 
expanded as a series of partial fractions involving its latent roots. This was 
done by FRAZER, DuNcAN & COLLAR [2]. We are to develop a rather different 
form for the inverse of (1) which has the advantage of being developed entirely 
in matrix form, and which is unaffected by the presence of a latent root of 
multiplicity « (>1) provided (1) has degeneracy « when evaluated at this 
latent root. 

The method consists in first finding the inverse of D@)(A), in which case 
the matrix elements are linear in 4; and then for />1 we transform the problem 
to that of inverting a matrix of order /m whose elements are again linear in A. 
The application to the theory of linear vibrations arises when /= 2. 

We shall assume throughout this paper that the leading matrix, Ag, is non- 
singular. Consider the set of homogeneous equations 


(Aja +4, 2*4+---+4))q=0, (2) 
and the set obtained from the transpose of (1), 
(Aga + AA 4+.---4+ A) r=0, (3) 


where q and r are (m X1) column vectors. The latent roots of matrix (1) are 

the values of 4 for which its determinant is zero and for which (2) and (3) have 

non-trivial solutions, known as the latent vectors. If we count each root accord- 

ing to its multiplicity, there will be /m latent roots, 4;, in all, each with associated 

latent vectors gq; and r;. We shall always assume that if a root of multiplicity 

a occurs (1Sa<n), 4,, say, then D”(A,,) has degeneracy «. That is, there are 
Arch. Rational Mech. Anal., Vol. 6 8 
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a linearly independent solutions of equations (2) and (3) when / is replaced by 
1,,*. We arrange the latent roots and vectors in some fixed order. Define the 
(n xln) modal matrices, Q and R, as the matrices whose columns are the latent 
vectors of (2) and (3) respectively. We also define the (ln xln) diagonal matrix A, 
so that its 7 diagonal element is 4;, the 7” Jatent root, and note that the elements 
of Q, R, and A may be complex numbers. 


2. Inversion of lambda-matrices 
The following theorem is not new in principle, but the formulation and proof 
given below are the most appropriate for the present work. 
Theorem 1. If A and C are real (nxn) matrices, of which A is non-singular, 
then the latent vectors oj the equations 


(AA+C)q=0, and (A'A+C')r=0 (4) 

can be normalised in such a way that 
R AO =1)and ~R CO=—A, (5) 
in which case (Ad+C)4=O0(A— AAR, i# ALD, 6) 


and provided that for any latent root i,, of multiplicity a (1Sa<n), the matrix 
(AA,,+C) has degeneracy a. 


Proof. We notice that in this case Q, R, and A are all ( xm) square matrices. 
It follows from our definitions that 


ACA COO 


ad AURAL GOR So: 


Premultiply these equations by Rk’ and Q’ respectively, and after transposing 
the second we obtain 


(R'AQ)A+R'CQ=0 and A(R’AQ)+R’CQ=0. (7) 


aa (R'AQ)A=A(R’AQ). (8) 


If no two of the latent roots are equal, this can only be true if (R’AQ) is a 
diagonal matrix. If there is a single root of multiplicity « (1<a<m) among 
the latent roots, we order the m latent roots so that the « equal roots take the 
first « positions. From (8) we then deduce that the leading (« x«) minor of 
(R’ A Q) is the only part of that matrix which is not in diagonal form. However, 


for arbitrary constants f,, By,...,B, and y,,72,...,Y,%, we have 
‘ (AA, + C) (B19 + BoQe+-:- + Bade) =0, 
an , , 
(APA + C') (11. + 29a +++: +a) =0. 


* The problem which stimulated this work was that of linear mechanical vibrations. 
The two strong conditions I have made — that A, is non-singular, and that DM (1,,) 
is of degeneracy « — will invariably be satisfied in these applications. The inherent 


simplicity of this important case seems to be lost in the treatment of FRAZER 
DuNcAN & COLLAR. ‘ 
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Thus any linear combination of the first « latent vectors is also a latent vector 
with latent root 4,. By hypothesis the matrix (4/,+C) has degeneracy «, and 
this implies that the vectors q,, Qs, ..., Y,, are linearly independent, and so are 
T,,1%g,...,%,. We may therefore use the indeterminacy of these vectors to 
complete the biorthogonal system which will reduce (R’ AQ) to diagonal form. 
The extension to cases with more than one multiple root follows immediately. 
We then obtain the first of equations (5) by an appropriate normalisation of 
the latent vectors. The second of equations (5) then follows from (7). 


Using equations (5) we now see that 
R'(AA+C)Q=IA—A. 
If 2+/,;, we may invert both sides of this equation and obtain 
O4(AA+ OC) 4h) 4] 714A 
which implies equation (6) and completes the proof. 


Theorem 2. If A, is non-singular and D(/;) has degeneracy equal to the 
multiplicity of A; for 1=1,2,...,1n, then the latent vectors of equations (2) and 
(3) can be normalised in such a way that 


AD (a)j47= QA (LTA—A)AR’, for 7:=05192; 2.25 (1 — 4), (9) 
and 


#[D(A)}4 = QA (LA — A)4 RB’ + Ag (10) 
provided A= /,. 


Proof. (We assume from this point onwards that /22.) We note that Q 
and R are now (m X/n) matrices, and A is an (Jn x/n) matrix. For any latent 
vector # (of either (2) or (3)) we denote by x the product A*a, A being the 
corresponding latent root. It can easily be verified that (2) is equivalent to the 
following partitioned matrix equation (suggested by Foss [3]) whose elements 
are linear functions of A. 


0 O 0 Ay 
Oy wipers Ay A; 
: At AltA. 
: hes ire ta 
Gan) Ay 
0 A, A, 
A, Ay.A A,_ 
Oneeimet 2 I-1 (11) 
0 0) 0 mae gn ete 
0 EA Ay q 
(/—2) 
A aay Gea ostilla| oe 
* ; an hes 
—A, = Ap 0 gq 
Sov — Ay = ay 3 ee 0 q 
0 0 0 Of ay 


8* 
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Define the two partitioned matrices appearing in this equation as VY and @ 
respectively. Since A, is non-singular, we see immediately that . is non- 
singular, and that @ is non-singular if Ag and A, are non-singular. Further, 
if all the matrices A, are symmetric, then so are and @. Define the (Jn x1n) 
modal matrices associated with 7%, @ and ./’, @’ to be 2 and & respectively; 


where, from our definitions 


O Aa” eA 
Q= OA and #= cae (12) 
Q R 


If A,, is a latent root of (2) with multiplicity « (1<a#<m), then by hypo- 
thesis there are « linearly independent latent vectors q; associated with 4,,. 
It is clear that the corresponding columns of 2 are linearly independent also, 
so that we have « independent latent vectors of (9) corresponding to the latent 
root J,,. That is, (WA,,+@) has degeneracy «. 

The conditions of Theorem 1 are therefore satisfied, and # can be defined 


7 x 
in such a way®* that (MALC)A= DIAL AAS’. (13) 


For convenience, we define (7A+@)1=@ and denote the (” xm) partitions 
of B by B;;. On substituting from (12) into (13) ** we find that 
By=OAG LAS ATA Hig 

OA (TA Ay ee 


Thus B;; depends only on the sumz-+j. That is, there are only (2/—1) distinct 
partitions of @, given by (¢-+7)=2, 3,..., 2l. 
Now by our definition of Z we have 


(JA+¢6)B=TI,,, 


where we use J,,, for the unit matrix of order /m, reserving J for the unit matrix 
of order ». Writing the left-hand side of this relation in partitioned form we get 


0 0 : 0 =e A ah 
0 0 mY Pd wal hyde) AYA 
0 (AgA—A,) (AjA—A,) AGA 
0) : x 
0) —Ay (Ag 4—A)) ; : ; 
—Ay (A, A—A,) (A,A—A)) - i (A;_34—Aj)_») A;_24 
Ags AA A,A Aj24 (4,1 4+4)) 
 15;,] =f), a 


* The choice involved here concerns the definition of columns of 2 and @& 
corresponding to a multiple root. If we take any linear transformation of these 
columns, the form of equations (12) will not be destroyed. 

** In this and subsequent manipulations of partitioned matrices the matrix 


partitions can be treated in much the same way as scalar element 
s. See e.g. 
Duncan & COLar, loc. cit., § 1.7. e.g. FRAZER, 
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We use this equation to evaluate the partitions B;,. The partition of the 
product on the left appearing in the (1, /) position gives 
— A, B,1,+AA, B,,=0, 


and since A, is non-singular, 
- By_1,;=AB,}. (16) 


The partition of the product in the (2, /) position gives 
= Ag Bria (Ag A A;) Bp 752A, A By = 0, 


and using (16), we find that 
( ) By = 2 By). 


Continuing in this way up to the partition in the (/—1,2) position, we find that 
Bian OBS tO ds Dee Ll) (17) 
From the partition in the (/,/) position we obtain finally 
Ay By, A+ A, ByjA+ +++ +A, 2 B41 A+ (Ap A+A) B, =I, 


and using (17), 
(Apa + A,A> +---+ 414+ A) By =. (18) 


From (18) and (14) we get 
B,, = [D(a] *=Q(LA—A)*R’; 
and from (17) and (14), 
B,, =A [DO A)]*= QA (LA A)AR’ 
for 7=1, 2,...,(J—1). These two results establish equation (9). 


To prove equation (10) we consider the partition of the left-hand side of 
(15) in the (1, 1) position. We obtain 


Se ats Anh Bae Ee 
A By y= By11+Ao'- 

But B,,=B,,=/'"B,,, using (14) and (17); so that substituting for By_,, 

from (14), we obtain By = QAMLA— AAR’ + Apt, 


or, 


and (10) follows from equation (18) *. 

Equation (9) with y=/—4 (or even equation (10)) may be valuable in obtain- 
ing [D(A)|+ for large values of A. We shall see below that some interesting 
results are obtained from these equations when A is zero. 


3. Properties of the latent vectors 
Applying Theorem 1 to the matrices , @ again, equation (5) becomes 


R' AQ=1,,, and #CQ=—A. (19) 


* These results can be continued to values of 7,7 for which 7+j7</, but they 
rapidly become more cumbersome and are of doubtful value, so they are not included 
here. 
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These are the orthogonality conditions under which (13) is true. We see from 
equation (7) that if one of these is satisfied then so is the other. We define 


(In x1n) matrices A¥* by 

Av = KAO. 
Now substitute for Z #, from (12) and for ~%, @ from (11), and it is found 
that we can rewrite conditions (19) in the form 


A#t,+(AAhe+ Abed) + (A? Af s+A Af 344 AfsA%)+--: (20) 
(At Ak EAM AR AS AAS A AG A ie ie 


and 


A* — (A A* A) — (A? A? 5 A+A AP AY — --- 


ras 
1. — (AMT AB A+ AR AS AR. + AAS AMY) = — A, ie 


and, as noted above, either one of these implies the other. The complication 
of this normalising condition detracts from the apparent simplicity of equations 
(9) and (10). 

If we pick out the 7 diagonal element from (20), we see that 


W{4,4+24,A;2+°°° +1" A, $q;= 1. 


Thus the normalisation condition can be written in the scalar form; 


p @ 


r=, | 


D(A) =n G=1, for 7=1, 2,i-50e 


Of the sets of /m vectors obtained from (2) and (3) at most can be linearly 
independent. Relations between these vectors are provided by the (2/—1) 
matrix conditions obtained by equating the expressions given in (14) to the 
corresponding expressions obtained by analysis of (15). The statement of 
Theorem 2 contains the most simple (J+ 1) relations obtained in this way. 
However, to be of any value (from this point of view) these conditions should 
be independent of A. We obtain (/+-1) such conditions by considering the limits 
of equations (9) and (10) obtained when A->0. This may lead to difficulties 
if there are any null latent roots, for these will give rise to at least one element 1// 
appearing in the matrix ([A—A)1. 

In spite of this we can prove the following: 


Theorem 3. Under the conditions of Theorem 2 
QA R= Aj’, (22) 
OAS R= 08 for u = 0,4, 2,2. b= 25 (23) 
and 1f in addition A, is non-singular, then 
QATAR =—A;. (24) 
Proof. Putting y=0 in equation (9), we have 


[D9 (4 = OTA AAR” (25) 
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If we now substitute this back into (9) and (10) (with r+0), we get 
MQ(LA—A)1R’=QA(IA—A)7R’ 1Ole t= 1 eens el) 
AOLA—A)*R’ = OAT A =A RR’ + Ay. 

Rearranging these equations, we get 


Q(I# —A’) (IA —A)1R’ =0, fore — 4 2s (t =r) 


and 


ga QU # — A) (1A — A)7AR’ = Ap. 
Th 

us, Q(L ar +4 A+ ..-+A’—1) R’=0, 1Ol Fao mee (be 1g 
and 


OU RA Se + A’) R’ = Aj’. 


We may now put A=0 in these equations and obtain the results (23) and (22) 
respectively. 

If A, is non-singular, then equations (2) and (3) have no solution when 1=0. 
That is, there are no null latent roots, and the matrix A+ exists. We may 
therefore put A=O in equation (25) (remembering the definition of (4)) and 
obtain the result (24). This completes the proof. 


4. Ordinary differential equations with constant coefficients 
Consider now the set of ” simultaneous differential equations 
D9 (DQ) q = (An>B+tAD* ++ +4, 149+4) q=8(), (26) 


where Y denotes the differential operator d/dt, and q, &, are vector functions 
of the independent variable ¢. We see immediately that if A; is a latent root 
of D(A) with associated latent vector q;, then e”’ q; is a solution of (26) when 


E (2) =0. For DY (2) eit qi aes ert DY (A;) q; == (0h, 


Under the assumptions made throughout this paper there exist /m such 
solutions which are independent vector functions of ¢. The complementary 
function for the system of differential equations is therefore 


G0 242; (27) 


where € is a column vector of /m constants to be fixed by the boundary con- 
ditions. A particular integral of equations (26) can be written in the symbolic 


aan q = [D9 (D)]7E (0. 


If the latent vectors are normalised according to (19), (20), or (21), then we 
may apply Theorem 2 with r=0, and the particular integral becomes 


q=QITD—A)*R’'E(i). (28) 
Hence t 
Jail Gut tds &(x).4.vs (29) 


where {, is an arbitrary constant. On combining this with (27), the lower limit 
of integration can be absorbed in the vector ¢, and the complete solution can 
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be written as 


q=Qet{et fe RE() ath. (30) 


If any of the first (/—1) derivatives of the solution are required, we can 
use equation (9) to obtain 


t 
B q(t) = QA et! {e+ f e4* R'E (x) de}, (34) 
for 4=0, 1, 2,0, (0-4): 


5. Linear vibration problems 
These problems often give rise to the equations 
D® (2) q=AoG+419+429=5(). (32) 


In this case the normalising conditions (20) and (21), on which Theorems 2 and 3 
depend can be written 


A¥ + (A Ag + AG A) = Len, (33) 
and 
JE ate EN AS, Wh eee (34) 


where A*—R’A;Q are (2m X2mn) matrices. We then have from Theorem 3; 
QAR'’=Aj', QR’=0, and if A, is non-singular, QA71R’=— <3". A general 
solution of (32) is given by equation (30). 

However, another interesting decomposition of D®)(A) is possible, and this 
gives rise to a rather different particular integral which can be obtained without 
normalising the latent vectors. 

Divide the 2” latent roots into two sets of m roots, and form with these 
the (m x) diagonal matrices A, and A,, any set of multiple roots being contained 
entirely in A, or in A,. Suppose that the corresponding latent vectors form 
the non-singular (7 xm) matrices Q,, Q. and R,, R,. Then in partitioned form, 


A, 0 
Q=[01,2.], R=[R, RJ], and A=| ; i (35) 
Oig Ae 
In equation (32) make the transformation to coordinates p; given by 
q = QP ’ 
and premultiply the resulting equation by Rj. Then putting R;§—T, we have 
(4) 2? + 4,9 +A4,)p=b, (36) 


where A;,=R2A;0,, for i=0,1, or 2. 


The result we are now to obtain (equation (37)) is a decomposition, or 
factorisation, of the matrix operator in brackets in equation (36). The same 
sort of reasoning as that used in Sections 2 and 3 is required, but in order to 
avoid tedious repetitions I shall merely outline the argument. 

If we refer back to Theorem 1, we see that equations similar to (5) can be 
established if we relax the normalisation condition, but replace the right-hand 
side of these equations by unspecified diagonal matrices (although (6) will no 
longer be true). The same is true of equations (19), (20) and (24): 
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If we now substitute from (35) via the definition of A* into (33) and (34), 
supposing that the above normalisation has not been performed, then we obtain 


Ry A, 0O|[R Ri v4 PAG) 
gf 4nl2s Cel | || ee] Aad On) +] op 


= (diagonal matrix), 


| 4ol0. 0.) | 


and 


R gee Otis YE 70 
2 AiO, Qa) — | is | | BI 4o(0, al | 2 i = (diagonal matrix). 

If we pick out the partition in the (2,1) position from each of these equations, 
we obtain 


A, aD St) (A, Ay pai A, At) , 
and 


“a 


A, =AyAy Ait 
Substituting these expressions in (36), we obtain 
[Ap B— (A, Ag+ A, Ay) D+A, A, Ay] p=, 


and the left-hand side factorises to give 


(IZ—A,) A,\ID—A,) p=. (37) 
Thus, 
p=(I2—A)* Ay I D—A,)7E, 
t Lo T (38) 
et [eh An es [eA Goya. 
to os 
On inverting the order of integration the t integral can be evaluated, and in 
terms of the unnormalised modes we obtain the particular integral: 


P= OLD QJ [B ett") — eB C(y) dv, (39) 


where f is the (n xm) matrix with elements «,,/(A;—4;), «;; are the elements 


ij 
of Ag", and A;, A, are the 7" and j diagonal elements of A, and A, respectively. 
If &(¢) is periodic, then it may be possible to represent it by a Fourier series. 
For such a case, and for the practical problem of steady sinusoidal motion 
occurring in resonance tests, it is important to consider the case in which 


E(t) = & ds 


where &, is independent of ¢. (Note that i is now /— 1.) In these cases one is 
often interested in the particular integral part of the solution only; the com- 
plementary function then corresponds to a transient motion which eventually 
becomes negligible. 

We can deduce the two following results from (29) and (39), or directly 
from (28) and (38). If (33) and (34) are satisfied, then 


Y= O(Liw —A)*R’ oe”, (40) 
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(for iw =:/,). This result includes the special case considered by GuRTIN [4], 
in which A,, A,, and A, are symmetric and all the latent roots are complex. 

Otherwise, the above normalisation condition need not be fulfilled in order 
to use the result: 


gq = Q(Lo + 4A) 4 Ap (Lo + 6 An) Ro Gy et (41) 


The power of the results (40) and (41) lies in the fact that once Q, R, and A 
have been evaluated, q can very easily be calculated for a large number of 
values of w. I.e. it is only necessary to invert diagonal matrices for each new 


w to be considered (having once inverted Ay for (41)). 


6. Significance of latent vectors 


It is worthwhile to demonstrate that motion can occur in any one mode 
of free vibration (& (t) =0) without any tendency to disturb the system in the 
other modes. This is a well known property of conservative systems but does 
not appear to have been clearly stated for non-conservative systems. Such a 
mode of free vibration is described mathematically by a latent vector q,, of (2) 
with J=2. Physically, of course, we should only use the real part of q,, e’™’ to 
obtain the displacement of the system due to this mode at time f. 

In order to demonstrate this we suppose that the system is first brought 
to rest and held in the position given by the real part of q,,. It is then set in 
motion by a set of impulses which give rise to a velocity distribution given by 
the real part of A,,q,,. These one-point boundary conditions are to define the 
(2m x1) column vector e of equation (27). On putting ‘=0 in equation (31) 
with y=1 and r=0O, we obtain 

eel 


Q. Yn 
which has the obvious solution in which € has a one in the m™ position and 
zeros elsewhere. On substituting this vector € in equation (30), we see that 


at any subsequent time the displacement of the system is given by the real 
part of 


Y=nem’. 


Acknowledgement. It isa pleasure to acknowledge my debt to Prof. A.M. OSTROWSKI 
for his valuable criticisms of the first draft of this paper. In particular, for his neat 
proof of Theorem 3, which I have presented above. 


References 


[1] TurnBu.ti, H. W., & A.C. ArrKEeNn: Canonical Matrices, p. 160. Blackie 1932. 


[2] Frazer, R. A., W. J. Duncan & A. R. Coxrar: Elementary Matrices, jos slie/i5. 
Cambridge 1955. 


[3] Foss, K. A.: J. Appl. Mech. 25 E, 364 (1958). 
[4] GurtTIN, M. E.: J. Appl. Mech. 26E, 461 (1959). 


Department of Mathematics 
University of Malaya 
Singapore 


(Received August 3, 1960) 


Stability in the Large of Autonomous Systems 
of Two Differential Equations 


I. H. Murti 


Communicated by E. LEIMANIS 


1. Introduction 
In this paper’ we discuss the stability in the large of some special systems 
of two differential equations with the help of qualitative methods (of course 
in combination with the construction of the Lyapunov function). 


In Section 2, we consider the system of equations 


ae =f(o) where o=cx—dy 
and obtain certain sufficient conditions which guarantee the stability in the 
large of the trivial solution of the above system. 

In Section 3, we solve the problem of AIZERMAN for the following systems of 
two equations: 


& =ax+i(y) 
(1) 
dy _ 
aR 
and 
@ =} (x) ay 
(2) 
ied 
fg OH OY. 


For system (1) MALKIN [JI] showed that the trivial solution is asymptotically 
stable in the large under the conditions a+c<0, (acy—bf(y)) y>0 for y0 


and 


a bianco as |y| >+ ©. (3) 


1 This paper forms a part of my doctoral thesis submitted in the Department of 
Mathematics of the University of British Columbia, Vancouver, B.C., Canada. 
I express my thanks to Professor E. Leimanis for his help. 
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We prove a similar theorem without the requirement of (3). In the case of system 
(2) we give a new proof of a theorem which asserts that if c2+ab=-0, then the 
trivial solution is asymptotically stable in the large under the generalized Hurwitz 
conditions. This theorem was first proved by ERUGIN. 

In Section 4, we discuss the stability in the large of the systems 


Gp ey) = ie 2(%) +ey 
d 
Sah) +ho), Gpabetey 


2. Stability in the Large of the System da/dt = F(x,y), dy/dt =f (¢) 
In this section we shall consider the system of equations 


Oy s. ny ; 
SF eee (o), T=C4— ay 


where c and d are constants, c=-0; the functions F(x, y), f(a) are continuous 
and F(0, 0)=0, /(0)=0. Besides, the fulfilment of conditions of uniqueness of 
the solution x=0=y is assumed. 

The above system was considered in the works of Ersnov ([3], [4]) and 
Krasovskil [10]. Following KRASovSsKI, we transform the system (2.1) to the 
following form (2.2) by the change of dependent variables expressed by the 
relations 

o=cx—dy 


as Vx 


The above transformation is non-singular because c is assumed not equal 
to zero. We then have 


do 
cyt p(s, y) 

(2.2) 
= 10) 


where 


o(0,9) =eF (422, ») — aft. 


KRasovskii constructed the following Lyapunov function for the system (2.2): 


Vie.y) = I Ho) do — f p(0, ») ay, 
and using Theorem 4 (BARBASHIN [2]) proved the following theorem: 
Theorem 2.1. J} the conditions 
of(c)>0 for 3G14- 0 (2.3 
y p(0, vy) <0 for “ys=0 3 
olp(s,¥)— 90, ¥)]<0 for «+0 (2.5) 
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and 


+00 
J f(a) do 


are satisfied, then the trivial solution x=O0=y of system (2.1) 1s asymptotically 
stable in the large. - 


=e | F 910.9) dy| = (2.6) 


It may be remarked that conditions (2.3) and (2.4) can be replaced by the 
following conditions: 
of(c) <0 for o +0 (2.3’) 


yo(0,V) >O0 tor yy =450. (2.4’) 


Ersuov [4] claimed that Theorem 2.1 holds without the requirement of 
conditions (2.6). In fact he stated the following theorem: 

Theorem 2.2. I? conditions (2.3), (2.4) and (2.5) are satisfied for the system 
(2.1), then the trivial solution of system (2.1) ts asymptotically stable in the large. 

The following example shows that conditions (2.6) cannot be removed in 
general. 

Example. Consider the system of equations 


Gay Hx) 
Ne (2.7) 
ye (x) 
where /() is defined as below 
e—2% 
I( ees re | oral 
es 


Obviously, xf(x)>0 for x=+-0 and f(0)=0 


y 7(0,¥)=y(—y—#(0)) =—y2<0 for y +0 
o[y(c,¥) — 9(0, y)]=x(—y—f(™) +y)=— xf (x) <0 for x0. 
Moreover it is not difficult to show that f(x) is continuous and satisfies the 
Lipschitz condition. Thus all the conditions of Theorem 2.2 are satisfied. We 


show that the trivial solution of this system is asymptotically stable in the 
sense of LyAPUNOV but not in the large. The stability in the small follows from 


the following Lyapunov function, 


(2,9) = Jie) de + by 


We now show that there exist trajectories going to infinity for 7>-+ oo. 

It is easy to verify that y= —e~“* is a particular integral of the system on 
the interval 1<x<oo passing through the point (1, —¢7) at t=0. We show 
that along the curve y=—e~”, a >0O, 7.e., x increases with the increase of 


time. 
dx - adi ¢ af 3 
dt =—y—](x)= + é@ 4 "eran se age for xz 1. 
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5 aa e-* : = eng ee 
We integrate cere along the trajectory y=—e~* and have 


x 


t x t 
[Eeas= fat or e* + %| =4| or x+e*—e Vth 
1 0 
0 


e7* 

1 

From the last equation it follows that as >-+- oo, x—>-+ oo, t.e., the positive 

half trajectory y= —e~“* of the system (2.7) tends to infinity as >-+ co. Hence 
it follows that the trivial solution is not asymptotically stable in the large. 


2.2. Let us consider the system 


ax 
Se) 
ot (2.2.1) 
%Y = fo), =ex—dy, e+0 
under the conditions: F(0,0) =0, #(0) =0 (2.2.2) 
of(o)>0 for o+0 (2.255) 
o[g(o,y) —v(0,¥y)]}<0O for o==0 (2.2.4) 
yo(0,y)<0 for y+0 (2.2.5) 
where Par 
(oc, ¥) =cF( : ,y) —df(o). (2.2.6) 


Besides, the fulfilment of conditions of uniqueness of the solution x=0=y 
is assumed. As before we reduce the system (2.2.1) to the following system: 


ites (o, y) 
iy (2.2.7) 
age or AC) 


We consider the positions of the curves represented by the right-hand sides 
of (2.2.7) on the (a, y) plane, 1.e., of p(o, y)=0 and f(a)=0. Since of(o)>0 
for o+0 and f(0)=0, f(¢)=0 only when o=0, i.e., {(o)=0 represents the 
y-axis. We now turn to the curve represented by y(a, y)=0. We observe that 
y(0, 0)=0 and on the o-axis y(o, y)>0 for o<0, p(o, y)<0 for o>0; on the 
y-axis y(o, y)>0 for y<0, p(c, y)<0 for y>0. From these facts it follows 
that p(o, y)<0, o>0, y>0 and y(o, y)>0 for «<0, y<0. In deriving these 
conclusions we have made use of the conditions (2.2.4) and (2.2.5). Thus it 
follows that y(o, y) changes sign in the second and fourth quadrants, and hence 


the curve p(o, y)=0 lies in the second and fourth quadrants. It is easy to see 
that 


ape p(o,¥) >0 for the points lying to the left of the curve y(o, y)=0 
an v(o,¥) <0 for the points lying to the right of the curve y(o, vy) =0 
dy ; : : 

Fr =l(o)>0 for the points lying to the right of y-axis 

dy 


Fis f(o) <0 for the points lying to the left of y-axis. 
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The curve p(o, y)=0 and the co-ordinate axes decompose the plane (o, ) 
into six regions. y(¢) is maximum on the y-axis for y>0O and minimum for 
y <0; o(t) is maximum on the curve (o, y) =0 for ¢>0 and minimum for ¢<0. 
The direction of motion is indicated in Fig. 1. We introduce polar co-ordinates 


o=rcos?, y=r7sin?, 
then Seas a 
7=ocos?+ysin# 
and ° or ‘ 

73 =—csin¥+ycos@. 

The signs of 7 and § in different regions are given as below: 
(1,4) #% may be 0, d>0 Y 
(2,5) # may be 20, #>0 
(3,6) *#<0, may be 20. 


The Lyapunov function for the 
system (2.2.7) is 


V(o,9) =J f(@) do — J 90,9) dy. 4 


Its total time derivative in view of 
equations (2.2.7) is 
V =f(e) [y(o, y) — ¢(0, 9)] 
=O 10is=o-e 0 y' 


=O ior C=O- Fig. 1 


oe. 


Obviously, V(o, y) is positive definite and V is of negative sign. Now ape (0, y) 


for o=0 and it is different from zero unless y=0. This means that o=0 does 
not contain any other positive half trajectory of the system (2.2.7) than the 
origin. Thence follows the asymptotic stability in the sense of LyAPUNOV of 
the trivial solution of system (2.2.7). Thus we have shown 


(i) the point (0, 0) is the only point of equilibrium, 
(ii) the unperturbed motion is asymptotically stable according to LYAPUNOV, 


(iii) there are no periodic motions, since for (2.2.7) a Lyapunov function has 
been constructed. 


Since the o positive half axis is intersected by the motions in one direction 
only, it can be taken for the straight line L(0, co) appearing in Theorem 1.3 
(ErucIN [5]). Thus all the conditions of Theorem 1.3 (ERUGIN [5]) are satisfied 
except the fourth, 7.e., the motions having bounded polar angles are bounded. 
We now indicate what additional conditions are to be imposed in order to realize 
the fourth condition of Theorem 1.3 (ERUGIN [5]). 

Let a motion M(t) start in a region, say (6). Any motion M(?) started in the 
region (6) or entering this region either tends to the origin or goes out of this 
region and enters the region (1). This follows from the fact that 7<0 in this 
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region. After entering the region (1) the motion M(t) either crosses the y-axis 
or tends to infinity along the y-axis, but it cannot go to the origin since o>0: 
We are thus led to impose the following 


Condition A. We assume that (a, y) satisfies such conditions called A that 
the motion entering or starting in the regions (1) and (4) leaves these regions 
with the increase of time. 

This condition guarantees that there are no motions with bounded polar 
angles in the regions (1) and (4). We impose another condition called B. 


Condition B. We assume that g(a, c’)=0 has a solution for all c’. 

If Condition A is satisfied, then the motion M(¢) enters the region (2). We 
now show that the motion M(t) leaves the region (2) with the increase of time 
if Condition B is satisfied. In fact if g(o,c’)=0 has a solution for all c’, 2.e., 
if yc’ intersects the curve g(a, y)=0, the motion after entering the region 


(2) can neither cross the line y=c’, since a< 0 in this region go to the 


origin, because 8>0 in region (2). Therefore the motion must intersect the 
curve y(o, y)=0. Similar reasoning can be carried out in the regions (3), (4) 
and (5). The above analysis shows that there are no motions with bounded 
polar angles in the regions (1), (2), (4) and (5) and that the motions with bounded 
polar angles can occur only in the regions (3) and (6) and as proved above they 
are bounded. Thus any motion with bounded polar angle is bounded. 
Condition B was imposed with a view to ensure that in the regions (2) 
and (5) there are no motions with bounded polar angles. This can also be achieved 


if we assume that in these regions 0 > >0. This condition we call Condition C. 
We now collect these results in the following 
Theorem 2.2.1. If Conditions (2.2.2)—(2.2.5) are satisfied and if either Con- 


ditions A and B or Conditions A and C are satisfied, then the trivial solution of 
(2.2.1) 1s asymptotically stable in the large. 


It may be remarked that Conditions B and C are independent, 7.e., if Con- 
dition B holds then Condition C may or may not hold or vice versa. In the linear 
case, 7.e., when the right-hand sides of system (2.2.1) are linear, however, both 
hold. 


Let us consider a few examples now. 


Example 1, Consider the system 


a (2.2.8) 


Obviously, this system satisfies the conditions (2.2.2) —(2.2.5). We verify 
that in this case Condition A is satisfied. Let us for the sake of definiteness 
assume that the motion M(j) is in the region (1). We assume that the Condition A 
is not satisfied. Then the motion M(t) goes to infinity as t+-+co and during 
this y(t)>-+ 0c. From the first equation of (2.2.8) follows that dx/dt becomes 
infinitely large but negative, which means that MW (¢) cannot remain in the first 
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quadrant, and hence our assumption is not true. It is not difficult to verify 
that the straight line yc does not intersect the curve 
— Ve —=1—(0 for all’e*(see Fig. 2). 


Hence, the Condition B is not satisfied. However, the Condition C holds. In 
fact rO=— xsin 0+ ycosd=—sin#(—y+e-*—1)+xcos® and 


= ee y+e-*—1)+ cos? O>e>0 
in the regions (2) and (5). This es- J 


tablishes the asymptotic stability in 
the large of the trivial solution of 
(222N8). 


Example 2. Consider the equation 


fa +9(S)e0) +70) =0. 229 


This equation can be thrown into 
the system 


dx ‘ v. 
So 4) BAY) TY) 
dy 
px =*% 
dt y' 
by writing x for dy/dt. Fig. 2 


We assume that the following conditions are fulfilled: 


yi(y)>0, 7f(0)=0 


and either 

h(x) >0O for x+0 and g(y)>0 forally 
or 

Wiajy—20" for x =-0" and ¢(y) <0" for ‘all y 
where 


hy (x) P(*) for x=+0. 


The Lyapunov function for this system is 


V(x, y) =[xdx+Ji0)dy. 
Then 
V =x (— y(x) g(y) —f(y)) +4(y) * 
= — x 9(x) g(y). 


It is not difficult to see that in this case Condition C of Theorem 2.2.1 is 
satisfied. Thus if Condition A is also satisfied then the trivial solution of (2.2.9) 
is asymptotically stable in the large. 
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Example 3. Consider the system 


© = (x) + fal) 
; (2.2.10) 
d 
LEeP UAE) 
For the system (2.2.10) the V function is 
=Sh( )dx — Stel 
Th 
* V(x, 9) = A(x) pla). 
We subject (2.2.10) to the following conditions: 
g(0)=0, hg(x)<0O for x=+0, where hg(x) = (*) for x = 0) 
- 2.2.11) 


f,(0) =O and xf,(*) >Oforx +0, /f,(0) =Oand yf,(y) <0 for y +0. 
The trivial solution of (2.2.10) is thus asymptotically stable in the large if in 


addition to (2.2.11) the Conditions A and B or A and C are satisfied. 
Example 4. Consider the equation of the second order 


dy 
dp 7 Aer "2 ar VE (2.2.12) 


We write it in the form of a system of equations 
Sait = [ry) a= ey) 


We assume that 
g(0)=0, yeg(y)>0 for y=+=0 


f(a O1niors x == OF 


and 


It is easy to verify that asymptotic stability in the large holds if Condition A 
is satisfied. 
3. The Problem of Aizerman 
In 1949 AizERMAN [1] proposed the following problem. 
Let there be given a system of linear differential equations 


dx “ d fe : 
Fr = DL BiH OM, = av F (¢ SSeDN3 Sexgiteyis (3.1) 
j=1 
Suppose that for the given constants a;; (t=4,...,”, 7=1,...,) and for an 


arbitrary value of a from the interval oe B all a roots of the characteristic 


equation of (3.1) have negative real parts. Let ax, be replaced by /(%,) in (3.1). 
We then have 


dt = Di GAT), SE = Say 3 (1=2,...,). (3.2) 
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It is required to find out whether the trivial solution x,=*,=---=x,=0 
of system (3.2) is asymptotically stable in the large or not, for an arbitrary 
choice of the continuous function /(x,), which reduces to zero for x,=0 and 
which satisfies the inequality 


OX < Xf (%) <B x for x, 0. (3.3) 


The answer to the above problem is in the negative (PLiss [12]). The interest 
in the problem is revived if we slightly change the above problem and ask our- 
selves the following questions. For what values of a,; is the answer to the above 
problem in the affirmative and for what values in the negative? If the trivial 
solution is not asymptotically stable in the large under the generalized Hurwitz 
conditions (i.e., the condition (3.3) and /(0)=0), what additional assumptions 
should be made on /(%,) so that the trivial solution becomes asymptotically 
stable in the large? We shall discuss here the case when n=2. The case n=3, 
k=2 has been completely solved recently by Piiss [13]. First of all we take 
n=2=h, 1.e., we consider the system 


sapite OM lV) ry 
ay ; 
rapt aor be 
under the assumptions that 
HOV Onna ret, 0, 20 -On(yi> 0, Vy = 0 (3.5) 
where f(y) 
h(y) = 12 for V0: 


Besides, the uniqueness of the trivial solution is assumed. System (3.4) was 
considered by N. P. Erucin [6] and I. G. MarkIn [1/1]. 

MALKIN [11] showed that the trivial solution of system (3.4) is asymptotically 
stable in the large if for sufficiently large values of |y| the inequality ac — bh(y) >e 
holds. This condition can be relaxed to the condition that 


i ides eo ine Rees as |y| +. (3.6) 


0 


We show that asymptotic stability in the large of the trivial solution of 
(3.4) holds without the requirement of condition (3.6). For b=0, from (3.5) 
follows that a<0, c<0, and an immediate integration of the system (3.4) shows 
that the equilibrium is asymptotically stable for arbitrary initial disturbances 
and for an arbitrary choice of the function /(¥). 


Let b-+0. We introduce new variables defined by 


x’ =bx—ay 


yy; 
then 
Ce dx ay _ sh b = — y’(ac— dh(y’ 
# = bg DY —b(ax+fly)) albu toy)=— 9 (¥’)) 
YY ba teyau tay toy =x t(at+o)y’. 


9* 
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Thus the system (3.4) reduces to the following: 


4% = —y(ac—bh(y)) (3.7) 
arg 2). 


We represent the curves —y(ac—bh(y))=0 and x+(a+c)y=0 obtained 
by putting the right-hand sides of (3.7) equal to zero on the (x, y) plane. The 
first of these represents the straight line yO, and the second the straight line 


y= — au It is easy to see that 
ee =— y(ac—bh(y))>0_ below the x-axis 
oe =—y(ac—bh(y))<0 above the x-axis 
DY oat (ate) y>0 below the straight line %*+(a+c)y=0 
a =x+(a+c)y<0 above the straight line x+(a+c)y=0. 


The direction of motion is indicated in Fig. 3. The straight line x+(a+c)y=0 
and the co-ordinate axes divide the plane (x, y) into six regions. We introduce 
polar co-ordinates x*=rcosy, y= 
ysing. Then r7=x cos g+ysin 
and rg=—xsing+ycosg. It is 
not difficult to verify that in differ- 
ent regions the signs of 7 and @ are 
given as below: (1,4) 7 may be =0, 
p >0; (2,5) 7<0, @ may be = 0; 
(3,6) * may be £0, p>0. 
The motion started in the region 
(6) must cross the axis of x. This 
follows from the inequality 
dy 
ae 
in the region (6). After entering the 
region (1), it cannot remain there 
Fig. 3 and must cross the straight line 
; x-+ (a+ c)y=O with the increase of 
time, since x=/ intersects the straight line x-+(a+c)y=0 for all 1. The Lya- 
punoy function for the system (3.7) is easily seen to be 


=%*+(atc)y2d>0 


ts 
2 Vise wits 2) (a oteNe ENO Nianey 


Repeating the same argument as in Theorem 2.2.1, we arrive at the following 


Theorem 3.1. The trivial solution of system (3.4) is asymptotically stable in 
the large under Conditions (3.5). 
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It may be noted that we may or may not have ¢>e>0 in the regions (1) 
and (4). Before studying the case n=2, k=1, we discuss the system 


raetiene soit 
(3.8) 
G=— 82). 


This system is a particular case of the system (2.2.1). Since it is quite an 
important system, we discuss it independently. The conditions to which (3.8) 
is subjected are 


y 8) 
« F(x) >0, y 
wen) OL 10r 4 0 (3.9) 
F(O}=0 = 2 (0), 


The V-function in this case is 


2V = 97 +.2.f g(x) d%003.10),, % 
and ° 


a = 9(—8(x)) +8(x) (y Fl») 


=p 9 (x) bs )i< 0. fOr aX =-0 


—\() OLN Oe 


The straight line «0 does not 
contain a positive half trajectory of Fig. 4 
the system (3.8) except the origin. 
In fact, dx/dt=y for x=0. Hence it follows that the trivial solution is asymptoti- 
cally stable in the small. Since there exists a Lyapunov function V for the 
system, there cannot be any limit cycle. Also we note that the origin is the 
only point of equilibrium. It is not difficult to see in virtue of the Conditions 


(3.9) that 
— =y—F(x)>0 to the left of the curve y = F(x) 


Tae ee F(x) <0. to the right of the curve y = F(x) 


ee =— g(x) >0 to the left of the y-axis 
= — g(x) <0 to the right of the y-axis. 


The curve y=F(x) and the co-ordinate axes divide the plane (x, y) into six 
regions. The direction of motion is represented in Fig. 4. We introduce polar 
co-ordinates x=r7 cos y, y=rsin y. Then 


7=%cosp+ysingp 
rp=—xsing+ycosg. 


The signs of 7 and @ in the regions are given as below: (1,4) 7 may be = 0, 
@ <0; (2,5) *<0, p may be £0; (3,6) * may be £0, P<0. 
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Let us follow the motion M(t) after it intersects the negative half y-axis. 
Let us assume that y= A intersects the curve y=F(x) for all A, 7.¢., F(x) =A 
has a solution for all A; then, since in the region (4) ¢<O and y is increasing, 
the motion M(é) cannot remain in the bounded region OL M. It definitely can- 
not cross LM and OL. Therefore it must cross the curve y=F (x) and enter 
the region (5). In the region (5) it either enters the origin or the region (6). 
This follows from the fact that 7 <0in this region. Since dx/dt= — F(x) + y2d>0 
in the region (6), the motion must leave this region and enter the region (1). 
Similar argument holds for the regions (4), (2) and (3). This shows that there 
are no motions with bounded polar angles in the regions (1), (3), (4) and (6). 
The motions with bounded polar angles can only occur in the regions (2) and (5), 
and as proved above they are bounded. Thus it follows that all motions with 
bounded polar angles are bounded. For the straight line L(0, oo) can be taken 
the positive half y-axis. Thus all the conditions of Theorem 1.3 (ERuGIN [5]) 
are satisfied, and we have the following theorem: 

Theorem 3.2. If the Conditions (3.9) are satisfied and F(x)=A has a solution 
for all A, then the trivial solution of (3.8) 1s asymptotically stable in the large. 

We remark that by using Theorem 4 (BARBASHIN [2]) the stability in the 
large holds if 


Sg(x)dx—>oco as |x| >. 

0 
This is required in order to make the V function in (3.10) infinitely large. We 
further remark that the boundedness of the solution in Theorems 3.1 and 3.2 


could be shown by making use of Lyapunov functions V, as has been suggested 
by Professor J. P. LA SALLE. 


The widely discussed equation 
d* x dx 
Fg ted eg doe) (3.11) 


can be dealt in the same way as the system (3.8), since (3.114) can be transformed to 


OF aay — Fie) 
dy - 
ap = 8(*) = where F(x) = f f(x) dx. 


0 


We shall now discuss the Aizerman problem for the case n=2, k=1, i.e., 
the system 
d 
4 (3.12) 
ad 
under the conditions 


f(0)=0; c+h Osa nee Ris 
ee e+h(x)< ch(x)—ab>0 for x+0 (3.13) 
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This system is discussed in the works of Erucin ([4], [7]) and Markin [17]. 
It was proved by ERUuGIN that if c?-+a@b-+0, then the trivial solution of (3.12) 
is asymptotically stable in the large under conditions (3.13). The proof by 
EruGIN of this theorem is quite lengthy, so we give here a short proof of the 
theorem. It will become clear from our proof why the stability in the large does 
not hold under conditions (3.13) in the case c?+ab=0. ErucIn also showed 
that for the stability in the large in the case c?+ab=0, it is sufficient that 


Tim | f (eft) —abx)dx+c f(x) —abs| == sLico 

3 (3.14) 
lim LS (cHx) —abx)dx—cf(x) +abr| = od, 
>= OO 0 

KRASOVSKII [8] showed that con- 
ditions (3.14) are necessary as well 
as sufficient. The question of the 
region of stability in the case where 
the stability in the large does not 
hold is discussed by Puiss [12]. 

We assume that c?+ab=+0. Let 
a=0; then immediate integration of 
the system yields the stability in the 
large of the trivial solution x=0=y 
of system (3.12). Now let a=+-0. We 
introduce new dependent variables 


age 
y =ay—cx. 
Then 
dx’ ax h ' , 
fae te wea 
CS OESy hee = 
ee, = 4 (D K+ CX) c(f(x) +ay) 
=abx' —cf(x’). 


The system (3.12) is reduced to the form 


Fm A ate ed 
= — x(ch(x) —ab). 


Comparing it with system (3.8), we have 
F(x)=—x(c+h(x)) and g(x) =x(ch(x) — ad). 


The curve y+-cx-+/(x)=0 and the co-ordinate axes divide the (x, y) plane 
into six regions. The direction of motion is represented in Fig. 5. We consider 
the following cases. 
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Case i. c<0. Consider the straight line 


Pbatinida (3.15) 
C 


If c2-+-ab>0, then the straight line (3.15) lies in the first and third quadrants. 
Let us see how it is situated on the (x, y) plane with respect to the curve 


ytex+ f(x) =0. (3.16) 
Let y® and y®) denote the ordinates of (3.16) and (3.15) respectively; then 
cttab 
Mtb a tm tag CA sey 
—c?x—cf(x)+(c?+ab)*  —x(ch(x)—ab) >0 
a C C 


according as x 20 


i.e., the curve (3.16) lies above the straight line (3.15) in the first quadrant and 
2 
below in the third quadrant. Obviously, the straight lines y= A and y= — Bae % 
intersect, whence follows that y=A and the curve y=—cx—f(x) intersect 
for all A. Hence the trivial solution of (3.12) is asymptotically stable in the 
large in this case according to Theorem 3.2. 
eee 


If c?+ab<0, then the straight line y= — 7 TOF lies in the second and 


fourth quadrants, and no such conclusion as above can be drawn. However 
we can prove that in this case ¢ < — e<( in the regions (1) and (4) which ensures 
the stability in the large. 


According to Conditions (3.13) 
e+h(x)<0, ch(x)—ab>0 for *+0 
1.€. 


h(n =o. Te) Bus 


but since c?-+-ab<0, the above inequalities are satisfied if we take h(x)<—ce. 
/(x) can then be written as 
f(x) =—¢ x —a(x) where «(x) 20 according as 
%~Z0. 
Let us calculate @: 
ry =—sin 9 (f(x) +ex+ y) — cos p(ch(x) —ab)x 
=— sin p(y — «(x)) — cos p[e(—¢ x — a(x)) — ab x| 
= —sin p(y —a(x)) + cos p[(c2?+ ab) «+ cw (x) | 
sin yp 5 
aeons (v — a(x) +e [(c2-+ ab) x + ca(x)|<—e<0 


in the regions (1) and (4). 
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Case it. c>0. In this case c?-+ab<0, because from Conditions (3.13) 
a h(x)<—c. 


Let h(x) be taken as h(x) = —c—a(x). Then 


or 


O<a(x)<— 
Let us consider the region (1). If g¢ « —e<0O, then 


lim —* [(2 + ab) x texa(s)|=0, te, 


X%—>0O 


lim x«(x)=co, whence it follows that the straight line y=A intersects the 


curve 
y=—cx—f(x) forall A. 
In fact, 
—A=cxux+f(x)=cx+(—cx —xa(x)) =— xa(x); 


therefore xa(x)=A. If xa(x)=A does not have a solution for all A, then 


lim x a(x) = finite = D 


x0 


lim — [(c? + ab) x+c¢xa(x)| 2 lim [— (c+ a0) x] + lim (— c x«(x)) 


X%—> 0O %—> 00 


= lim — x(c? + ab) — lim (eazon(a) ys 


X—> Co X—>CO 


whence follows that p< —e<0. 
Case wit. c=0. The System (3.12) reduces to 


oy wi a 
dy _ 
ee 


and the Conditions (3.13) reduce to /(x) <0, ab<0. 
For this system we take the following Lyapunov function 


2V=—abx+ y?, 
the derivative of which is 


This V function satisfies all the conditions of Theorem 4 (BARBASHIN [2]). Com- 
bining all these results, we have the following theorem: 
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Theorem 3.3. I} for the System (3.12) c?+ab=50, then the trivial solution 1s 
asymptotically stable in the large under C onditions (3.13). 
In the case c2-+ab=0, we remark that c<0 and from the proof of stability 


; : b 
in the Case 7, we find that the straight line y= — ss ae x coincides with the 


x-axis, and hence we can neither say that the straight ie y=Ai intersects the 
curve y+cx+}(x)=0 for all A nor ¢<—e<0O in the regions (1) and (4). 


4. A Generalization of the Problem of Aizerman 
In this section we discuss the system of differential equations given by 


9% ant f(y) 
* (4.1) 
Y= f(z) + ey 


under the conditions 
abheox0, ac—hly) dls) S00 flor ey 0) VEO, Oy — 6 
= f,(0) =0 


where h,(y) and h,(x) are defined by f,(v)=yvAy,(y) and fy (x) =x hy (x). 

The above system was first discussed by KRAsovsk1 [9] who obtained certain 
theorems regarding the stability in the large of the trivial solution x =0=y 
of (4.1). We prove here the following theorem: 


Theorem 4.1. If either h,(y)>0 for y==0 and h,(x)>0 for x0 or h4(y) <0 
for y=+-0, hg(x)<0 for x=0, then the trivial solution of (4.1) is asymptotically 
stable in the large under Conditions (4.2). 

Proof. We assume /,(y)>0, y=+0, and h,(x)>0, x=+0. From Conditions 
(4.2) follows that a<0, c<0. We consider the curves 


ax+f(y)=0 and f,(x)+cy=0 (4.3) 


Since /,(y)>0, y=-0 and a<0, the curve ax+f,(y)=0 lies in the first and 
third quadrants. The same is true for the curve f,(x)+cy=0. We further 
note that the curve ax-+/,(y)=0 lies above the curve /,(x)-+-cy=0 in the first 
quadrant and below in the third quadrant. In fact, if y, and y, denote the 
ordinates of the curves in (4.3) respectively, then 


(4.2) 


ae Ae ue ED ee as a(*). 20 accordingas x20. 
It is easy to see that 
a =ax+f(y)>0 tothe left of thecurve ax +fi(y) =0 
oy =a%+h(y)<0 to the right of the curve ax+f,(y)=0 
co =f(x) +ey<0 tothe left of the curve cy+f,(x) =0 
dy 


as =f,(*x) +cy>0 to the right of the curve CY +fo(x) = 0. 
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The function x(/) is maximized on the curve ax+f,(y)=0 for y>O and 
minimized for y<0; and the function y (¢) is maximized on the curve f.(x) +cev=0 
for y>0 and minimized for y<0. The curves and the direction of motion are 
shown in Fig. 6. As before we introduce polar co-ordinates x=r cos Y, V=rsin @. 
Then 7=%cos y+y sin 9, 


71p=—xXsing+ycos@. 


The signs of 7 and @ in different regions are (1,5) * may be =0, ¢ <0 (2,6) *<0, 
@ may be £0 (3,7) * may be £0, p>0 (4,8) *<0, @ may be Zs 

From Theorem 2.1 of ERUGIN’s work [5] it follows that there is at least one 
integral curve going to the origin, the only point of equilibrium, in each of the 
regions (4) and (8). Other motions 
started in the regions (4) or (8) either 
go to the origin or enter the regions 
(3), (5) or (4), (8) (since *<0, p may 
be = 0 in (4) and (8)). Let us suppose 
that the motion enters the region (1). 
We show that the motion leaves this 
region with the increase of time. Let 
Cy. be the least upper bound of h,(x). 
Cy is finite, otherwise ac — hy, (y) ho (x) 
>0 will not hold. We consider the 
straight line c.x+cy=0. It is easy 
to verify that this straight line lies 
above the curve f,(x)+cy=0 and 
below the curve ax+f,(y)=0 in the 
first quadrant. In the third quad- 
rant the positions are reversed. We 
consider the region bounded by y=, the curve ax+/,(y)=0 and the positive 


Fig. 6 


half y-axis. The motion cannot intersect the straight line y=/ since <0. 


It cannot go to the origin since g¢<0. Therefore it must intersect the curve 
ax+f,(y)=0 and enter the region (2) where it goes to the origin with the 
increase of time. For region (3) we take the straight line x=/’ and similarly 
show that the motion crosses the curve /,(%)-+-cy=0 and enters the region (2) 
where it goes to the origin as ¢+-++ co. Similar reasoning holds for the regions 
(5) and (7), which completes the proof of Theorem 4.1. 


In the same way we can prove the following theorem for the system 


(4.4) 


under the conditions 


h(x) +ce<0, *=+0, chy(x)—ha(y)b>0, x40, y+0, 
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where hy (x) = A), ply clones bi fale y +0. (4.5) 


Theorem 4.2. If bh.(y)>0 for y+-0 and the Conditions (4.5) are satisfied, 
then the trivial solution of system (4.4) is asymptotically stable in the large. 
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I. Introduction 
Much is known about the behavior of the solutions of the differential system 


An +H w= (0%), 160) =0, 

in a neighborhood of « =0 provided that the eigenvalues of the constant matrix A 
have nonzero real parts and f(t, x) 1s small compared with x when x is small 
(see, e.g. [5, Ch. 13]). However, if some of the eigenvalues of A have zero real 
parts, then the theory is not so well developed, although some great contributions 
have been made by A. Liapunov [/4] and more recently by DyHMAN [8], MAL- 
KIN [16], PERstDskiI [18] (see also [/1, Ch. 6]). The purpose of the present paper 
is to discuss a particular class of equations of the above type (and even more 
general ones) where some of the roots of A have zero real parts. The relationship 
of this system to the ones considered by the above mentioned authors is given 
in Section III. The system of equations considered arises in a natural way as 
variational equations for invariant surfaces which consist of families of periodic 
solutions. 

As a motivation for what follows, consider the system of autonomous differ- 


ential equations 
(1.1) te Xx) 

dt 
where x, X are n-vectors, X € C2 and suppose that * = °(¢) is a periodic solution 
of this system. It is known [5, p. 323] that this periodic solution is asymptotically 
orbitally stable (7.e., the closed curve Cy in x-space defined by x = x°(t) is asymp- 
totically stable) if 7 —1 of the characteristic exponents of the linear variational 
equation 


dy Ox (*°) 
(1.2) Geter 


have negative real parts (necessarily one characteristic exponent is zero). If 
x(t) is a solution of (1.1), then x°(¢-+ gy), where y is an arbitrary constant, is 
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also a solution and, therefore, the function x°(¢-++q) defines a cylinder, C, in the 
(x, t)-space. If the closed curve Cy in the x-space is stable, then clearly the cylinder 
C in the (x, é)-space is stable. But, by discussing the solutions in the (x, t)-space 
more information can be obtained. In fact, under the same hypotheses on the 
linear variational equation (1.2), not only is the cylinder asymptotically stable, 
but for any solution x(¢) of (1.1) with x (0) sufficiently close to Co, there exists 
a MY such that x(t) — x°(¢ +) +0 as > 0 (see, e.g., [5, p. 323]). This is called 
asymptotic orbital stability with asymptotic phase. 

Now, suppose that (1.1) has a (k +1)-parameter family of periodic solutions 
x = x°[w(b)t-+,b], where b =(b,,..., ,) and @ are real constants. In the (x, ¢)- 
space, such a family of solutions defines a k-parameter family of cylinders C, 
or a (k-+2)-dimensional manifold Nt of solutions. For a fixed value of b, the 
variational equation (1.2) can be defined and, in general, has k +-1 characteristic 
exponents equal to zero. Is it possible to obtain some type of stability for this 
manifold IN if m—(k-+4) characteristic exponents of (1.2) have negative real 
parts for every value of the vector d? Under very weak conditions on Jt this 
question is answered in the affirmative in Theorem 5.1. In fact, it 1s shown 
that for any solution x(t) of (4.1) with (0) sufficiently close to WW, there exist 
Qo, 09 such that x(t) — x°[w (bo) £+ Mo, 29] > 0 as t> ov, 1.€., x(t) approaches the 
cylinder C, with asymptotic phase. 

Another question that naturally arises is the following: suppose that system 
(1.1) has an asymptotically orbitally stable periodic solution x°(f) and the right- 
hand side of (1.1) is perturbed by a function X*(t, x) which approaches zero 
as t-> oo uniformly for all¢=2T and x€ U, a bounded set in x-space. Is it true 
that the solutions of the perturbed system approach the cylinder defined by 
x°(t-+q) even though x(t) will not, in general, satisfy the perturbed equation ? 
Asymptotically autonomous systems have been previously discussed by L. Mar- 
kus [17, p. 27], and he has answered the particular question posed here in the 
affirmative provided the vector x has dimension two. The technique used by 
MARKUS seems to depend very strongly on the fact that the phase space is two- 
dimensional. This question is further investigated in this paper for vectors of 
arbitrary dimension and (k +1)-parameter families of periodic solutions of the 
unperturbed equation (see Theorems V.6, V.7 and V.15). As can be seen from 
these theorems, the rapidity with which X*(t, x)>0 determines the behavior 
of the solutions in the neighborhood of the manifold. 


All of the results mentioned above are obtained by first introducing local 
coordinates (s, a, h), where s is a scalar, ais a k-vector, h is an n — (k +1) vector, 
in a neighborhood of the manifold of solutions of (1.1) in such a way that the 
manifold is given by h =0, a=constant, s =w (a)t-+q. The differential equations 
in the new variables s, a, h have the general form given by 


d 
> — @ (a) ++ Sy (S a, h) sac S2(s, a, h, t), 
da 
(4.3) Fr = Axls, a, h) + Ag(S, a, h, t) , 
dh 


“ap = Hs, a) h + y(s, a, h) + Hy(s, a, h,t), 
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where each of the functions S,, 4,, H, is O(||A||?) as ||2||->+0 and the functions 
Sz, 42, H, behave essentially the same as the perturbing function iX* (tt) 
mentioned above. In general, the most that one can hope to obtain is that a 
solution of (1.3) has the property that h->0,a—a constant, c,, s — [w (c,) t+] 0, 
as t-> oo, where @ is a constant. 

The discussion of the behavior of the solutions of a system like (4.3) is diffi- 
cult since each of the functions s(¢), a(t), h(t) is expected to behave in an entirely 
different manner as a function of ¢. Although it may be possible to discuss all 
of the solutions of (1.3) at the same time, we have chosen to look at the equations 
separately. For example, consider the equation 


(1.4) Be = H(s(2), a(t))h + H, (s(t), a(t), h, t) + Hy (s(t), a(d), h, 2) 


where s(t), a(¢) are given functions of ¢. For given functions H, H,, H,, how large 
can the classes of functions .Y, # be so that for every s€ 4%, aC AW, the solution 
h(t) of (4.4) approaches zero as t-> oo? A result of this type is contained in 
Lemma II.9. This seems to be a more natural way to attack the problem since 
the results so obtained may be useful in other applications than the ones consid- 
ered in this paper (see, e.g., Remark II.11). The other two equations in (4.3) 
are discussed in the same way, and all of Section II is devoted to a statement 
of these results. In the proofs of the results of Section II, the fixed point theorem 
of TycHONOvV is used systematically. 

The results of Section II are then used to obtain theorems about the behavior 
of the solutions of the complete system (1.3) in somewhat the same manner 
that the method of elimination is used to solve a system of linear algebraic 
equations. More precisely, the second equation in (1.3) is solved for a(h, s) (t) 
as a function of h,s, hC#, s© YF, where W, FY are given classes of functions. 
It is then shown that a€.V, where is a given class of functions. Then, in 
the first equation, a is replaced by a(h, s) (t), and h/ is replaced by hE # to obtain 
a solution s(h)€ Y. Finally, in the third equation, s, a are replaced by s(h) €-Y, 
a(h, s(h))€.V, respectively, to obtain a function h€H#H. The three functions 
3=s(h), a(h,3).h yield a solution of (1.3) belonging to Yx x H. These 
results are contained in Section III. One reason that the proofs of these results 
are difficult is the fact that soo as t->oo, and the asymptotic behavior of 
s depends upon the asymptotic behavior of a. If either of the variables s or a 
is absent in (1.3), then much more elementary arguments can be used (see, for 
example, J. K. Hate & A. Stokes, A Discussion of Differential Equations on 
Product Spaces, Proceedings of Symposium on Ordinary Differential Equations, 
University of Mexico, 1959). The relationship of some of the results in Section III 
to the general problem of stability in “product spaces” (see S. Lerscuetz [11, 
p. 117]) is given. 

A key requirement in this approach is that the functions a(h,s) and s(h) 
be continuous in the topologies of the various function spaces. Lemma II.14 
shows that this requirement can be fulfilled. The proof of Lemma II.14 is such 
that it is clear that this method can be applied to a system of equations, for 


any n. 
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In Section IV a transformation, valid in a neighborhood of a family of per- 
iodic solutions of (1.1), is given which transforms systems of type (4.1) or per- 
turbations of (1.1) into special cases of systems of differential equations considered 
in Sections II and III. In Section V, some theorems are given concerning the 
application of the theorems of Section III to the stability of surfaces of periodic 
solutions of (1.1) or perturbations of (1.1). Numerous examples are given to 
illustrate the theorems of Section V. 

The applications in Section V have been the motivation for writing this 
paper. The systems considered in Section V together with the transformation 
in Section IV lead directly to systems of the type considered in Section III. The 
theorems in Section III may be read independently of the rest of the paper, 
but the proofs depend ‘very strongly upon the results in Section II. 


II. Notations and basic lemmas 


Let R” be n-dimensional Euclidean space, and let £7 be the space of all 
continuous functions x(¢) from [7, co) +R”, T 20. The topology on E7 is the 
compact open topology, 7.e., a subbase of neighborhoods in £7 is the collection 
of subsets of E% of the form M(L, W) ={f€ E7|f(L) < W} where L is a compact 
subset of [T, co) and W is an open set in Rk”. Since R” is a metric space, conver- 
gence in the compact open topology is equivalent to uniform convergence on 
compact subsets of [T, oo) (KELLEY [J0, p. 230}). Also, since R” is a metric 
space and [7’, co) is the sum of a countable family of compact subsets having 
the property that any other compact subset of [7, oo) is covered by some finite 
subfamily, the space £7; is a metric space and, as R” is complete, E> is complete 
with respect to its metric (ARENs [/]). Furthermore, £7 is a linear topological 
space, and to define the topology it is sufficient to give a system of neighborhoods 
of the identity. Now, consider the family of pseudo-norms {f,(x)}, p,(”) = 
upiS UE l|~ (||, where «€ E> and ||y|| for yC R" is any vector norm. It can be 
shown that the family of sets {AV,}, A=0, » =1, 2,..., where V, = {x € E7| ,, (x) 
<1}, considered as a subbase for a system of neighborhoods of x =0, generates 
a topology on E> which is equivalent to uniform convergence on compact subsets 
of [T, oo). Finally, the sets V, are convex, and, summarizing the above, we have 

7 ts a complete, locally convex, linear, topological space (for more details, see 
Boursaki [3]). 

A set, A, in E% is bounded if the image of A under each pseudo-norm, Py (*), 
is bounded. Notice that this does not imply that any function x¢ A is bounded 
on [T, cc). An operator on E7 is a continuous function mapping E% into itself. 
An operator ¥ is compact (completely continuous) if for every bounded set 
A CE* the closure of its image under ¥%, (A), is compact. Let 


t 
I(x) (t) = f F(x(s), s) ds 
dt 
where Bf : R"x[T, co) +R" is continuous. By using Ascort’s theorem, the diago- 
nalization procedure and the fact that a sequence {x,} converges in the compact 
open topology if and only if the sequence {x,(¢)} converges uniformly on every 
compact subset of [7, oo), it is easy to show that ¥ is a compact operator. 
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In the following pages, the fixed point theorem of TycHonov [2/] will be 
used systematically. 


Theorem (Tycuonovy). Let E be a complete, locally convex linear topological 
space, and let §: E—>E be a continuous operator. If A is a closed, convex subset 


of E, ¥(A) CA and _.$(A) is compact, then there exists a fixed point of ¥ in A. 


Note that if % is also a compact operator and A is bounded, then .f(A) is 
necessarily compact. The following sets of functions are useful: 


(21) Hg, T) = \hE EF |All Sg (t), £= T SO, where T is a given constant 
and g is a given function, g€ Ef, g(t) +0 ast>o. 


(2.2) Sy, T) = {s€ EF |\lsOllSv(@),t=T =O, where T is a given constant 
and @ is a given function, p€ Ef. 


(2.3) &W&(m, T) = \acEt | lal Sm), t= T 20, where T is a given constant 
and m is a given function, m€ E>, m bounded on 


kgrco\t 


(2.4) WV'(c,y,T) ={a€ EF |\la@) —cl S y(t), t= T =0, where T,c€ R’ are given 
ue 8 
constants and p(t) is a given function, pC E>, p 
bounded on [T, ov). 


(85 Ad Use a i {SC EF ils —2@ Sy(t), ¢2T=O0, where T is a given 
constant, and J, y are given functions, /€ E4, pe Ef. 


Notice that if w(t) approaches zero as ¢—> ov, then in (2.4) each member of 
' approaches c as tov, and in (2.5) the difference between each member of 
Y' and the function /(t) approaches zero as t— ov. 

The sets defined above are closed, convex and bounded in the compact open 
topology. The fact that the sets are closed and convex is clearly true, and they 
are bounded since the definition of a bounded set only requires that the image 
of the set under each pseudo norm, #,,(x), is a bounded set on the real line. This 
does not require that the functions be bounded on [T, oo). 


Lemma II.1. Consider the equation 


(2.6) ) £4 = A,(a, h(t), 80,4) +Ao(a,h@, 8), ), 


where A,, A, are continuous on S={(a,h,s,t),aCUCR,hEVCR, se Re VenOs 
and U,V are open bounded spheres with centers at the origin of R’, R" respectively, 
and assume that 


(2.7) || Az (a, 2, s, é)|] S||Al]- Li (|Je|]) on 2, 

where L,(r) is a given continuous, positive, non-decreasing function for r = 0. 

(2.8) Fort20,a€U,hEV and se Ef, assume that ||A,(a, h, s, t)|| <K(t), where 
KEE} is given, K(t) 0 as too and f K(t) dt< oo. 


0 
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Tpige Eo is a given non-negative function, g(t)< vadius of V, +20, g(t) 0 
as tov, ft t)dt<oo; if ajCU is given; if C rest is any constant such that 


}20I ad oe radius of U; and if T,=0 is such that fel ce fi (u)) du, 
Qo||-+- 
PR) t) dt <C, then, for every T=T,, every s€ Ep. hCH(g, T), there exists a 


T; 
function m(t), mC EP, m(t) bounded for t=T, such that a solution a(t) of (2.6) 


satisfying a(t)= ay exists for t=T, a(t)€U, t2T, and a(t) les m A(m, T). 
Also there exists a function f(t), fC Ey, f(t) +0 as t—> co such that 


|| a(t.) — a(4)|| S74), eat aad © 


and hence there exists a constant c,€U such that a(t)—>c, as too. Also, the 
constant T,, and the functions m(t), f(t) depend only upon the functions K(t), g(t) 
and the constants a, and C. Furthermore, the bound on m(t) can be made as close 


to ||a9|| + C as desired by taking T sufficiently large. 

Proof: For any a€ U, let C> 0 as above be any constant such that ||a@ || + C 
is less than the radius of U. Let =O be such that D(Z, = SK) thdt<C 
and i [L,(u)]4du> fe g(u) du, where the integral of Ly' may be infinite, in which 


Ilao||+ 
case the inequality is Se The constant 7, exists by the hypotheses on K(t) 


and g(t). Consequently, there exists a solution m(t), mE EZ, with m(T)= 
||@o|| + C, T = Z,, to the equation 
d 
(2.9) Fr = 8 (t) Ly (m() 
for every T2T,. Clearly m(t)20, m(t)<|la)||+C+M(T), t=7, where 


M(T)—>0 as To. Now choose T = J, such that ||a,|| ++ C + M(T) < radius 
of U, and then a(t)€U for t=T, if aC. A(m, A; 


With this choice of T and m(t) we define I*: E4—>E% by the relation 
(2.10) I%(a) (2) = ay + f {Arfa(u),(u), s(u), 1] +A,[a(u),h(u), s(u), u}} du. 
From (2.7), (2.8) and (2.9) 

[C2 © Slag + Lee 1, (Ja) du f Kw) du 
S ||| +C4S gu Ly (1 (u)) du = m (8) 
for every sC ES, hE H(g,T), a€ sd (m, T). 


Thus, by TycHonov’s theorem, there exists a solution, a(t), of (2.6),a€ A(m, T) 
and by our choice of T, a(t)€U, for t= T. ; 
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Since a(t) = I%(a) (¢), it follows that 


lla (¢e) — a (&)| SJ Axfa(), WC), s(x), du + 


(2.11) | + fala le), 40), 50), wlll du S Ly (MA) F ge) due + 


1 


ae K(u) du = f(t) 
for tj =t, =T, where M, is a bound for m(t). Clearly {(t) 0 as ft > ov, and the 
lemma follows. 


Corollary 11.2. If, im addition to the hypotheses of Lemma II.1, the function 
A, 1s identically zero, then for every «>0, there exists a constant By(x) > 0 such 
that for all BS By(«) and hE# (Bae, 0), sc Ef, and a,€U, there exist a con- 
stant C, and a function m(t), mE E>, m(t) bounded for t=0, depending only on 
B,«, such that a solution a(t) of (2.6) with a(0) =a, exists for t=0, a(t)CU, 
t=0, aC A(m, 0) and a(t) satisfies the relation 


|| 2 (¢2) Pa(h) (pss Cs Bah, 2420; 


and hence there exists a constant c,€ U such that a(t)—>c, as t > 00, and a(t)>c, 
exponentially. Furthermore, the bound, M,(B), on m(t) ts given by ||aq|| + Cy B. 


Proof lu the ‘proot of ‘Lemma IT-1 let 7,= 7 = 0, €=.0,sand choose 98, (a) 
so small that «By («) < radius of V and for all BS By («), 


flo] dus f[pae“*du=BZ. 
I|-|| 0 

Then m/(t), the solution to (2.7) with m(0) = ||a ||, is defined and bounded for 
t=O, and, by further restricting f, the bound on m(t) may be made so close 
to ||a@9|| that a(¢) € U fort =0. Furthermore, f(é) in (2.41) is equal to L,(M)) Be™, 
and the conclusion follows. 


Remark II.3. It is clear from the proof of Lemma IJ.1 that, for an arbitrary 
aC A(m, T), hCA(g, T), sc EG, the function a,=J%(a), where I*(a) is defined 
by (2.10), satisfies the property ||a, (¢.) — 4, (4)|| S/(4), t2 24, = T, and f is given 
by (2.11); 7.e., a,(¢) approaches a constant c, as too and a4,€¥#'(c,,f, T), 
where .%’ is defined by (2.4). Furthermore, the function f depends only on 
T,g, K and a). Now consider the larger class of functions B ol iil le 


with the set U given in Lemma II.1, and U is the closure of U. If a€ B, the 
lim a(t) =c, exists, and the convergence of a function to its limit is uniform 
too 


with respect to a€ B. This uniformity implies that the map from B into k%, 

taking a into c,, is continuous on B, where the topology on B is that induced 

by the compact open topology of the function space E%. In fact, for a fixed 

a€ B, let a,€ B,n=1, 2,..., be such that a,—>a in the compact open topology. 

We wish to show that ¢,,—>c, as m—> oo. For any ¢ >0 choose 4,27 so large 

that f(t)<2e, and N so large that ||a,(f)—a(t)||<3e for »2N. Then, 
10* 
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for every n =N, we have 
l!eon— Call Ig — x (to)l] + |] (€o) — 2 (¢o)|] + I] (0) = al 
<2f(t) tte<fe+tge=e. 
Lemma II.4. Consider the system of equations 
(2.12) a2 (a(t), h(t), t) + Si (a(t), A(t), , t) + S, (a(t), A(t), s, 4), 


where S, S,, S, are continuous on X={(a,h,s,t), @€ UORt, EVR, sek. 
t=0} and U,V are open bounded spheres with centers at the origin of R* and R’, 
respectively. Suppose there exists a constant L such that 


(2.13) I|Si(@, 4, s, 4) || SLIIAl| 


for all (a, h,s,t)€ 2X, and also there exists a function K,(t), K,(t)€ E), Ky (t) 0 
as t—> oo, such that 
(2.14) || S2 (a, 2, s, t)|| SK (2) 


for all ac U, hE V, sCR®, t20. 

If gC EQ is a given non-negative function, g(t)< radius V, t20, g(t)>0 
as t->0o, then for any constant T=0, sC R?, hCH(g,T), and any function 
ac E} with a(t)€U for t=T, there exists a function pcEt, independent of the 
particular functions a, h, such that there exists a solution s(t) of (2.12), with s(T) =s9 
which belongs to S(y, T). 

Proof: Let I?: E%—+>E*% be defined by 


(2.15) 1*(s) (f) = 8 i [S(a(u), h(u), w) + Sy (a(u), h(t), s(u), «) + 
as 


u 
+ Sy(a(u),h(u), s(u), u)| du. 
Let S*(é) =sup(||S(a, A, t)||, a€ U, he Vv), where U, V denote the closure 
t 
of the sets U, V. Let p(t) =so+ J [S*(u) +Lg(w) + K,(u)|du. It is immediate 
T 


that for a, h as given, I?(S(y, T)) < F(y, T), and by the TycHonov fixed point 
theorem the conclusion follows. 

Remark II.5. In the proof of the preceding lemma the transformation J? is 
defined even if a(#) is a continuous function of s,, namely a(t) =a(s, h) (t) 
provided that, for sc AY(p, T), hE#(g, T), we have a€ Eb, a(t)€U fort=T. 
The proof above shows that I?(s)€ Y(g, T) if se (9, T) even in this more 
general situation for defined as above. 

By TycHonov’s theorem, for each hE W(g, T), there exists a function 5= 
S(A)ES(p, T) such that 5= J*(S), and this in tum implies that ds(t)/di= 
S(a(h, 8) (t), h(t), t) +S, (a(h, 3) (d), h(t), 5(d), t) + So(a(h, 8) (t), h(t), S(t), 2). 

Corollary I1.6. In addition to the hypotheses of Lemma II.4, suppose that 


ih 
yh 


dlyfed (u) du< Pel g(u) du < co, and also let f be a given non-negative function, 


[oe) 
fE Er, f()>0 as tc sisal {(u) du < oo. Further, assume that there exists a 
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constant L, such that for a,,a,CU, hy, hye V, t=T, 
(2.16) I| S(@1, Ma, t) — S(aq, he, t)|| S Ly (|]a, — a5|| + |] — All). 
Let F(t) be defined by 
(2.47) F(t) =f (Lif(w) + (L +11) g(u) + Ky(u)] dw. 
t 
For a€@'(c, f, T), hE#(g,T), se S (yp, T) there exists a og Sapte 4 C= 


Co(4,h, s) such that if l=I(a,h,s)C E% is defined by I(t) = Cot f SI (ciOntdd, 


t=T, then we have I*(s)€ F'(l(a,h,s),F, T), where I?(s) is eee by (2.45). 

Furthermore, Cy(a, h, s), l(a, h, s) are continuous in a,h,s foracC UM'(c, f, T) 
c€U 

A(m, T), he#(g,T), se S(p, T), where mt) is a given function, m(t)<M,< 

vadius U, t=0, M, constant. 


Proof: If s)€ R? is given, if 1(a, h) (t) =sy+ f [s(a (u), h(u), w) — S(c, 0, u)|du 
t T 
+ f S(c, 0, u) du, and if s,(a, h, s) (t) = I?(s) (t) —1(a, h) (2), then 
a8 


co 


SildeH,S, 1) Saal [S(c, 0, w) — S(a(u), h(a), u)|du + 


vr) [Si (a (uw), A(u), s(u), «) + S,(a(u), A(u), s(u), uw) | du. 
By (2.16), (2.13) and (2.14), for every 4;24,=2T7, we have 
|| si (4, 2, 8) (t2) — 8, (4, 4, s) (4)|| SF) 
where F(t) is defined in (2.17), and, therefore, lim S, (a, h, s)(t) = d(a,h, s) exists 
uniformly with respect to a€ Uw'(c, f, T)N Y(m, T), he#(g, T), se S(q, T). 
_ceU 
If l(a, h, s)(t) =d(a, h,s)+1(a, h)(t), then the preceding inequality is clearly 
satisfied for all ¢, >t, =T if the function s,(a, h, s)(t) is replaced by I?(s) (t) — 
I(a, h, s) (t) = s,(a, h, s) (t) —d(a, h, s). Therefore, letting t,—> 00, we have I?(s) € 
SF' (l(a, h,s),F, T). Furthermore, it is clear that 


l(a, h, s) (t) =C,(a, h,s) +J Ste 0,u) du, 
Cy(a,h,s) =ad(a,h,s) + so + F [S(a(), (4), w) — S(c, 0, u)| du, 


and the first part of the lemma is proved. 

We now show that /(a, h, s), Cy(a,h,s) are continuous in 4, h, s for ae e 
'(c,f,T) 1 dm, T), he#(g, T), se SF (y, £). By Remark IT.3, the ae 
taking a—>c,, where c, is the limit of a as tov, is continuous, and it follows 
that ten 0, uv) du is continuous. Therefore, Cy(a,h,s) will be continuous if 


T 
I(a,h,s) is continuous. On the other hand, /(q, h, s) (4) =4d(q, h, s) +1 (a, h), 


where d(a, h, s) = lim s, (4, h, s)(t), s,(a, h, s) (¢) = I?(s) (t) — i(a, h)(t), and their 
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limit exists uniformly over the above sets. Consequently, by using a proof sim- 
ilar to that in Remark II.3, one shows that d(a, h, s) is continuous if s, (a, h, s) is 
continuous. But since J?(s) is certainly continuous, it remains only to show that 


t oo 
L(a, h) = b(a, h) + f S(c, 0, u) du, b(a, h) = So + if [S(a(u), h(u), u) — S(c, 0, u)] du, 
T i t 
is continuous. But we have already remarked that f S(c, 0, #) du is continuous, 
2: 


and we now show that b(a, #) is continuous. We need to show that a,—a, h,—>h 

as n —> co in the compact open topology implies that b,=0(a,,h,) +o = b(a, h) 

as m—>oo. Let c,=lima,(t), and note from Remark II.3 that c,—>c as n> oo. 
t—0o 


Given ¢ > 0, choose ft) > T so large that f L,[f(u) +g(u)]< Ze, and choose N 
t 


so large that |la,(¢)—a(é)|| +||2, (4) —2(@I| < e/[4L1(o — T)], for Tata 
n=N, ||e,—e|| < e[[4L1(f— T)] for all »2N. From the definition of 6(a, h) 
above and (2.16), we have 


\|o,, — || < aS ee — S(a,h, w)|| du fl See 0, u) — S(c, 0, u)|| du + 
ge g% 


mr nN? 


+ fsa In, 4) — S(c,, 0, 4)||du + f || S(a,h, 4) — S(c,0, w)|| da 
to ty 


5 ITs (a) — a(u)|| + |]2, (%) — h(w)||) du + 
+fhille— c||du + 2b, Sf (u) +g(u)|du<e 


for all mn 2N, which completes the proof of Corollary I1.6. 


Corollary 11.7. In addition to the hypotheses of Lemma II.4, suppose that 
relation (2.16) holds, S,=0, f(t) =C,Be~*, a>0, and sete T=0. Then for 
hC#H (Bae, 0), a€X'(c, CoB ee, 0), s€ L(y, 0), the conclusions of Corol- 
lary II.6 are valid, with the addition that the function F(t) =L' Be’, for some 
constant L’. 


Proof: The proof is the same as that of Corollary I1.6 with K, (¢) =0, and 
£=0. The conclusions of Corollary II.6 are immediate, and F(t) is obtained 
by replacing /(u), g(w) and Ky (uw) by C,Be~*", Bae *", and O respectively in (247) 


Remark II.8. Notice in the proof of Corollary II.6 that the condition f(t) 
be integrable is not needed if the function S is independent of a; i.¢., a€ E%, 


a(i)CU, t2T, hE#(g,t), f g(u)du<oo and the Lipschitzian condition on 
Z 
S imply the conclusions of the Corollary. 
Lerma II.9. Consider the equation 
dh 
(2.18) = = H(a(t),s(), )h + A, (a(t), h, s(d), t) + H, (a(t), h, s(t), t) 
where H, H,, H, are continuous on r= t(d. Wes, by aie UCR  HEVC Ra seh 


t= 0}, and U,V are open bounded spheres with centers at the origin of R? and R’, 
respectively. Suppose that the functions H, H,, H 2 Satisfy the following properties 
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on Dd: 


(2.19) Given cCU, s€Es, let W(c,s, t) denote the matrix solution to dW/dt = 
A(c, s(t), t) W, with W(c,s, 0) =I, the identity matrix. Let us suppose 
that for some cCU, ICE}, B>0,0>0, we have ||W(c, 1, t) W>(c, 1, w)|| 
= Ber 0S at 


(2.20) There 1s a constant L,>0O such that ||H(a,, s,,t)—H(az, sz, t)||S< 
L,(|]a.— || + ||si—sel|) for all a, a,€U, 5, 526 R®, t=0. 


(2.21) Given e>0, there exists a d> 0 such that ||H,(a, h, s, t)|| Se||h|| for 
all a€U, sER*, t= 0 and |h|| SO. 


(2.22) There is a continuous function K,(t)C E}, K,(t) 20, t= 0, K,(t) +0 as 
t—>oo such that ||H,(a,h, s, t)||S K,(é) for all acCU, hEV, sER?, 
t= 0. 


If the above conditions are satisfied, then there exist constants T’=0, y’,7’, 
such that for every T=T', ywSy', nS1’, there is a function g(t)C Ep, g(t) >0 
as t->oo such that for every aC W'(c,p,T), cCU, seSF'(l,y,T) and every 
hy ER’ with \\ho|| <n, there is a solution i ) of (2.18), existing for TSt<+o, 


BAD) = =h,, and hE#(g, T). Furthermore, I Ryle )du< co implies ic ) du <o, 
LIK ) du dt-<oo implies f J g(w LW HEA cos 


Proofs leet-y,. > 0,6 > 0 2 numbers such that B(Lyy+ ¢«) <o. eS c= 


o—B(Lyy +6) >0. By (2.21), there is a 6 >0 such that || Ay (a(t), A, s(4), t)||S 
e||h|| for all a(t)€U, s€ Ef, t=0, |lAl| SO. pe ealiee: we ii Lice ft) 
so small that the sphere ieee isin V. By (2.22), ||H,(a(), h, s(), )|| SA 


for every a(t)CU, sCE§, hEV, t=0. 

Let 7'> 0 be any number, 7’<0, 7’ <6/B. Since K,(t) +0 as t > 00, we may 
take T’=0 so large that By’ + (B/«) sup (K,(t), #27’) <6. For a given c€U 
in (2.19) choose the constant y’< $y so small that if a€.V'(c, yy’, T) for T2T’, 
then a(t)CU for all ¢=T. In (2.5), for any py’ and /(é) in (2.19), define 
il, 1). Therefore, for every ae (ex Ves eer (UN, Mle) mace eae 
y <y’ we have ||H(a(i), s(t), t) — H(c, l(t), t)|| S Ly. Now consider the function 
gCE}, defined by the differential eqtaon 


(2.23) ae Te es 87) = By. an, 
1.¢., by the integral equation 


t 
(2:23). gij=se “Bn t fe’ BK (udu, t=T. 
a 


Obviously g(t) 20, g(T) =By and g(t) < By +sup (Ky (i) te i) (Bla)<6. It 
is also clear that g(t) +0 as ¢->oo. Another integral equation equivalent to 


(2.23) is 
t 
(2.23)” g(t) = Bye "9 + f Be“ [(Lyy + 8) g(u) + Ky (u)] du 
iB 
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Now let #(g, T) be the class of functions in (2.1) defined by the function 
g(t),t 2T, above. For every hy € V, ||7o|| S 7, 7 above, define the transformation 
I’: E’-E’ as 

t 
0.24) (A) () — WW WA(T) hg +f WO) WA (w) [{H (aw), (0), 4) — 
T 
— H(c,l(u), u)} h(w) + A, (a(u), (uv), s(u), u) + 
+ H,(a(u), h(u), s(u), «)|du, 


where W(t) =W(c, lJ, t) is defined in (2.19). For every hEH, we have 
t 
(2.25) |lI"@) Os Bye? + f Be [ay + eeu) + Ki) de. 
T 


From (2.23) the second member is g(#), and hence ||J’(%) (‘)||Sg(), t2T. 
Therefore, I’(h) (t)€H(g, T) and by TycHonov’s fixed point theorem, there 
exists a solution h(t) of (2.18) belonging to #. To complete the proof, it is suf- 
ficient to recall that g(t) +0 as t->oo, and to observe, from (2.23), that 


(2.26) J g(u)du=o[g(T) —g)] + Bat f Ky(w du. 


Remark I1.10. Suppose all the conditions of Lemma II.9 are satisfied except 
that in (2.19) the bound on W(c, J, t) W4(c,/, uv), OS wu St, is supposed to exist 
for any cCU,1€ J CE§ for some set J, where B is assumed to be independent 
of c,/ and o=a(c)>0 is independent of / and is a continuous function of c. 
Given T,y and a function g€ E7, g(t) +0 as ¢t->oe, consider functions a= 
a(h,s)E'(c(h, s), yp, T) N L(m, T), s=s(h)E F'(L(A), y, T)N S(y, T), where 
a(h,s), c(h, s), s(h), l(h)€ J are continuous for hEA(g, T), se (p, T) and 
\|c(A, s)|| Sd < radius of U, for hE#(g, T), s© Y(p, T), where d is a constant. 
Then there exist constants JT’ 2 0, y’, yn’ such that for all T=T’, pSy’, nSX7’, 
there is a function g€ Ey, g(t) +0 as ¢— oo such that for every hy € R’, |Zo|| Sn 
and any functions a(h, s), c(h, s), s(h), and J(h) as above, there is a function 
hE Hg, T), h(T) =hy such that I’(h) =h. This clearly implies 


= H(alh, s(h) ), sh) (0,1) A()) +H, (a(h, s()) (1), 2(2), 8 (A) (8), 2) + 
+ H,(a(h, s(h)) (t), A(t), 8(A) (t), 2). 


The proof of this follows from the proof of Lemma I1.9, by first letting ¢o = 
mina(c), ||c|| Sd. If « is defined using this value of o, then the function g(t) 
in (2.23) is independent of a, s, I, c, |c|| Sd. Further, W(c,J, 2), as a solution 
of the equation, dW/dt = H (c(h, s(h)), L(h) (t), t) W, is also a continuous function 
of AEH (g, T) (see Lemma II.14 below). Thus the operator J’ is well defined, 


and the proof asserts that I’(#(g, T)) (#(g, T). The result follows by Tycuo- 
NOv’s theorem. 


This will be needed in Section 3. This same remark applies to Corollary II.2 


and Lemma II.13 below. For simplicity, these lemmas are not stated in this 
more complicated form. 
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Remark II.11,. Let us analyze some of the implications of the statement 
and proof of the above lemma. First, if a, s do not appear in (2.18), then Lemma 
II.9 coincides with a theorem of E. A. Coppincton & N. Levinson [5, p. 327]. 

Second, in the proof of Lemma II.9, the only restriction on the constant T 
for the given constants B,7',«,6, was that By’ + (Bla) sup (K,(d),t=T") <6. 
However, if we suppose that the perturbations are of a more restricted character, 
namely, that the function K,(t), not necessarily approaching zero as ¢ > oo, can 
be made as small as desired for all t= 0, then the above inequality can be satisfied 
for all? 20, and so the function g(t) in (2.23), for T = 0, is well defined and < 6 
provided K, (é) is sufficiently small. Then the remainder of the proof of the lmma 
is exactly the same in order to show that there is a solution h(t) of (2.18), with 
h(0) = ho, ||%ol| Sn and ||h()||S g(t), t= 0. 

Notice that the equation (2.18) has the solution h=0 for H,=0 and the 
solution h=0O is asymptotically stable. However, for H,=0, the function 
h =O is not generally a solution of (2.18), but with the restricted class of per- 
turbations above, the ¢-axis in (h, t)-space is stable in the sense of LYAPUNOV. 
This contains a result of MALkIN [16, p. 302]. If we further require that the 
function K,(t) +0 as foo and that it also can be made as small as desired 
for ¢= 0, then the ¢-axis is stable under this type of perturbation, and h(t) +0 
as [> oo, 


Corollary 11.12. In addition to the hypotheses of LemmaII.9 assume that 
the function H, is identically zero. If y>0, e>0 are any numbers such that 
Bl,yy+te)<o and a=o—B(Lyy-+e), then there exist constants wy’ and 
"=O, such that, for psy’, ns, scFiLy,0), a€V(c,y,0), Mek, 
I|2o|| Sn, there is a solution h(t) of (2.18), h(0) =I, and hCH(Bne™, 0). Thus 
ASB nen, 1=0. 

Proof: It is clear that T may be chosen equal to zero in this case. Further- 
more, in (2.23), K,(t) may be chosen equal to zero, and the conclusion follows. 


Lemma II.13. Suppose that system (2.18) satisfies all of the conditions of 
Lemma II.9 except that (2.19) ts replaced by the following stronger hypothesis: 


(2.27) For c€U, there exists B>0, o>0, so that for every sC€Es, we have 
\|Wic, s, 4) WA(c, s, w)|| S Be) 0S ut, where Wc, s, t) is defined 
im (2.19). 


Then there exist constants T’ = 0, y'>0, 7'>0 such that for T2T’, yn S71’, 
psy’, there is a function gCEr, g(t) +0 as t—+co such that the conclusions of 
Lemma II.9 are valid with the set S'(l, p, T) replaced by EF. 


Proof: The proof is almost identical with the proof of the preceding lemma 
except that a new operator Jj: E’—> E” is defined by 


Ij (h) (t) = W(t) W1(T) ho +f W(t) Ww) [{H (a(u), s(u), 4) — 
(2.28) —H(c,s(u), u)}h(u) + Hy (a(u), h(u), s(u), u) + 
+ H,(a(u),h(u),s(w), u)|du, 
where W(t) =W(c, s, t) is defined in (2.19). 
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The following lemma is basic for the proofs of all of the theorems in Section 3. 
Lemma II.14. Let f,(x, v, z, 4, 4), fo(%, ¥, 2, 4, t) be continuous on R?x RIX 
R'XRXR, with f, having values in R?, fy in R1. On thas domain, assume for 


eats 2, 
: IFC Vas Za Mb) — K(X Vos 22, % 2) 


<L;(u,6, 0, 4, t) 2 x. || -F I| 4 = yall at l|z1 — 22 


ie 


where = max ||x,||, o = max ||y;||, @ = max I|z, ||, 7= 1,2, and L; 1s continuous 
in all of tts arguments. 
Let y, : R->R? be continuous, and define I, : Ef x E§ x Ey > EG by I, (x,y, 2) (t) 


t 
=) +S A(*(u), v(m), 2(w), 4, t) du. Let A,, Ay, A, be bounded subsets of Ef, 
0 


E% and E}, respectively, and assume that for (y,2)€A,XA,, there exists at least 
one x =x(y, 2)C A, such that I, (wig) 2) Ve 2) SAN ae 
Let y,: RR! be continuous, and define I,: Ej x Eo > EG by Ip (y, 2) () = 


gall) +f fol=(o, 2) (4) 9), (0), af) ds 


Assume that for z€ A, there exists at least one y = y(z) CA, such that I, (y(z), 2) 
= y(2). 

Under these conditions, the function x(y, 2) 1s well defined (1.e., unique) and 
continuous on A,X A,, and y(z) 1s well defined (t.e., unique) and continuous on Aj. 


Remark II.15. The operators I, and J, above are defined in this general form 
so as to include types of integral operators ordinarily encountered in differential 
equations, that is, those of type J? or J’, defined in (2.15) and (2.10) respectively, 
or the operators J” or Jj, defined in (2.24) and (2.28). 

Proof of Lemma II.14: For n a positive integer, let ||x||,=sup||x(4|I, 
0575, and for y,,%,¢A,, 4, 22€4,, let %—%(y,, G), Ha—* Va naa), Delany. 
functions such that Ty (%;, V3, 2;) = x;,]7 = 1,2. We shall show that there exists 
a constant C(m) such that ||x,— %9||,< C(x) [|]v:— yell, + |]a—%ll,]- For 
V1= Vo, %=%, it follows that x;—~*,, 1.e., the fixed point is unique. The same 
indequality clearly implies that this function x(y, z) is continuous. 

Let K;=sup L,(u, 6, 0, u, t) for OS u, 0,0 SB, u,tC[0,n] where B is the 
bound over [0, | on the elements of A,, A,, A,. B exists, as all of these sets 
are bounded. 


Since %,, %) are fixed points of J,, it follows that 
(2.29) ||“ — »||() S Sa [Il 1 — all (#) + [] 1 — vel] (#) + lla — 20] ()] ae, 
for all ¢€ [0, n]; where || %,— xg|| (#) =|| x(t) — x9(¢)||, etc. 

If r(Z) i, || %1— x9|| (4) du, g(t) ae [|| 41 —ve|| (2) + |lai—zel]()] du, then 


Bad ant) <= K,g(t). Since 7(0)=0, we have r(t)< f ell) Ko (u) du 
< K,e™" f g(u) du, or : 
6 


(2.30) J || 1 — *9|| (w) du < K, KF Ella — Yell (t) + |]. — 29]| (t)|dtdu. 
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From (2.29) and (2. oy it is now clear that the desired inequality is readily 
obtained, with C(m) = Kin 4+ Ki n®e%", 

Now take %, 2,€A,, and let y,;=y(z,), x, ;=%(y;, 2;), be any functions such 
that I,(y;, 2;)=4;,7=1, 2. For a positive diiterer it suffices by the previous 
remarks to show that there exists a constant C’(m) such that ||y,— yell, < 
C’ (n) || 2 — 2||,- 

As before, for ¢ € [0, 2] and using (2.30), we obtain 


Il 4 oi Yol| (t)S os) [| «1 = X»|| (u) + ll ae Yell (u) + || 41 - 2|| (1) | du 
(2.31) KJ [lll (u) + ||41 — 2a] (#)] du + 


+ K, Kye f fll — yell (© + lla —all(o)) dean. 


tu 
Now let 4,(t) Py ieee Fede gandh are oe pil K, ae 


— K, Kye", (t), h(t) =K 2 lla all H) du + Ky Kye f flia,— all (0 dt dau. 


The above inequality is now equivalent to k(t) < A(t) for ¢€ [0, ]. To estimate 
p(t), ap, (t)/dt = p,(t), consider the following vector differential equation: 


P — Ap) +k), ~(0) =0, 
AS 


where p =col(py, po), (4,4); Ga On G1, dey — BR Kes, Gos Ka, 


k =col(0, k). Now, 
(2.32) p(t) =f AKMR(u) du. 
0 


Since a; ; 2 0, it follows that every element of the matrix e“' is also non-negative 
for £0. Thus, letting 4 = col(0, h), we see that (2.32) implies 


(2.33) p(t) <f (uw) au, 
0 


since k <h, where the order ‘‘<” is vector ordering, 7.e., x <y if and only if 
%, Sy, for all 1. 
So, if A is a bound on every element of e#’, for te [0, »], (2.33) and (2.34) 


imply that ||¥4— ell (2) h(t) + A (K,Kye%"-+ K,) fixe ) du, and since h(t) < 
(Kan + K, K,n?e%") eae, the desired baeadlity can be obtained easily. 
Remark II.16. Note that if the Lipschitz constants L;(u,o, 0, u, t) do not 
depend upon one of the arguments yu, 0, @, then the corresponding set A,, A,, A, 
need not be bounded. 
Remark II.17. The inversion of the inequality used in the proof of the first 


part of Lemma II.14 is GRONWALL’s inequality, and the proof given here is the 
standard one (see, e.g., [2, p. 35]). The proof of the second part of Lemma IT.14 
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is a generalization of GRONWALL’s inequality. The same result could also be 
obtained by integrating by parts twice in the right-hand side (2.34) and then 
applying GRONWALL’s inequality. However, the above proof is very easily gener- 
alized to the case where more than three integral equations are solved simul- 


taneously. 


III. Stability theorems for integral manifolds 


In this section we consider systems of equations of the form 


qs _ S(a,h,t) + S,(a,h,s,t), 


ae 
da 

(5.4) di = A,(4,h,s,1), 
oh = H(a,s,t)h + H,(a,h,s,2), t20, 


and of the form 
4s _ Sia, h,t) 4 S,(a,h,s,t) + S,(a,h,s,t), 


dt 
(3.2) 29 == Ay(a,h,s,1) + Ao(a,h,s,2), 
GY =H(a,s,t)h + Hy(a,h,s,t) +H(a,h,s,f), t20, 


where s€R?, aC R’, hCR’, and the functions S, S,,A;, H, H;,] =, 2, satisfy 
certain conditions which will be stated precisely in the theorems below. One 
condition is that S$, =A,=—H,=0 for h=0, and, therefore, system (3.1) has a 
family of solutions h=0, a=c, ds/dt =S(c, 0, t), s(0) =s ), where cC R’, sy€ R? 
are arbitrary constants. Theorem III.1 below is concerned with the stability 
properties of this family of solutions of (3.1). System (3.2) is similar to (3.1) 
except for the “‘perturbation” terms S,, A,, H,, which are not necessarily zero 
for h=0. A natural question to ask is the following: What type of ‘‘perturbation”’ 
terms in (3.2) yield a family of solutions of (3.2) which has essentially the same 
characteristics as the family mentioned above for the unperturbed equation 
(3.1)? Theorem III.4 gives a partial answer to this question. In case the a 
vector is absent in (3.1), (3.2), then the properties of the systems are somewhat 
different, and this is discussed in Theorems III.7 and III.9. 


As mentioned in Section I, applications of these results are given in Section V. 
The proofs of the following results will be given after Theorem III.9. 


Theorem III.1. Consider the equations 


a 
s = S(a,h,t) + S,(a,h,s, t), 
d 
(3.3) <F = Ax (a, h, s, 1), 
dh 


nh H(a,s,t)h +H, (a, h,s,t), 


where s 1s a p-vector, a is a q-vector, and h is an r-vector. Let Uc R!, and VCR’ 
be bounded open spheres with center at the origin, and assume that the functions in 
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(3.3) satisfy the following properties on Y =U x V x R? x [0, ov): 


(3.4) In &, S, S,, Ay, H, and H, are continuous in t and Lipschitzian in a, h, s 
with Lipschitz constant L ; 


(3.5) 37'(4,10)-5, 6) = Agia, O,S,t) = 0, for a@eU, SCR, te 6, \ and ihus 
|S. (2,2, s,8)|| SLIAl|, ||Ar(a, 4, s, t)|| SS L||h||, where L is the Lip- 
schitz constant given in (3.4); 


(3.6) Given > 0, there exists a constant 6 > 0 such that || H, (a, h, s, t)|| Se||h||, 
for ||h||S 6, a€U, sER?, t=0; 


(3.7) Let c€U be constant, and 1C E$ be the solution of alldt—S(C 0, 1), 
1(0) =, for t=0. If W(t) is the matrix solution of the equation dW/dt 
=H (c, l(t), t)W, W(0) = the identity matrix, then there exist constants 
B>0,6(c) >0 for allc€ U, o(c) continuous inc, such that ||W(t) W(u)]|| 
<Be OC), for t= uz= 0, and all l,€ R?. 


Under the above assumptions, given a,€ U, there exist constants « >0, n’>0 
such that for all n<1', ssER?, agER’, |\ao||S||a||, and hoER’, ||Ao|| Sn, the 
solution (s(t), a(t), h(t)) of (3.3) through (sy, a), A) at t=0 has the following prop- 
erties: 

(3.8) FOllSBne™, t20; 


(3.9) lim a(t) =c,€U exists and |\a(t)—c,||Sc'ne-™, for some c'>0 and 
too 
alt=0; 


(3.10) ||s(4) —Z (|| SL’ ne~™ for some L'>0 and all t= 0, where I(t) is the 
solution of dl/dt=S(c,,0,t) with 1(0)=Cy, for t20 and some Cy= 
Ca lash,s). 

Remark III.2. A special case of Theorem III.1 has been considered by 
1. G. MALKIN [16] where the variable s in (3.3) is absent and all of the functions 
A,(a, h, t), H(a, t), H,(a, h, t) are analytic in a, h in a neighborhood of the origin 
in R’x R’. MALKIN obtains the same result as above for this special case, namely, 
in a neighborhood U of the origin in R’x R’, every solution [a(é), A(t)] with 
[a(0),2(0)]€ U has the property that h(t)->0 exponentially and a(t) c, ex- 
ponentially as t > oo. 

Remark III.3. Theorem III.1 as well as the theorem of MALKIN mentioned 
above has an interesting relationship to some results of K. P. Persrpsxr [18] 
(see also S. LEFscHETZ [11, p. 117)). 

Consider the system of equations 


a) dy/dt=Y(y, 2,2), 
O44) b) ‘dejdt =Z(y, 2, t), 
where vER!, zE€R’, Y(0, 0, t) =Z(0, 0,4) =0 for t20, and Y,Z are continu- 


ous and Lipschitzian in some neighborhood U of (0, 0)€R’x Rk’. Let C(t) bea 
continuous v-vector, and consider the equation 


(3.12) dyldt = Y(y, C(t), t). 
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According to PERSIDSKI, y = 0 is quasi-stable [quasi-unstable] for (3.11a) when- 
ever the following condition is [is not] fulfilled: Given e> 0 and t) = 0, there 
exists an 7/(é, t), 0< < « such that if ¢(¢) is any function such that IF (40) <1 
and ||C(#)||< ¢, #=t, then any solution y(¢) of (3.12) with || v (49)|| S 7 satisfies 
ly (Q)||<e,#2%. By symmetry, the same notions are defined for (3.11b). PER- 
SIDSKI then states that if the origin is quasi-stable for both (3.11a) and (3.11), 
then the solution y = 0, z= 0 of (3.11) is stable in the sense of Lyapunov. Also, 
if the origin is quasi-unstable for one of the two partial systems (3.114), (3.11 b), 
then the solution y = 0, z= 0 of (3.41) is unstable. The last part of this statement 
is false as the following example shows. 


Consider the equation 


(3.13) daldt =ah, dh|dt=—h, 


where a, h are scalars. This system can be easily integrated, or one can observe 
that it is a special case of the system (3.3) without the s variable and, therefore, 
there exists a neighborhood U of (0,0) such that the solution [a(¢), h(¢)] of 
(3.13) with [a(0), h(0)]€U implies h(t) 0, a(t) a constant ast — oo, and the 
constant may be made as close to a(0) as desired by taking (0) sufficiently 
small. Therefore, the solution a=0, h=O of (3.13) is stable in the sense of 
Lyapunov. But the solution a= 0 of the equation da/dt = a*h is quasi-unstable 
since there are solutions of da/dt=a?y approaching oo in a finite ¢ range for 
every 7 > 0. 

This example points out how difficult it is to determine conditions for the 
instability of a system of equations by examining a partial system. For some 
results along this line, see the paper of M. R. DyHMan [8a]. 


Theorem III.4. Consider the equations 


d 
<= S(a,h,t) + S,(a,h, s,t) + Se(a,h,s, 2), 
da 

(3.14) de = A1(4, 4, 5, t) + Ag (a, h,s, t), 
adh 


Gq =A s,t)h + H,(a,h,s,t) + H,(a,h,s,t), 


where S, S,,A,, and H, are as in Theorem III.1, and suppose (3.7) ts valid for 
all JE ES, c€U and B,o(c) are independent of 1CE?, with o continuous on U. 
If X=U x V x R? x [0, oc), U, V as in Theorem III.1, then assume that S,, Ay, Hy 
are continuous on 2' and Lipschitzian in (a,h,s), and there exists a continuous 
function K,(t)—>0 as t->0e such that, on &, 


(3.15) Il S2(4, 2, s, 4)|| SK, (0), I|42(4,h,s,t)|| < Ky (0), || Z2(a, h, s, t)|| S Ky(t). 


If, in LH) K, (t) dt < 00, then given a,€ U, there exisis ny’ >Oand T'=0, 
such that if T2T',n<n' and (s(2), a(t), h(t)) is the solution to (3.14) through 
(Ss, Aq, hg) a t=T where SpE R?, agER’, ||ap]| <|]a,|], AER", ||ho|| <n, then 
S(t) 1s defined for t= T, jim. a(t) =c,€U exists and h(t)>0 as t >. 
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Moreover, if fi if K,(u) dudt <0, then given a,CU, there exist constants 


Ha On, Lee) stay that if T2T', nS’, and (s(t), a(t), h(t)) is the solution to 
(3.14) through (S89, a),h 9) at t= » MER’, 
|2o|| <1, then h(t) > 0, lim Ae Bue U exists, and ||s(t) —1(t)|| +0 as too, 
where dijdt =S(c,, 0, t), 1(T) =Cy, t2T, Cy=C,(a, h, s). 

ree i It will be clear from the proof of this theorem and Remark II.8 
that ||s (é Bilal as t-> ce, with s,/ as above, if S(a, h, t) does not depend 


on a, provided re K, (t) dt < 00, even though f i K,(u) du dt may not be finite. 


Also, the last somelliien of Theorem III.4 ~ salle if we require that condition 
(3.7) is true only for those /(¢) for which dl/dt = S(c, 0, t), c€U. 


Remark III.6. Without the variable s in (3.14), some closely related results 
have been obtained by M. R. DyuMan [80]. 


Theorem III.7. Consider the equations 


4 _ Sth, t) + Sy(h, 8,0) + Se(h, 5,2), 
(3-16) hs 
a7 =H (s,t)h+ Ay (h,s,t) + Hy (h,s,2), 


where S, S,,S,,H, H,, and H, satisfy the conditions of Theorem III.4 with the 
variable a absent. 

Then there exists aT’ =O and an y'>0 such that for every T=T’, " =e 
the solution (s(t), h(t)) of (3.16) through (so, ho) at t=T, where syCR?, ER’, 
\|Ao|| Sy, ts such that s(t) is defined for t=T, and h(t) >0 as t>o. 


Remark III.8. This theorem is quite elementary and can be proved very 
easily by other methods than those used here. The following theorem is of a 
more difficult nature. 


Theorem III.9. Consider the system of equations (3.16) where all of the 
hypotheses of Theorem III.7 are satisfied except that rt 1s required that the matrix 
solution pa of dW/dt = H (I(t), t) W, with W(0) =I, the identity matrix, satisfies 
|W) W3(u)||S Be", OS u St, < only for those I(t) for which dl/dt = 


S(0,t) fort =T, and B, o independent of |. If JR) t)dt <oo, then there exista T’=0 


and an y'>0, such that for every T=T’, Pe <7’, the solution (s(t), h(t)) of (3.16) 
through (Sy, Mo) at t=T, where sy R?, hyER’, ||hol| Sn, ts such that h(t)>0 
as t > 00, and ||s(t) —1(t)||> 0 as t >, where L(t) is a solution of dl/dt = S (0, t), 
Pie Catt = 1. CoC, is, 1). 

Proof of Theorem III.1: Let o = inf o(c ) Mls < d, for a fixed d, ||a,||<d< 
radius of U. Using this value of o, let « =a— B(Ly + e) as in Corollary I1.12. 


(i) Consider the differential equation 


(3.17) de = Ala, h(t), s(t],  4(0) = 49, 
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where a,€U is a given constant, ||a)||S ||,||, s(é) is an arbitrary function in 
E’, h(t) is an arbitrary function in H (Bae, 0), B > 0, where # is defined in 
(2.1) and « given above. The conditions of Corollary II.2 are clearly satisfied, 
and thus there exists a constant B, > 0 so that for any B S fp, there exist a con- 
stant C, (in this case, C,=L) and a non-negative function m€ E>, m(t) S M,(B), 
t>0, M,(B) =||4|| + C28 < radius of U, such that a solution a(/) = a(h, s) (t) 
of (3.17), a(0)=4 , satisfies the properties a(h, s)EA(m, 0), |\a(h, s) (te) — 
a(h, s) (t)|| S CoB eth t= 2.0, for alluse Ee, HeZ (pu e7%'£0). yy bhat ie 
there exists c(h, s)€U, ||c(A, s)||S M,(B), such that a(t) >c(h,s) as t>o, 
a€ A’ (c(h, s), CoB e~*, 0), for hE# (Pa e 0), s€E%. Furthermore, from 
Lemma 11.14 and Remark II.16, a(h,s) (¢) is a continuous function of (h,s) for 
hC#H, sCE’, and so, by Remark II.3, c(/, s) is also a continuous function of 
(A, s) on these sets. 


(ii) Now, by Remark II.5, there exists a function gC E> such that for any 
HEH (Bae, 0), BX Bo, there is an s=s(h)€ (gy, 0) such that for t= 0, 


b.18) S!= sfa(h,s) ), Al), + Sila(h,s) OA), s@,41,  s(0) =50, 


where a(h, s) (f) is the function determined above satisfying (3.17). Furthermore, 
since S is Lipschitzian in a and h, it follows from Corollary II.7 that there exist 
a Cy=C,y(a(h, s(h)), h, s(h)) =Co(h) CR? and an l=1(a(h, s(h)), h, s(h)) = L(A), 
with dl(h)/dt=S(c(h, s(h)), 0,t)), L(h)(0) =Cy(h) such that s(h)C ¥'(1(A), 
L' Be, 0), since s(t) = I?(s(h)). Also, from Lemma II.14 and Corollary II.7, 
the functions s(f), Co(h) and /(h) are continuous in h, hE#(B ae, 0). 


(iii) Finally, we apply Corollary 11.12 and Remark II.10 as follows: In 
Remark I1.10, let J be the set of functions / for which the estimate in (3.7) is 
valid; that is, dl/dt=S(c, 0,2), c€U, 1(0)€R? is arbitrary. By (ii) above, 
clearly 1(h)€J, for he#(Bae, 0), BS By. If we further restrict B, so that 
My (Bo) = ||4o|| + C2bo<4, then from (i), ||c(h,s(h))|| Sd for BSB, and 
hCH (Bae, 0). If y’, yn’ are the numbers given in Corollary II.12 and Re- 
mark II.10, then we further restrict By in such a way that for B= By, we have 
C,psy, LB sy’, anda p/B <n’. Thus, if y, 7 are defined by y = max [C, p, L’B], 
7 = «B/B, then the conclusions of Corollary II.12 and Remark II.10 state that 
there is a function h(Z), h(0) =hy, he#(B ae, 0), for ||Ao|| <7 =« B/B, and 
h satisfies, for t = 0, 


G19) “70 = Hath, s@) (), 9), 1hW) + Hlalk, s) ,k(), 5 (8), 4), 


where a(h, s(h)) and s(h) are the functions defined above. The remainder of 
the proof of the theorem is now obvious since the functions 3 =s(h),@=a(h, 3) 
and A are a solution of (3.3) with the desired properties. The Lipschitzian bon 
ditions imply the uniqueness of the solution, and the theorem is proved. 


Proof of Theorem 14: Let T20, n> 0 be constants to be determined. 
Let o =info(c), ||c||<d, for a fixed d, |< d=-radius of Us let oe 
B(Ly + €), where y > 0, ¢ > 0 are any numbers such that « =>0. Let git) be 
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the solution of (2.23), dg/dt =— ag +BK,(t), g(T) =By. One can show easily, 


as was done in the proof of Lemma II.9, that g(t) +0 as ¢—> oo; that f K,(t) dt 
T 


< co implies f g(t) dt< oo; and that f f K,(u)dudt< oo implies [ [ g(u) dudt 
T ES Tj 


< oo. Also, since g(t) satisfies the above differential equation, it follows as in 
the proof of Lemma II.9 that g(t) <By-+ sup (Ky (é),t=T) (B/a), and, thus, 
there exist 7) 20, y)>0, such that for all T= Tj, 7 <p, we have g(t) < radius 
of V for i2T. Using T=, 7X7 and g(t) given above, define the class 
wt (g, L) by (2.4). 

(i) First, consider the equation 


(3.202) § “* = A,[a, h(t), 8,4) + Agla,h), s(0),4,  a(T) =a, 


where a,€U is a given constant, ||a9|| <||a,||, s¢ Zf is arbitrary, and / is an 
arbitrary function belonging to #(g, T). Since L, (7) in Lemma II.1 is a constant 


in the present situation, we have f [L,(u)]4du—=-+ for all C, and the second 
lleo|+O 
inequality on g in Lemma II.1 is satisfied. Take C >0, so that ||a,|| + C <d, 


and further restrict Tj so that f K, (t) dt <C forall T => 7. Thus, from LemmalI1.1, 
T 


there exist functions m, € EZ, m bounded for ¢ = T, / (¢) > 0 as t > , such that 
a solution a(t) = a(h, s)(t) of (3.20) exists for¢=T, a(t)€U, t=T, a€ A(m, T) 
and ||a(t.) —a(t)||=/(G), 42427. By further restricting J), m can be 
chosen so that m (t) < d for allt = T = 7T,. Now, by Lemma II.14 and Remark II.16, 
it follows that a(h, s) is continuous on W(g, T) x E%. In addition, a(h, s) (t)> 
c(h, s) as t->oo, uniformly on #(g, T) x EF, since ||a(h, s) (t) — c(h, s)|| S(O, 
t=T=T), and f(t) given in (2.11) is independent of the functions a, h, s over 
the above sets. Thus, by Remark II.3, c(h, s) is also continuous on HW(g, T) x E>. 

(ii) Applying Lemma II.4 and Remark II.5, one can choose a function 
@CEy such that for every hC#(g, T) there exists a function s(¢) = s(h) (é) 
waithis(/)==sy, se-/(, 7) ands satisiies, for ¢2 7 = 1,, 


ah oY) = Sfa(h, s) (t), h(t), t] + S,[a(h, s) (t), A(d), s(6), 4) + 
| + S,[a(h, s) (4), A(t), s(é), 4). 


Furthermore, g is independent of the particular function hE #(g, 7). 

(iii) Now, we apply Lemma II.13 and Remark I1.10 as follows: In Remark I1.10, 
let J] =E%. If y’ is the number given in Lemma II.13 and Remark I1.10, then 
we further restrict J) so that sup (f(t) <y’. Then if p= any (7 (t)), there exist 

t=T, t= 


‘<n, T’=Ty such that if g is chosen as indicated in (i) above wane Ih 2 
n <n’, g(T) =By, then there exists a function hc#(g, T), h(T) =ho, |%o||S 7, 
and h(t) satisfies, for t= T, 


dhi(t) s - = a a ee 
3.22) dt H[a(h,s(h)) (t), (2) (), a) + [a (2, 


”n 
poe 
(he el 
= 
— 
™/. 
SS 
>| 
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~~ 
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=< 
SS 
Js 
= 
Se 
SSS 
™~ 
(Wars 2! 


Arch. Rational Mech, Anal., Vol. 6 14 


154 J. K. Hate & A. P. STOKEs: 


where a(h, s()), s(Z), are given in (i) and (ii). The remainder of the proof of 
the first part of Theorem III.4 is now clear. 


CO CO 
To prove the second part of the Theorem, assume now that f f K,(u) dudt 
Cores 


< oo. This implies that f f g(w) du dt < oo, and, therefore, f f(t) dt < ce, from 
ras ze 


the definition of f(t) in equation (2.11). Thus, by Corollary II.6, there exist 
Cy(a(h, s(h)), h) = Co(h), U(a(h, s(h)), h) = L(A), with dl(h)|dt = S (c(h, s(h)), 0, ¢), 
1(T) = C,(h), such that s(h)€ -Y'(l(h), F, T), where F CE} is defined by equation 
(2.17) and F(t) > 0 as t+ oo. Thus, s(h) (t) —1(h) (f) + 0 as t + co, and the proof 
is complete. 

Proof of Theorem III.7: The proof of this theorem is almost obvious since 
it follows so closely the proof of Theorem III.4. Since the variable, a, is absent 
in (3.16), one need not repeat step (i) in the preceding proof. 

Proof of Theorem III.9: This follows by applying the same reasoning as in 
the proof of Theorem III.4 except omitting (i) and applying Lemma II.4, Corol- 
lary II.6 and Remark II.8 in (ii), and Lemma II.9 with Remark II.10 in (iii). 


IV. A convenient transformation of variables 


Let U? be the #-dimensional solid sphere in R? with center at the origin 
and radius @, and let U? denote an open set in R?. 


Consider the system of differential equations 


(4.1) o% = X(x) + X* (t, x) 


where x, X, X* are (n + k + 1)-dimensional vectors, and suppose that the auto- 
nomous system 


X€C(U"***1) has a (k + 1)-parameter family of periodic solutions 


x = x°[w(b) t+, 5], 


(4.3) XO (Stet, 0) == aoi(s 0). 


where b = (by, ..., x), p are constants and w(b)>0. Naturally x°(s,b) as a 
function of s belongs to C3. We assume that as a function of 0 it belongs to C2 
and also that w(b) € C2(U*)* , 

Furthermore, we assume that, for all }€ U*, sCR, 


0 *9(s, b) 0 x9(s, b) 
os” > ab =k+14. 


(4.4) rank 


ee (4.4) is equivalent to the statement that the function @x° (s, b)/s, 
[Ox(s, b)/Os] [Ow (b)/0b,] t+ Ax%(s, 6)/0b;, 7 =1, 2,...,k, are linearly indepen- 
dent, and, furthermore, each of these functions is a solution of the linear varia- 


* a URABE [22] has shown that these conditions are satisfied provided that 
 (b) is bounded in any Closed subset of U*. He also introduced a transformation of 
coordinates similar to the one discussed in this section. 
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tional equations 


(4.5) = hy sro) i+, by) 


for every fixed value of b. Therefore, there are k +4 chracteristic multipliers 
of (4.5) equal to one; or, k+1 characteristic exponents may be chosen equal 
to zero. Suppose that the remaining » characteristic multipliers are different 
from one, and denote the corresponding characteristic exponents by «,(2) 
a, (d). 

If s=q@(b) t+ g in (4.5), then 


ypeeey 


Gyan ant 0X [#°(s, b)] 

(4.6) Ue” Ay Ox > 
and if we let 

0*°(s,b 0m(b) [| s— 0x(s,b 
(4.7) Po (s) a a p; (s) me Pols) ae aad o > 
7 =1, 2,..., k, then it is clear that there is a fundamental system of solutions 
Y(s) of (4.6) of the form 
(4.8) Vis) =o (5), BS), Oris ble et 


where P(s) = (f,(s),..-, #,(s)) is an (m+ +4) x k-dimensional matrix, the 
functions #,(s) are given in (4.7), Q,(s, 6) = Qi(s+ 2, 0) is an (n+ k-+41) X n- 
dimensional matrix, Q,(s, b)€C?(R x U*) and H(b) is an m xX m matrix whose 
eigenvalues are «,(0),..., «,(b). Therefore, ep oe of (4.6) can be written 
in the form y(s) = fo(s) u# + P(s) v + Q,(s, b) eS? OD w, where uw, v, w, are con- 
stants, “is a scalar, v is a k-vector, w, 1s an 1-vector. Agen ae the transforma- 
tion y—> (u,v, w,) transforms the system of equations (4.6) into an equivalent 
system du/ds =0, du/ds =0, du,/ds=0. If we let we? OO, then the 
transformation (wu, v, @,) > (u, v, W2) leads to the equivalent system du/ds = 0, 
dvjds=0, dw,/ds = [w(b)|4H(b) w,. Or, returning to the ¢ variable by the 
relation s =qw (b) t+ @, we can sees eRe fis by saying that there exists a matrix 
Q(s, b) = Q(s+a, b) Q(s, 6) = Q(s, b) MOT’? such that the transformation 
of variables 


0x9(s, b) 0 x°(s, b) we ) 


(4.9) y(t) = as CPs Be ue Es 2 


v + Q(s, 6) w 


s=q(b) t+, transforms system (4.5) into an equivalent system 


du dv dw 
(4.10) Obie oe Wye: a = Ho) w. 
Cn 07) 0m (6) One 


Notice that det Ba as apn zap ,Q(s,b)| + 0 for all € U*, where 


s=(b)t+ 9, since all of the above transformations are one-to-one. 


The matrices Q(s, 6) and H(b) may be complex, but it is known from the 
PoincarRé-LyaPunov theory (see, ¢.g., A. M. Lyapunov [14, p. 399] or I. G. 
Makin [16, p. 213 ff.]) that it is always possible to choose both Q(s, 6) and H(6) 
real, O(s + 2, b) = Q(s, b), if none of the characteristic multipliers of (4.6) is 


1A 
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negative. If some of the characteristic multipliers are negative, then Q(s, 8), 
H(b) may still be chosen to be real matrices provided that it is only required 
that Q(s + 22, b) = Q(s, d), z.e., Q is periodic in s of twice the period of (4.6). 
These results are summarized in the following lemma. 


Lemma IV.1. If condition (4.4) is satisfied, if n of the characteristic multipliers 
of (4.5) are +: 41, and if none of the characteristic multipliers of (4.5) are negative, 
then there exist real matrices Q(s,b), H(b), QEC?(R x U*), H(b) ec? (U"); 
Q(s-+ 2, b) = Q(s, b) of dimensions (n+k-+ 1) x n and n xX n, respectively, such 
that the transformation of variables (4.9) is one-to-one and transforms system (4.5) 
into an equivalent system (4.10), where the eigenvalues of H(b) are the nonzero 
characteristic exponents of (4.5). If some of the characteristic multipliers of (4.5) 
are negative, then the matrices Q(s, b), H(b) may still be chosen real, provided that 
Q(s +22, b) = Q(s, b). Furthermore, Q(s, b) 1s defined in such a way that w= 
Q(s, b) &#) together with the functions in (4.4) yield a fundamental system of 
solutions of (4.5). 

For later reference notice that 


(4.11) CNS PY oy (b) = X[2%(s, B)], 


from (4.2) and (4.3), and that 


(4.12) oe  O(8,8) = OE a/b) + Q(s,0) HO), 


by differentiating (4.9) and using (4.10), (4.11) and the fact that the functions 
in (4.4) satisfy (4.5). 


Lemma IV.2. If condition (4.4) is satisfied, then there exists a matrix Q(s, b) 
(given in Lemma IV.1) such that the transformation of variables 


(4.13) x = x°(s,b) + Q(s, d/h, h=(h,...,h,), 


transforms system (4.1) into an equivalent system, 


d 

= =o (b) + Q,(s, b, h) — 2, (s, b, h, t) , 
(4.14) @ = Pi(s,b,h) + Pls, d, 2, ), 

adh 

ay = Hb) h + Hy (s,0,h) + Hy(s,b, ht), 


for all he U} for some o>0. Furthermore, if (0Q/0b) h=[(4Q/00,) h, ..., (@Q/0b,) hj, 


(4.45) A(s,b,h) =|F2 4+ 2p, 2 4 205, of 


then A(s,b,h)€C!(R x U* x U") and 


(21, BH) = Als, ®)|— X(08) — 250°) On +. XU + OM), 
(Q,, Py, Hz) = A(s, b, h) X*(t, x94 Qh), 


(4.16) 
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Proof: Let Q(s, b) be the matrix defined by Lemma IV.1. Then, for a fixed 
b, det [dx/ds, (Ax/As) (Aw/4b) t+ 0x°/0b, Q(s, b)] + 0, s=w(b) t+, for all ¢. 
But this implies that det[dx/ds, 0x°/4b, Q(s, b)] + 0, for all s. But since each 
of the elements of this matrix is periodic in s, we have 


0%%(s, b) @%°(s, db) 
as ? ab 5) Q(s,5)] 0. 


Now, the transformation x—>(s,a,h) given by (4.13) (notice that when 
h=0,s=at-+ o, b=constant) yields, from (4.1), (4.11) and (4.12), the equations 


0x 0d ds 0x Q 
ie r melee ya an dt roa Hh) 
eee Xo foe OX (#") On + X(x9 + Oh) + X*(t, x + Oh). 
Since (4.17) is satisfied, there exists a 9 >0 stich that the determinant of 


ds db dh : n 
i ads Hhis + 0 forallh€ UY, 


(4.17) min, det 


the coefficient matrix of the variables 


b€U*, sER, and the result follows. 


V. Applications to essentially autonomous systems 
Consider two systems of differential equations 


(5.4) oe ak (esie50: 
(5.2) SE = X(x) +X", A gicki=0s 


where x, X, X* are real m-vectors and X, X* satisfy conditions which assure 
the uniqueness of the solutions of (5.1), (5.2). 

Suppose that Pane (5.4) has a family of solutions x = x°[@, (61) t+ q,..., 

w,,(0”) t-+,,, b',...,6"] where b’ is a vector of dimension k; bE U;., where 
Ti is either an open set in R® or consists of a single point, @; is i real Eorcane 
—c<g,<-+ oo, and o,(b’) is a given scalar function, 7=1, 2,...,m. Also, 
suppose there exist real Honibers Piienop «suchsthat weiss. 2, aes Props 
Deer liao (ye Oe te, Om) 2 10r— all. O5C Ure mY Rares ee 
4,2,...,m. This represents a k-parameter family of solutions where k= 


> dim (UY) +m. 
j—1 
To understand a little better the nature of such a family of solutions, con- 
sider the case m = 1. If U, is an open set in R™, then the solution x = x° [w (b!) t+ 
@1,0'] is a (Rk, + 1)-parameter family of periodic solutions of (5.1), and if U 
consists of only a single point, then x = x° corresponds to a one parameter family 
of periodic solutions (the only parameter being the arbitrary phase constant 9,). 
In case m= 2, then x= x° is a family of doubly periodic solutions and by an 
appropriate mapping corresponds to a family of tori. 
In the following, we shall let b= (b},...,6") and U=Ux-::- xU,. If 

b€U is fixed, define C, as the set of points in (x, ¢)-space defined by x =x (Hj, 

sm, 0), 0 <— F< + 00,7 = 1, 2, 00, Mp Ot <a OO, andlet Mt= UC). 


For any subset U* of U, let I* Feed N,(IN*) be a o-neighborhood of m* 
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in (x, t)-space and 2, (M*) = [N, (M*) A (R*-1— )] U M*. Note that in general 
N, (M*) is not a neighborhood in the usual sense, but if M*=M, then XR, (M) = 
N, (M). 

The manifold IM of solutions of (5.1) will be said to be an asymptotically 
stable manifold of (5.2) if given e>0 and U*CU, where U* is a closed subset 
of U, there exist constants o>0, T =O such that every solution x(t) of (5.2) 
with («(T), T)€N,(M*) satisfies the relation (x(t), t) EN (M*) for t2T and, 
in addition d[(x(t), 1), M]—0 as t+ oo, where d(z, A) is the distance from a 
point z€ R"*! to a set ACR". If M is an asymptotically stable manifold of 
(5.2) and for every solution x(T) with («(T), T) €%,(M*) there exists a b)€U, 
which depends upon x(7), such that d[(x(¢), t), C;,| > 0 as t + oo, then we shall 
say Mis asymptotically stable with asymptotic amplitude. If M is an asymptotic- 
ally stable manifold of (5.2) and for every solution x (¢) of (5.2) with (x(T), Te 
N, (M*), there exist by U, pjo€R, 7 =1, 2, ..-, m, such that d[ (x(t), t), (y(), 4) ] 
—> 0 as t > co, where y(t) = x°[a,(bd) t+ Pros «++» Om (00) E+ Pmol, 00, +++» 00, then 
we shall say M is an asymptotically stable manifold of (5.2) with asymptotic ampli- 
tude and phase. In case U consists of a single point, then we shall omit the phrase 
“asymptotic amplitude’. 5 

Notice that the manifold 3% need not be an invariant manifold of (5.2) unless 
X*(t, x) is identically zero. 

To understand the significance of the above definition, let us compare it 
with well known examples. First, suppose that X*(¢, x)= 0 for all ¢, x and 
that 2% consists of a straight line in (x, t)-space, 7.e., equation (5.1) has a constant 
solution x°. The above definition of stability coincides with asymptotic stability 
of the solution x = x® of (5.1) in the sense of LiapuNov. As another example, 
suppose X*(¢, x) =O for all ¢, x and that (5.1) has an isolated periodic solution 
x = x°(t). Since the system is autonomous, x = x°(¢-++q@), where gis an arbitrary 
constant, is also a solution of (5.1). Therefore, in (x, ¢)-space, the family of solu- 
tions x = x°(t+ ) represents a cylinder MM of solutions of (5.1). The concept 
of asymptotic stability of § above is equivalent to the asymptotic stability 
of the orbit of x° in the x-space. The concept of asymptotic stability of I with 
asymptotic phase coincides with the results given in E. A. Coppincton & N.LEv- 
INSON [8, p. 323] for this case. The reason for introducing these more general 
concepts will be clear from the theorems and examples which follow. 


Theorem V.1. Suppose that the system of equations, 
(5.3) = X(x) 


where x, X are veal (n+ k-+ 1)-vectors, XEC2, has a (k + 1)-parameter family of 
veal periodic solutions, 


(5.4) x = x°[w(b)t+ 9, b], DONS rb) SAO 


where b = (by, ..., by), p are real constants, x°(s, b) € C3(s), x°(s, b) € C2(b), w(b) > 0, 
w(b)EC?, BCU, U is either an open set of R* or a single point®. If n of the char- 


* If U consists of a single puint, then the dependence on b will not be mentioned, 
and we take k=0, This remark applies to the remainder of the section. 
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acteristic exponents A,(b), ..., 2,(b), of the linear variational equation of (5.4) have 
negative real parts for all bC U, and 


DENG (0) eB 10) |) 
(5.5) rank a ae k+ 4 


for all bEU, sER, then the manifold IN defined by (5.4) is asymptotically stable 
with asymptotic phase and amplitude. 


Proof: By Lemmas IV.1 and IV.2 and assumption (5.5), there exists a real 
matrix Q(s, b) = Q(s + 22, b) (twice the period of x°) such that in a neighborhood 
of any periodic solution (5.4), say b=), the transformation x = xs, b) + 
Q(s, 6) h transforms system (5.3) into an equivalent system 


ds 

& =0(b) +2(s,6,h), 

db 

“ap = EMS, 8, P), 

O) = Hy(b)h + HE (s, 0, h), 


where the eigenvalues of H(b) are A,(b),...,A,(0), provided that ||b—dp|| is 
sufficiently small. If we let b= b)+a, then 


@ = w, (a) + S,(a,h, 5) 
d 

oa = A,(a,h,s) 

G. = Ha) h + Hy (a,h,s), 


and this system is obviously a special case of the system in Theorem III.1, with 
S (a, h, t) in (3.3) equal to w, (a) and all other functions in (3.3) independent of t. 


Remark V.2. From Theorem 111.1, it is clear that more has been proved than 
stated. In fact, not only does a(t) approach a limit c, but also this limit can be 
made as close to ay as desired, and s(t) —[@,(c) + g’]->0 where q’ is as close 
to s(0) as desired. 


Remark V.3. In case k = 0, condition (5.5) of Theorem V.1 is automatically 
satisfied. The case k = 0, 7.e., an isolated closed curve in phase space, is exactly 
the same as Theorem 2.2 of [5, p. 323]. 


Example V.4. Consider the sit of equations 


e 3 + 8(x) =/(%,9), 
(5.6) 
D + y(x) 9 = h(x, 9) 


where x, y are scalars, g(x), v(x), f(x, v), h(x, vy) € C? for all x, y; O=7(%, 0) = 
h(x, 0) = af (x, 0)/Ay = ah(x, 0)/8y, ¢(—x) =— g(a), x8 (x) >0, x +0, y(x)>0, 
for all x, y. System (5.6) has a two-parameter family of periodic solutions y = 0, 
a" (« (6) t+ 9, b), x°(s + a, b) = x°(s, 6) for all b. In fact, these solutions are 
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defined by the equations 
Oe a tel) =0, y=0. 


We assume that det (0x°/0q, 0x°/0b) +: 0 for all 6 and w(b)€C?. The linear 
variational equations for a fixed b are 

du ! eo’) Fay 

dt? Ox 


a + y(x°) 0 =0, 


(5.8) 


and from the Floquet theory one of the characteristic exponents of (5.8) is 
given by 


A(b rig [70 b)t + ,b)\dt 


where 7T(b) = 2/w(b). Since y(x) > 0 for all x, A(b) <0 for all b. Condition (5.5) 
of Theorem V.1 is satisfied. Therefore, by Theorem V.1, the manifold of solutions 
defined by (5.7) is asymptotically stable with asymptotic phase and amplitude. 


It is clear that this example could be made much more general. In fact, 
it is not necessary to suppose x g(x) > 0, but one need only suppose that system 
(5.7) has a two-parameter family of periodic solutions. 


This example could also be generalized to one of the form 


males at) = 1 (* a 
8 + 9(0, 9 <i 889 


where y, g,/,4€C?; f,A and their first derivatives with respect to y vanish 
for y=0; y(x, x’) > 0 for all x, x’; and the equation 


60 sf ele f= 


(5.9) 


has a two-parameter family of periodic solutions. Conditions under which this 
last requirement is satisfied have been given (see, e.g., [4)]). 


Example V.d. Consider the system of equations (' = d/dt) 


(5.11) Kl Alm axle h(x,’ ,'s) 

Where 4 =(%4,.../%,), A= diag’ (a7; ...,'o7), FEC™ and ae this system 

has a two- ene family of periodic solutions x = x°(wt-+ @, b, €), where 9, 6 

are scalars, wo = o (b, 2 w(b, 0) =0,, 4 (wt+ g, b, 0) =b sin (ot ‘a Q), x (wt+ 

y, 6,0) =0, 7 =2,3,...,. For a fixed b, the linear variational ere is 
yl! 2 OL, oe 8) Af (a 48", , 

(5.12) +Ay= — Pre auaes SUT: 


This equation has two characteristic exponents equal to zero, and if 


(5.43) 26,30, o,+0o,+#0(modo,), 7+, 1 RiD Bow 


ey , 
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then it is clear that one can satisfy condition (5.5) for ¢ sufficiently small. Also, 
if (5.13) is satisfied, it is known [9a] that the real parts of 2 — 2 of the charac- 
teristic exponents 7 3,..., T:, of (5.12) are given by 


T 
TOPE ee eee Ee ; 
(5.14) (T2;—1) aT | ax (x°, x®’, 0) dt + O(e?), 


Voz =T27-1> 7 = 253,...5M, T =2a/o,. 


Consequently, we have the following result: Jf (5.13) is satisfied and 
ia 
(5.45) EN OTM, xo ONO NAL O17 2= 2,5, 5., 8, 
0 


then the manifold defined by the two-parameter family of solutions x = x°[wt + ,b, €] 
of (5.11) ts asymptotically stable with asymptotic phase and amplitude. 

A particular case in which system (5.11) has such a family of periodic solutions 
is the following. Let x, = 4, (x.,..., %,) =v, and then f(x, x’, ec) =f (u, v, v’, v’, €). 
Iff(u, 0, u’, 0, e) =0, then system (5.11) hasa family of solutions x, =) sin (o,¢ + q), 
AO — hae 

Also, one could even assume either that /,(u, 0, — uw’, 0, €) =f, (u, 0, wv’, 0, €) 
or f,(— 4, 0, wv’, 0, €) =— f, (u, 0, uw’, 0, e) and f,(u, 0, wv’, O, €) = 0, 7=2,3,...,% 
(see, ¢.g., [4]). 

An interesting property of the above result is its relationship to a previous 
result on stability. In case (5.11) has an isolated periodic solution x (¢) = x°(wt + 
Y, €), @ = a(e), w(0) = 01, 8 (wt + g, 0) = dosin (o,¢ + ¢), R(wt+ 9, 0) =0,k = 
2,3,...,, then it is known [90, 20] that this periodic solution is asymptotically 


T 

orbitally stable if e¢ f [2f;(x°, x°’, 0)/@x;]dt<0 for 7=1,2,...,”. The above 
0 

result is exactly the same as this one except the index, takes only the values 


25), BM 


Theorem V.6. Consider the system of equations 


(5.16) , x) + X*(t, x), 


where the system dx|dt = X(x) satisfies the conditions of Theorem V.1, and suppose 
that there exist continuous functions h(y), g(t) and a T=0 such that g(t) >0 as 
t—>-+ oo and 


(5.17) ||X*(@ *)l] Se A (ll) 


for all x,t =T. If, in addition to these conditions, 
(5.18) “f g(t) dt< +o, 


then the manifold of solutions (5.4) of (5.3) is an asymptotically stable manifold 
of (5.16) with asymptotic amplitude. Finally, if 


(5.19) J Jam u)dudt<o, 


t 
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then this asymptotically stable manifold of (5.16) has asymptotic amplitude and 
phase. 

Proof: The proof is the same as the proof of Theorem V.1 except that one 
applies Theorem III.4 instead of Theorem TT. 

In case the manifold of solutions (5.4) is only a one parameter family, 7.¢., 
an isolated closed curve in phase space, then the following result holds. 


Theorem V.7. Consider the system of equations (5.16) with X* (t, x) satisfying 
the condition (5.17). Then the cylinder in (x, t)-space defined by x°(wt+q) ts an 
asymptotically stable manifold of (5.16). If, im addition, g(t) satisfies condition 
(5.18), then this asymptotically stable manifold of (5.16) has asymptotic phase. 

Proof: The proof is the same as the proof of Theorem V.1 except that one 
applies Theorem III.7 instead of III.1. 

Let us discuss in more detail the implications of the two previous theorems. 
Consider first the case where system (5.13) has an isolated asymptotically orbitally 
stable periodic solution x°(wt+ ¢). In the (x, t)-space this means that the cyl- 
inder C of solutions x°(wt-+ q) is asymptotically stable and, in addition, any 
solution x(t) of (5.3) with x(0) close enough to the cylinder approaches a parti- 
cular solution on this cylinder as t > oo. The first observation to be made about 
the above two theorems is that the cylinder C defined by the solutions x° (wt-+ q) 
of (5.13) need not contain any solutions of (5.16). Secondly, the solutions of 
(5.16) may not be defined for all f= 0, but only for all ¢=7, where T is suffi- 
ciently large. Third, the first part of the conclusion of Theorem V.7 states only 
that there exists a JT = 0 such that every solution x(t) of (5.16) with x(T) close 
enough to the cylinder C approaches C as t->o, but this solution does not 
necessarily approach any x® as ¢->oo. On the other hand, if condition (5.18) 
is satisfied, then every solution x(t) of (5.16) satisfying the property above 
approaches a particular one of the functions «°(wt+ q) as t>o. 

Now suppose that system (5.3) has a two-parameter family of periodic solu- 
tions, 7.e., the solutions generate a one-parameter family of cylinders C, in the 
(x, t)-space or a three-dimensional manifold $M solutions. The first part of 
Theorem V.6 states that I is asymptotically stable with asymptotic amplitude; 
that is, the solutions “‘close to” 3 approach one of the cylinders C, as t>o, 
but do not necessarily approach any «° as {> oo. On the other hand, if (5.19) 
is satisfied, then the solutions do approach some x° as t-> oo. One interesting 
relationship to notice between Theorems V.6 and V.7 is that the function X*(t, x) 
must approach zero faster in Theorem V.6 than in Theorem V.7 to obtain essen- 
tially the same qualitative behavior near the manifold of solutions. These re- 
marks are illustrated by the following examples. 


Example V.8. Consider the second order equation 


6a i bal $8) olen tn 


where ?, wy € C? with respect to x, dx/dt; there exists a T > 0 such that w is con- 
tinuous in ¢ for = 7; and there exist continuous functions g(t), h(x, y) such that 


(5.21) p(t, % ¥)|<e(d) h(\|x 


Vella es eh 
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for all x, y and g(t)—>0 as t->oo. Also, suppose that the autonomous system 


2 
(5.22 2 + 9(s, 48) =o 


has an isolated periodic solution x°(wt-+ ) whose nonzero characteristic expo- 
nent is less than zero. This implies that the cylinder of solutions of (5.22) gen- 
erated by x°(wit-+ q¢) is asymptotically stable with asymptotic phase. Many 
papers have been written concerning the existence of such a periodic solution 

f (5.22) (see, e.g., [4, chapter 9]). System (5.20) satisfies the conditions of 
Theorem V.7, and, therefore, the cylinder generated by x°(wt-+ q) in the (x, #)- 


space is an asymptotically stable manifold of (5.20). If, in addition, [ g(t) dt 
Te 


<-+ oo, then this cylinder is an asymptotically stable manifold of (5.20) with 
asymptotic phase. 

Now consider ae iene case where 9(%, y) =mMof (x) y+ G(x), Uy >, 
p(t, x, v) = [My — H(t)] f(x) y, where p(t) is continuous for ¢ 20, u(t eel as 


a cay onc (4) <O'for 4 S=\0 and F(%) =F 10) 


is negative for 0<x%<% and positive for x > x for some finite ¥ and F'(x . —>+ oo 
monotonically as || x|| > oo. The equation «’’ + fof (x) x’ +G(«) = 0 has a unique 
solution x°(¢) which has a negative characteristic exponent (N. Levinson & 
O. K. Situ, [13]). Consequently the manifold generated by x°(wt + @) in (x, #)- 
space is an asymptotically stable manifold of the equation x’ + w(t) f(x) x’ + 
G(x) =0. This result has been previously obtained by L. Markus [/7, p. 27]. 
The result of MARKUs is actually more general than the one here since he shows 
the solution is globally asymptotically stable. This can be done in two-dimensional 
phase space, but as we shall see in Example V.17, a result similar to the one 
given here is valid in higher dimensions. Also, we get the further result here that 


if f | w(t) —wo| dt <-+ ce for some T= 0, then the above manifold is an asymp- 
T 
totically stable manifold with asymptotic phase. 


Example V.9. Consider the system of equations 


T* + g(x) =H(x,9) + pl) Flx,y), 


(5.23) 
Ge +(e y = h(x, y) + p(t) A(x, y) 


paras J gta F, HEC, y, pec, 


(t)| dt < +, 
i | w(t)| dt <<-+ co and suppose that the limiting equation 


(5.24) igi dC? 
& + y(x) y = h(x, 9) 


satisfies the same conditions as system (5.6). Then the asymptotically stable 
manifold of solutions of (5.24) is an asymptotically stable manifold of (5.23) 
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with asymptotic amplitude. This follows immediately from Theorem V.6. If, 
in addition f [| p(u)| du dt <oo, f f |p(u)| du dt < oe, then, from Theorem V.6, 
ey i is 


this asymptotically stable manifold has asymptotic amplitude and phase. Ob- 
viously, this example could be generalized as in Example V.4. 


Theorem V.10. Consider the system of equations 


(5.25) 


where x,y are vectors; P,X,X,, Y;CC? with respect to x,y; and X,, Y, aré 
continuous and bounded in t for all x, y in some bounded region, and X,(t, x, 0) = 9, 
Y(t, x, 0) =0, OY, (t, x, 0)/Oy = 0. Furthermore, suppose that the system x' = X(x) 
satisfies all of the properties of system (5.3) of Theorem V.1; in particular, there 
is a manifold, M, of solutions of this equation, defined by x = x°[w(b) t+ Q, bd], 
which is asymptotically stable with asymptotic amplitude and phase. If, in addition, 
the characteristic exponents of the periodic system 


(5.26) << = P[x°(w(b)t + 9,5)] y 


have negative real parts for every b, then the manifold IN, defined by x= 
[w(b)t +, 6], y=0 ts an asymptotically stable manifold of solutions of (5.25) 
with asymptotic amplitude and phase. 


Proof: The functions x = x°[w(b)t-+ q, 6], y = Oareasolution of the equations, 

dx|/dt = X(x), dy/dt = P(x) y, and the linear variational equations are given by 
GH OA (A) | ay. 

(5.27) Ee RCO pe Pla?}"y., 
where the coefficients are periodic of period a in the variable s=w(b) t+ 9. 
Let the real matrix Q(s, ) =Q(s +22, 6) associated with the first equation in 
(5.27) be defined by Lemma IV.1. In the same way, let R(s, b) be the real non- 
singular matrix associated with the second equation in (5.27) defined by 
Lemma IV.1. 


Observe that 


*" «o (b) = X [x%(s, d)], 
(5.28) SALE 9 (5,8) = 209) 4 () + O(s,8) HO), 
* w() + R(s, b) L(b) = P[x9(s, b)] R(s, 6), 


where the eigenvalues of the matrices H(b) and L(b) are the nonzero character- 
istic exponents of the first and second equations (5.27), respectively. These 
matrices are determined by Lemma IV.1. Now, if we apply the transformation 
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x = x°(s,b) + Q(s,b)h, y= R(s, 6) z to system (5.25) and use (5.28), we obtain 
ax | 80 ,\(ds __ ax  0Q,\db | dh 
ee + 3s \ (ai «) ! lee ey i) di + Ols.6) er 
= X(x°+ Qh) — X(x9) — 2) On + x, (6, 294 Qh, Re), 
oR (ds oR . db dz 
Seta; 8) tap PRG Li 


= P(x®+ Qh) Rz—P(x°) Rz + Y, (i, x°+ Qh, Rz). 


Now, by assumption (5.5) concerning the system dx/dt= X(x) and the 
fact that R(s, bd) is nonsingular, the determinant of the coefficient matrix of the 
variables ds/dt—(b), db/dt, dh[dt— Hh, and dz/dt— Lz in the preceding 
system is different from zero for h= z= 0. Therefore, for h, z sufficiently small, 
one can solve for these variables to obtain 


d 
= — w(b) = S,(b, h, 2, 8,2), 
db 
ae = A,(b,h; 2):s3t)s 
adh 
(5.29) 7h TO) P= 17 (b, WS; 1), 
& —1()2 =Z,(0,h,2,,2) +Rik, S,(b, h, 2, 8,4) + 


+ RAS 7A, h, 2, 8, t,) 


where S,, Ay, H, are O(\\A|+||2])) while Z, is O[flz||(\\*l+llzl)] as [ll 
||z|| +0. It is now clear that this is a special case of the system considered 
in Theorem (3.1) since the functions on the right-hand sides of (5.29) have con- 
tinuous first derivatives with respect to s, b, h, z and are continuous and bounded 
in ¢t. The remainder of the proof follows immediately by applying Theorem III.1. 


Remark V.11. It is not at all obvious when one can be assured that the 
characteristic exponents of system (5.26) have negative real parts for every Db. 
One would hope to give conditions for this which depend only upon the matrix 
P(x) in (5.25). Using a result of A. WINTNER [23] (see also R. Conti [6]), one 
can prove the following result: I} the eigenvalues of the matrix B(x) = 27[P(x) + 
P*(x)] (P? is the transpose of P) have negative real parts for every x in a neigh- 
borhood of the manifold IN, then the characteristic exponents of (5.26) have negative 
real parts. In fact, if A (t,b) is the maximal eigenvalue of the matrix B[x°(w(b)t+ 
gy, b)|, then the result of WINTNER implies that every solution x(t) of (5.26) 


t 
satisfies the property ||«()||? < || *o|?exp 2 f A(z, 6) dt. But A (é, b) is negative 
0 


for ali ¢, b and periodic in ¢, and the result follows. In particular, the above result 
implies that the characteristic exponents of (5.26) have negative real parts if 
P(x) is a symmetric matrix whose eigenvalues are negative for all x. 

Remark V.12. Theorem V.10 can obviously be generalized along the same 
lines as Theorems V.6 and V.7 but will not be discussed here since the details 
can be so easily supplied. 
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Example V.13. Consider the system ('=4d/dt) 
x + g(x) =f (t, %¥), 
y ty(ex)y =h(t %, ¥) 


where x is a scalar, y is a scalar, g,y,/, hE C? with respect to x, y, f, 4 are con- 
tinuous and bounded in ¢ for all x, y in some bounded region; Tne), = Oe 
h(t, x, 0) =0, Oh(t, x, 0)/2y = 0 for all ¢, x and the equation, x’ + g(x) =0, has 
a family of periodic solutions *=x«°[w(b)t+g, 6], where 9, b are scalars, 
det [x/p, 0x%/0b] += 0, and w(b) CC. If y(x, x’) >0 for all #2 Jor, ii arie 
a vector and y(x, x’) is a matrix which satisfies the properties of Remark V.11], 
then (5.30) is a special case of system (5.25) if we put %1—=~, %2 =x". [heretore, 
for f, h bounded functions of #, the manifold of solutions defined by x = x°[P (b)t + 
gy, b], y=0 is an asymptotically stable manifold of solutions of (5.30) with 
asymptotic amplitude and phase. This is a generalization of Example V.4 in 
that the variable ¢ is introduced in the functions 7, h and df(t, x, 0)/¢y may be 
== 0: 

Example V.14. As an illustration of the case where y(x, x’) in (5.30) may be 
a matrix, consider the system (’=d/dt) ; 


(5.30) 


x + g(x) = f(x, x, 9,9’) 
yi + u, (x, x") Vj + 0, (x, x) ¥, = h(x, Y) Wi oy i ae ee 


where y= (¥1,---5 Vn)> Ff, Uj, 0j;,4;EC*, F(x, x’, 0,0) =0, A;(x, 0) =0 for each 
7=1,2,...,”, and suppose that the equation x’’+g(x)=0O hasa family of 
periodic solutions x = x°[w(b) t+ g, b], where 0, are real scalars, det [0 x°/dq, 
dx i¢b) += 0,,and w(6)€C*. The manifold += x, y = yi== 0, is assolutionser 
(5.31), and (5.31) will be a special case of system (5.25) of Theorem V.10 if we 
can determine conditions under which the solutions of the linear equations, 


(5.34) 


(5.32) yj +; (x, %°") yy Os COE Vey aceasta Aue ence 


approach zero as ¢—> oo. Since this is a linear equation with periodic coefficients, 
the approach will necessarily be exponential. Many papers have been devoted 
to giving conditions which imply that the origin is an asymptotically stable 
solution of (5.32) (see, e.g., [4, p. 63], [19, p. 210]). In particular, if 

(5.33) vj (x9, x9) > 0, 204; (x®, 29”) v, (29, x0”) + “7 (40, x0) 0, 7 = 1,2... 


dt “ 


md, , 


for every b, where x° = x®[w(b) t+, b], then y¥;=0,7 = 1, 2,..., mis an asymp- 
totically stable solution of (5.32) [12]. Therefore, the manifold x = x®, y=0 
of solutions of (5.31) is asymptotically stable with asymptotic amplitude and 
phase if (5.33) is satisfied. Condition (5.33) is certainly satisfied if v; = constant 
> Oand u,;(%, x’) >0,7=1, 2,..., ”, for all x. This last result has been obtained 
by E. W. THomrson [20] for the case where g(x) in (5.31) is linear in x. 

So far, the theorems of Section 3 have not been used in all their generality ; 
namely, the variable s in the previous applications has always been a scalar. 


The following two theorems are given to illustrate how situations may arise 
where s is a vector. 
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Theorem V.15. Consider the system of equations 
(5.34) 


where x,y are vectors; X,Y, X*, Y*EC?; each of the systems AX dt—X (x), 
dy/dt=Y(y) has a manifold of solutions defined by x= x®[w(b) t+, b], y= 
y[v(c) t+, c], respectively, where p is a constant scalar, b,c are constant vectors 
and each of these manifolds satisfies the conditions of Theorem V.1. Also, suppose 
there exist continuous functions g(t), hy (u,v), he(u, v) and a T = 0 such that g(t) >0 
ast >+ cand 


(5-35) |X*@&x, | Se@ (lal lyll), WY*G*% Nl Se@ Ao((la lb yl), 


for allt=T and x, y in some neighborhood of the manifold MN, defined by x =x, 
Veil} 


(5.36) Sel dt< +o, 


then IN 1s an asymptotically stable manifold of (5.33) with asymptotic amplitude. 
If, in addition, 
(5.37) f fe(u)dudt<o, 

an 


t 
then M also has asymptotic phase. 

In case x°, y® are isolated periodic solutions, then M can be mapped onto a torus 
im (x, )-space. In this case, this torus 1s an asymptotically stable manifold of 
(5.34) even af (5.36) ts not satisfied. If (5.36) 7s satisfied, then this torus also has 
asymptotic phase. 

Proof: From Lemma IV.1, we may define matrices Q,(s,, 0), Qe (Sg, 0) asso- 
ciated with the variational equations a= Xiao (s,, 6) ]/e x} x, yi= OV [iy? (Se, €)]/ 
dy} y, 8, =o(b)t+ , s,=v(c)i+g, respectively. Then the transformation 
x% = x9(s,, b) + Q1(s,, 6) yy, y= y9 (So, ¢) + Qo(Se, ¢) Ag transforms system (5.34) 
into an equivalent system 


ae = (a) +2, (a, h, s) + Q, (a, h, S, t), 
ee = Pans) + P,(a,h,s,t), 
dh 


a = H(a)h+ H,(a,h,s) + H,(a,h,s,t), 

where s = (s,, 53), @= (b, ¢), = (iy, hy), 10(a) =[e(b), »(c)], H(a) = diag [H¥ (b), 
H*¥(c)], where the eigenvalues of Hj(b), H¥(c) have negative real parts for every 
b,c. The functions 2,, P,, H, are O(| h||?) as ||h|| > 0 and 2,, P,, Hy, satisfy 
the conditions of the functions S3;, A,, H, in system (3.14) of Theorem III.4. 
Since by a linear transformation, a@ may be taken in a neighborhood of zero, 
the above theorem for the case where the periodic solutions are not isolated 
follows from Theorem III.4. If the periodic solutions are isolated, then the result 
follows immediately from Theorem III.7 and III.9. 
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Remark V.16. The above theorem may obviously be generalized to systems 
where more than two vector variables have the same behavior as those in (5.34). 


Example V.17. Consider the system of differential equations ("= d/dt) 

a Th (HK OO Oe 9, 9) 

y +e. y') =n G(% #97) 

where /, g, F, GEC?, w,7€C! and each of the equations x + f(x, x’) =0, y+ 
g(y, y')=0 has an isolated periodic solution x= *°(wt+@), y= y° (vt +9), 
respectively, and also each of these periodic solutions has a negative charac- 
teristic exponent. If u(t) +0, y(t) +0 as toe, then from Theorem V.15, the 
torus defined by x= x°(wt+q), y=y(vt+¢) is an asymptotically stable 


manifold of (5.38). If f|m(é|dt<-+ 0, f|n(t)|dt <-+ 0, then this torus is 
0 0 


(5.38) 


asymptotically stable with asymptotic phase. This example generalizes the 
result in Example V.8 and in particular the result of L. MARKus mentioned 
there. 

Theorem V.18. Consider a system of the form 


= X(x) + X*(x, 9,2), 
(5.39) = = Y(y) + Y*(%,¥,2), 
42 _ Plx,y)2 + 2(%,9,2) 


where'x, y, z are vectors; X,Y, P)X*,Y*) and Ze C*; 0 =" X* (x, ¥, 0) = Y*(2,9,0) 
= Z(x, y, 0) = 0Z(x, y, 0)/0z; each of the systems dx/dt= X(x), dy/dt=Y(y) 
has a manifold of solutions defined by x = x®[w(b) t+ q, bj, y= y® [v(c) t+, c], 
respectively, where py 1s a constant, b,c are constant vectors and each of these mant- 
folds satisfies the conditions of Theorem V.1. If there exists a function o(b, c)>0 
and continuous for all b,c such that the fundamental system Z,(t), Z,(0) =I, of 
the system of equations 


(5.40) a = P(x, y°] z, 
satisfies a relation of the form 
(5.41) |Z. Zt (w)|| < K e721") ty S0, K>o, 


then the manifold of solutions of (5.39) defined by x = x°[w(b) t+, b], y= y° [v(c)t+ 
p, ¢], 2=0 1s asymptotically stable with asymptotic phase and amplitude. 
Proof: As in the proof of Theorem V.15, the transformation x = x° (sy, 6) + 


Q1(S1, 4) 4, ¥ = Y%(S2, ¢) + Qo(sz, c),2—= 2 transforms system (5.39) into an equi- 
valent system Ae 
dt Ha) + Q\(a, h, s) + Q3(a, h, 2,3), 


d 

_ —- Bia; h, s) +P, (ah, 2s) 

dh 

Gp = A (a)h + A, (a,h,s) + He(a,h,z,s), 
We 


ae ~ PL (%,5), ¥°(se,0)]2 +Z(a,h, 2,8), 
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where s = (sy, Sy), @=(0, c), h = (hy, hg), w(a) = [w(b), v(c)], H(a) = diag [H¥ (6), 
H¥ (c)| where the eigenvalues of H# (b), H}(c) have negative real parts for every 
b,c. The functions 2,, B, H, are O(|h|"); 2,, P,, H, are O(|lz||) and Z is 
O(\|z|| (lzll+[lAll)) as ||2|| +0, ||z|| 0. Also by a linear transformation, we 
may assume a to be in a neighborhood of zero. The conclusion of the theorem 
now follows from Theorem III.4. 


Remark V.19. It is generally a very difficult problem to decide when condition 
(5.41) can be satisfied, since P[x®, y°] is an almost periodic matrix. In case 
P(x, y) were only a function of one of the variables x or y, then the problem 
would be completely solved by the FLoguer theory. However, in the general 
case, only fragmentary results have been obtained (see, eg, f%, 134): 


Remark V.20. Generalizations of systems (5.39) along the iines of the pre- 
ceding theorem could be given. Also, one could suppose there were more than 
two vector variables which had the same behavior as %, y in (5.39). 


This research was partially supported by the United States Air Force through 
the Air Force Office of Scientific Research of the Air Research and Development 
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part is permitted for any purpose of the United States Government. 


References 


[1] ArEens, R. F.: A topology for spaces of transformations. Ann. Math. 47, 
480—495 (1946). 
[2] Bettman, R.: Stability Theory of Differential Equations. New York-Toronto- 
London: McGraw-Hill 1953. 
[3] BoursBakt, N.: Espaces Vectoriels Topologiques, livre 5. Paris: Hermann 1953. 
[4] Cesari, L.: Asymptotic Behavior and Stability Problems in Ordinary Differen- 
tial Equations. Ergebn. d. Math., N. F. H.16. Springer 1959. 
CoppincTon, E. A., & N. Levinson: Theory of Ordinary Differential Equations. 
New York: McGraw-Hill 1955. 

[6] Conti, R.: Limitazioni ,,in ampiezza“‘ delle soluzioni di un sistema di equazioni 
differenziali e applicazioni. Boll. Un. Mat. Ital. (3) 11, 344—349 (1956). 

[7] DitrBErRTO, S.: On systems of ordinary differential equations. Contributions 
to the Theory of Nonlinear Oscillations, vol. 1, pp. 1—39, 1950. 

[8] DyuMan, E. I.: a) On a reduction principle (Russian). Izvestiia Akad. Nauk, 
Kazakh. SSSR. 4, 73—84 (1950). — b) Some stability theorems. Izvestiia 
Akad. Nauk, Kazakh. SSSR. 4, 85—97 (1950). 

[9] Hatz, J. K.: a) On the Behavior of the solutions of linear differential systems 
near resonance points. Contributions to the Theory of Nonlinear Oscilla- 
tions, vol. 5, 1960. — b) On the stability of periodic solutions of weakly non- 
linear periodic and autonomous differential systems. Contributions to the 
Theory of Nonlinear Oscillations, vol. 5, 1960, 

[10] Kettry, J. L.: General Topology. New York: Van Nostrand 1955. 

[11] Lerscuetz, S.: Differential Equations, Geometric Theory. New York: Inter- 
science 1957. 

[12] Lronov, M. J.: On quasi-harmonic oscillations. Prikl. Mat. Meh. 10, 575—580 
1946). 

[13] enna N., & O. K. Smiru: A general equation for relaxation oscillations. 
Duke Math. J. 9, 382—403 (1942). 

[14] Lrapunov, A.: Probleme général de la stabilité du mouvement. Ann. Math. 
Studies 17 (1949). 

Arch. Rational Mech. Anal., Vol. 6 42 


[5 


Le) 


170 


[15] 


[16] 
(17) 
[18] 


[19] 


[20] 


[21] 
[22] 


[23] 


J. K. Hare & A. P. Sroxes: Differential Systems 


Litto, J. C.: a) Linear Differential equations with almost periodic coefficients. 
Amer. J. Math. 81, 37—45 (1959). — b) Continuous matrices and the stability 
theory of differential systems. Univ. of Kansas, Department of Mathematics, 
June, 1959. 

Mackin, I. G.: Theory of Stability of Motion [in Russian]. Moscow 1952. 
Translated by Atomic Energy Commission, AEC-TR-3352. 

Markus, L.: Asymptotically autonomous differential systems. Contributions 
to the Theory of Nonlinear Oscillations, vol. 3, pp. 17—29, 1956. 

PEerRsipskul, K. P.: Some critical cases of denumerable systems [in Russian]. 
Izvestiia Akad. Nauk, Kazakh. SSSR. 5, 3—24 (1951). 

STARZINSKU, V.M.: Survey of works on conditions of stability of the trivial 
solution of a system of linear differential equations with periodic coefficients 
fin Russian]. Prikl. Mat. Meh. 18, 469—510 (1954). Translations Amer. 
Math. Soc. (2) 1, 189—238 (1955). 

Tuompson, E. W.: Asymptotic behavior and stability for nonlinear differential 
systems. Ph. D. Thesis, Purdue University, June, 1959. 

Tycnonov, A.: Ein Fixpunktsatz. Math. Ann. 111, 767—776 (1935). 

URABE, M.: Geometric study of nonlinear autonomous oscillations. Funk. Ekv. 
1, 1—83 (1958). 

WIntTNER, A.: Asymptotic integration constants. Amer. J. Math. 68, 553—559 
(1946). 


RIAS 
7212 Bellona Avenue 
Baltimore, Maryland 


(Received May 4, 1960) 


On the Singularities 
of Generalized Axially Symmetric Potentials 


ROBERT P. GILBERT 


Communicated by R. FINN 


I. Introduction 


Recently G. Szec6 [1], Z. NEHarrI [2], and A. ErDELy1 [3] have obtained 
some interesting results connecting the singularities of axially-symmetric po- 
tentials with those of analytic functions. It has been shown by the author 
[4], [5] that a similar connection also exists between the singularities of a three- 
dimensional potential and a function of two complex variables. This was done 
by considering the Whittaker-Bergman operator [6], [7]* B3(f, &, Xo) which 
transforms functions of two complex variables /(é¢, w) into ree! dimensional 
harmonic functions, and by using a modified form of an idea previously employed 
by HADAMARD in the proof of his theorem on the multiplication of singularities 
[2], [8]. In this paper we shall use these means to investigate the singularities 
of solutions to the partial differential equation 

2 
Lwl= Se +5 to ae = (1) 
(where k=O), referred to by WEINSTEIN [9] and ERDELYI [3] as the differential 
equation of generalized axially symmetric potentials. This equation is arrived 


2 
at by considering those solutions of the m-dimensional Laplace equation, ay +t... 


Ox a 
- f= =0, which depend solely on the variables z= %,, pe (ees js, 


Solutions to equation (1) may be constructed by a modified form of the 
Whittaker-Bergman operator. However, in this case the operator transforms 
functions of a single complex variable f(o) into axially symmetric potentials 


u(r, ?) (P=2+¢, and #=cos? =) that is 


u(r,8) =An(f, LX) = f Ho) (w— 2) (2) 
aie | 
La 
: H(X)=Bylf, 2X), Bel % X= ats fey, 
AE 
t= [—t inn te+ tin) S, 


|X,>—X|<e, X=(*, y, 2), Xp=(%o, Vo. 20), Where Y is a differentiable arc in the 
u-plane, and ¢ >0 is sufficiently small. 
Dre 
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where o=2-+ 5 o(u+ “IX —X1<e, X se (Ks, Hay hap My ee a ee ee xa)a 


é>0 is sufficiently small, and the integration path is the upper semi-circular 
arc*, Now, if f(o) has a Taylor expansion convergent in some neighborhood 


of o=0, that is 
f(o) = 24,0, (3) 


then the operator A,,(f, &, Xo) defines an axially symmetric potential 


“\ Te : 2 co ! os 
u(r, 8) = (4%) : (/2)) yo ies 7" CH? (cos 8) (4) 


for r sufficient!y small. This may be seen by considering the classical identity [70] 


7" Ch(cos®) = 2" FUT H) [ [z+ igcosg]" (sing) dp. (5) 
0 


Occasionally it is convenient to continue the arguments of w(7, #) to complex 
values, for instance if we introduce the particular continuation 
r= + (22+ 0%), 
z 
i 
which reduces to €=cos# for real z,0, we may obtain an inverse operator. 


Lemma. Let u(r, 0) be an axially symmetric potential regular at the origin, 
that 1s ; 


u(r, 8) = [ f(o)(u— 4)", 


uU 


where 


fo) => ayo"; 


1 _ gn) 234 jie 
Ko) = ghee f Wed) — : ale 9 
zi (ee 


where W(r, cos?) =u(r, 0), and C is the real axis. 


Proof. Let us define 


a 


(rl) = ae te enen(spere: a 
2 2 2 


* Our integrand is an analytic function of @; hence our integral is a Cauchy 
integral, and we need not restrict our path of integration to the upper semi-circular arc. 
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it may then be readily seen from the orthogonality relation for the Gegenbauer 
polynomials 

k+1\ 72 
ri ea 
2) | Dn m (8) 


2 D(k+n) 
(2n+~k) n! 


T(k) 


+1 
fa —2)2 °C CH (6 aE = 
al 4 
(where C is the real axis) that 
fo) = 3,0" = f K(2.8)(a 3 tem ci e)ae 


zi 2 “ 
(where a=(4i)* Ir ( Ey /4}. We note that K ( < j é) may be formally summed 
to give 


72 
K(2 é) = ye—th | T(k) a—ey2 Ja —k @ lt Deer }} 
ge ag k k+1 
Te nes | 43 
; ee 2 
ied rh) [ a—e-3 (1 ma (9) 
et aw Da. SE E21 Lee TP 
r(Z)r( >) (1-22 + 5)? —1 
providing *| <1; this follows from the classical identity 


~ 

Dec (é=A1—2€t+#)-, for |t| <4. 

nD 
The constant coefficient may be reduced by using LEGENDRE’s duplication 
formula, 


L( pre e 3 
siete ee aE 


one then obtains for the kernel 

2 
ke (1—#)" a 2 (1— : #) 
(24)? a ) 


K ( = ) -y oe 
y | Wy 8 Ray 

nae + a) 

We now note that & ( <, é } is analytic in é, and also, since u(r, #) is harmonic 


W(r, &) is analytic in € except for the points =+1. We may conclude from 
this that the integral involved in our representation for /(a) is a Cauchy-integral 
and that we need not restrict our integration path to the real axis. 


II. Singularities of Generalized Axially Symmetric Potentials 
In a previous paper [4], the author proved the following theorems concerning 
the singularities of three-dimensional harmonic functions. 
Theorem 1. If Z*=E{S(X|u) =O}, is the singularity manifold of 7 i (t,), 
then 
ne! du 
H(X) - Qni fil u) eae 
g 


174 RoBERT P. GILBERT: 


(where L is the unit circle) is regular at X=(x, y,2) providing this point does 
7) 
not lie simultaneously on the two surfaces S(X|u)=0, and ae S(X|u) = 
: phe aes Bs 4 
Theorem 2. Let H~(X) be a harmonic function regular at infinity ~ (x1 = + . 
and let V(r, &, @) be the function obtained from H™(X) by replacing x, y, 2 by 


= (0+ 0%) fi—#, 
(oc on 
Beetle. 


Furthermore, let Z = E> {v= @(E, o)} be the singularity manifold of V(r, €, @) = 
H®(X), X€¢% (where {3 is a region in the ideal, Euclidean three-space, that is 
the three-dimensional space of complex points), then the functions of two 
complex variables g(s,) which may be generated by the integral operator 


~ v(s+t) 
[AEB ve.80 


+1 


“2 dé 
0 


r 
(the integration path in the &-plane is the real axis connecting —1 to +1, 
and the path in the g-plane is the unit circle) is regular at (s, w) providing (s, ) 
does not lie on the “‘envelope”’ of the three-parameter family 


p(s,u|r.é,.0) =t—s=r]é 4 -yi—a(= | al Se¢ 


2 Q U 


subject to the auxillary condition r= D(é, 0). 


We now state the following two theorems concerning the singularities of 
generalized axially symmetric potentials without proof, since their proofs are 
essentially identical in structure to those of Theorems1 and 2. We mention 
in passing, however, that the proof of these theorems depends solely on the 
concept of attempting to enlarge the original domain of definition for a function 
given in terms of an integral representation by continuously deforming the 
integration path, subject to the condition that in this process of deformation 
the integration path at no time crosses a singularity of the integrand. 


Theorem 3. Let Z?=E{S(r, | u)=0} be the singularity manifold of —f(o o) X 


k— 
( a a) i ; then the axially-symmetric lacus 


. k-1 
u(r, 0) =A, (f, L, Xo) = [ro to -4) or 
sis a u = 0) 
is regular at (rv, 0) providing (r, 9) does not lie simultaneously on the two surfaces 


a) 
S(7,3|u)=0,) and ay 3" Olu) = 


a £0 ( foe) +n 
Je X) => fen — “ andg (t, u) has the form g(¢,u)= 3) Di ay,t-"—1u™. 


n=0 m=—0 
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Theorem 4. Let u(r, 8) be a generalized axially symmetric potential, and let 


Wir, &) be a function obtained from u(r, 3) by replacing 3 by cost (=). Further- 
more, let Z=E*(r=y(€)} be the singularity manifold of W(r, &); then the function 


1 ga)" ae 
eat pig ait (7 ae 
a (1-28 o +=)’ 


(C is the real axis between —1 and +1) is regular at o providing this point does 
not le on the ‘‘envelope’”’ of the two parameter family w(o|r, £) =r —20é+07?=0, 
subject to the auxillary condition r=y(é). This may be restated by saying o must 
not lie simultaneously on the two curves 


yr (o|€) = 72(€) — 20€ 7 (€) +0? =0, 


To illustrate the use of Theorems 3 and 4, we consider the case where c= 


is a singular point of = f(o) (u SS a. Here the singularities of the integrand 
occur foro =z+ 4 0 [u + =| =a,and we may represent the singularity manifold as 


a 


S(z, |) He) 5 Ol) 0. 


Eliminating « between S=O and = =0, we obtain the locus 0?+ (z—«)?=0. 
If u(r, #) is regular in the unit sphere but has singularities on r=1, they are 
located at z= 5 («+ =} When k=1, C4()=P,(é), and this is equivalent to 


NEHARI’S result concerning Legendre expansions [2]. 

We next consider the possible singularities of an /(a) generated from a W(r, &) 
with the singularities @?+ (a —z)*=7?—2aér+«a?=0. This is done by obtain- 
ing the “envelope” of the two parameter family y(o|7, &) =r? —20&+07=0 

iis a 


subject to = S (= ait =} Eliminating 7 between y, =r? —o (= + a} = O-==()) 


th 


and vs =2r(1 =| =0 yields o=a. We shall now attempt to formulate a 
YA a 


necessary and sufficient criteria for the occurence of singularities from the above 
results. 

The points =+1, and w=+1, because they are the terminal points for 
the paths of integration, have a special status; indeed, a branch of W(r, €) may 
have a singularity at £=-++1, which does not correspond to f(a) at o=+r7. 
f(a) on the other hand will not be singular at w=+1 for all X, hence all the 
singularities of f(a) will correspond to singularities of W(r, €). We realize, there- 
fore, that W(r,&) is singular if and only if the point (7, €) lies on the locus 
v2 —2arE+o2=0 (here o=« is a singularity, of f(c)), the one exception being 
points on the z-axis (7, +1), which may be singular points of W(r, ) without 
corresponding to singularities of f (0). 


ae 
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